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Ogo 5 Magnetopause Structure and Classical Reconnection 

B. U. 6. SONNERUP 

Dartmouth College, Hanover, New Hampshire 03755 

B. G. Ledley 

Goddard Space Flight Center, Greenbelt, Maryland 20771 

Thy; observations made during one unusual crossing of the magnetopause by the satellite Ogo 5 are 
compared with the classical reconnection model developed by Levy, Petschek, and Siscoe. The magnetic 
field observations appear to be generally consistent with this MHD model, although allowance must be 
made for the fact that the estimated magnetopause thickness was no more than 3.5 ion gyrodiameters. The 
nature of the finite gyroradius eflects in such thin structures is discussed. 


Introduction 

More than a decade ago, Levy et a/. [1964] proposed a 
magnetohydrodynamic (MHD) model for steady state mag- 
netic field reconnection at the magnetopause which has served 
as the focus for most subsequent attempts to observationally 
verify or deny the occurrence of this process near the subsolar 
point. A great deal of indirect evidence has been found which 
is compatible with the assumption of dayside reconnection 
Burch, 1974}. However, on the whole, direct evidence 
bearing on the occurrence of dayside reconnection has been 
more negative than positive. Minimum variance determina- 
tions of the magnetic field component perpendicular to the 
magnetopause [Sonnerup, 1976} show that a measurable com- 
ponent is present in a significant number of magnetopause 
crossings, in agreement with the reconnection model. How- 
ever, other aspects of the magnetic field structure of the mag- 
netopause predicted by the MHD model are seen only rarely. 
In general, magnetopause conditions appear to be highly time 
and/or space variable making the comparison with a steady 
state theoretical model difficult. The most negative piece of 
evidence so far has been the absence of energized protons in 
the plasma boundary layer observed just inside the magneto- 
pause [Haerendel et al., 1978]. As first emphasized by Heikkila 
[1975], energization of the plasma at the magnetopause is an 
unavoidable consequence of all reconnection models. 

It is the purpose of this paper to present the magnetic field 
record of a single magnetopause crossing which agrees in 
essential parts with the predictions of the MHD model by 
Levy et al. and where deviations from that model may be 
explained in terms of the finite gyroradii of protons. We find 
the agreement between theory and observation to be suf- 
ficiently remarkable to warrant a detailed discussion of this 
single crossing even though it is an uncommon case. However, 
it is to be emphasized at the outset that we do not consider this 
crossing to provide incontrovertible evidence for the occur- 
rence of reconnection, and certainly not for the importance of 
that process, at the magnetopause, Interconnection of field 
lines on the two sides of the magnetopause in this crossing is 
established by the presence of a nonzero normal magnetic 
field. But reconnection requires in addition that an electric 
field E|j tangential to the magnetopause be present. To prove 
the occurrence of reconnection, either Eh must be measured 
directly or else the particles energized i n this field m ust be 
detected. Afozer et aL [1978J claim to have detected a sub- 

Copy right © 1979 by the American Geophysical Union. 


stantial nonzero value of E||, but it remains to be established 
whether these results are valid and reproducible. 

The Levy-Petschek-Siscoe Model 

A brief review of the main features of the steady state MHD 
reconnection model is in order. A schematic of the flow and 
field geometry is shown in Figure 1. Except in the immediate 
vicinity of the reconnection region, at the center of the 
figure, the magnetopause consists of a large-amplitude inter- 
mediate wave front (Alfven wave; rotational discontinuity) 
across which the component of the magnetic field tangential to 
the front changes direction abruptly. When the magnetosheath 
field is antiparallel to the earth’s field, as assumed in Figure 1, 
the direction change is by 180®, but in principle the angle 
change is arbitrary. Simple MHD theory of a one-dimensional 
steady wave-front structure indicates that the tangential field 
component rotates with constant magnitude from the magne- 
tosheath to the magnetosphere direction. Since the field nor- 
mal to the front remains constant, the total magnetic field 
magnitude is preserved during the field rotation. 

MHD theory makes no prediction concerning the sense of 
rotation of the tangential field. However, kinetic models of the 
rotational discontinuity, to be discussed at a later point in this 
paper, indicate that structures of thickness comparable with an 
ion gyroradius have the electron whistler polarization. But the 
normal magnetic field component, along which the wave prop- 
agates, has opposite sense north and south of the reconnection 
region, pointing toward the earth in the former case, away 
from it in the latter (see Figure 1). Thus the electron-whistler 
polarization implies that when observed by a satellite crossing 
from the magnetosphere into the magnetosheath, say, the ac-* 
tual sense of rotation of the tangential field vector in a thin 
magnetopause structure would be clockwise north and coun- 
terclockwise south of the reconnection region, assuming the 
observer is facing the sun. 

The magnetosheath plasma flows across the magnetopause 
at a speed equal to the Alfven speed based on the normal field 
component, preserving its density and temperature. As the 
plasma crosses the magnetopause, the intense magnetic forces, 
j x B, increase its velocity component tangential to that surface 
by an amount approximately equal to twice the Alfven speed, 
based on the tangential field component. Thus the plasma just 
inside the magnetopause flows nearly tangential to the magne- 
topause in two high-speed jets directed away from the recon- 
nection region. On the side facing the earth, these jets are 
terminated by narrow expansion fans of the slow MHD mode 
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Fig. 1 . The Levy-Petschek-Siscoe model of magnetopause recon- 
nection. Inbound crossing of Ogo 5 just south of the reconnection line 
is shown. 

and centered at the reconnection line. As the plasma in the jets 
flows across these fans, it expands isentropically, i.e,, its pres- 
sure, density, and temperature gradually decrease to zero, 
while the magnetic field and the flow speed tangential to the 
magnetopause increase. The model which in the simple form 
described above assumes a vacuum magnetosphere, has been 
analyzed quantitatively by Yang and Sonne nip [1977]. 

Figure 2 shows schematic polar plots of the predicted behav- 
ior of the tangential and normal magnetic fields during cross- 
ings of the magnetopause and slow expansion fan north and 
south of the reconnection region when the magnetosheath field 
is antiparallel to the magnutospheric field. The orientation of 
the mutually orthogonal a/^es, labeled 82, and ZI3, is as 
follows. The axis is approximately due north, the axis is 
due west, and the axis is normal to the magnetopause and 
pointing away from the earth. 

Observed Magnetopause Structure 

The magnetic field observations to be discussed here were 
made from the satellite Ogo 5 during the large magnetosphere 
erosion event on March 27, 1968 [Aubry et aL, 1970; 1971], 
during which the magnetopause traveled inward with the 
satellite for a period of approximately 2 hours. The total 
number of clearly identifiable magnetopause crossings during 
this event exceeded 20, occurring at low solar magnetospheric 
latitudes in the midmorning sector of the magnetosphere. The 
structure of these crossings varied greatly, usually bearing little 
resemblance to the theoretical structure described in the pre- 
vious section. The crossing to be discussed here is therefore 
exceptional. It was also studied by Aubry et al [1971; Figure 
1 5], who arrived at conclusions rather dififerent from those to 
be given here. 

A time plot of the magnetic field data from the GSFG 
magnetometer during the crossing is shown in Figure 3. The 
magnetic field is presented in terms of its magnitude |fi|, solar 
magnetospheric latitude, 0, and longitude, The crossing 
from the magnetosheath into the magnetosphere was ex- 
tremely rapid, the entire magnetopause region being traversed 
in about 4 s. Presumably this unusually short time is caused by 
rapid motion of the magnetopause past the satellite. It gives 


confidence that the satellite may have nearly captured a snap- 
shot of the local magnetic structure of the magnetopause. 
The magnetic data for this crossing have been subjected to 
minimum variance analysis \Sonnemp and Cahill, 1967] in 
order to determine the principal axes of the variance ellipsoid, 
defined by the matrix 

^a/3 “ ^0(8^ — 8d8^ 


In this formula an overbar denotes an average over the data 
set. Also, 8d and 8^ (a, = I, 2, 3) are the Cartesian com- 

ponents of an individual measured field vector. The principal 
axes and the associated unit vectors Ni, Nj, and N3 correspond 
to the directions of maximum, intermediate, and minimum 
variance (Ai, X2, and A3) in the corresponding field com- 
ponents, 8xy 82, and respectively. When the component 82 
of the field along N3 remains nearly constant during a crossing 
and when, additionally, A3 « A2, one may interpret N3 as the 
vector normal to the magnetopause and 8^ as the normal 
magnetic field component. Such an interpretation is not en- 
tirely unique but it is consistent with a locally one-dimensional 
model of the magnetopause structure. In such a model the 
normal magnetic field component must remain constant as a 
direct consequence of V*B = 0. Adopting this interpretation, 
the orientation of the right-handed orthogonal triad (Nj, N2, 
N3) relative to the magnetopause is then approximately as 
described in the caption of Figure 2, 

Results from the minimum variance analysis of five nested 
data segments within the magnetopause crossing in Figur^3 
areshown in Table 1, along with the estimated error in 
the average value of 82, and AN3 in the orientation of N3, as 
discussed by Sonnenip [1971], Sonnenip and Ledley [1974], and 
Sonnemp [1976], Note that the error estimates do not include 
systematic errors caused, for example, by spacecraft fields or 
by changes in the orientation of the boundary during the 
crossing. 

It is seen that the results from all five data segments are 
consistent in all essential respects. As expected, the estimated 
errors are large for the innermost data segment 1. However, 
for the outermost segment 5 they are very small. The high 
accuracy results principally from the large separation between 
the smallest and the intermediate eigenvalue, A3 and Ag, respec- 
tively. If the estimated vector error in the magnetopause nor- 
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Fig. 2. Predicted polar plots of the magnetic field during magneto- 
pause crossings north and south of the reconnection line. The field 
components Bu 82, and are directed along the orthogonal right- 
handed unit vector triad Ni, N2, N3. The vectors Ni and N2 are 
tangential to the magnetopause and due approximately north and 
west, respectively, while N3 is normal to the magnetopause and points 
away from the earth. 
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Fig. 3, Time record of the magnetic field during Ogo 5 inbound 
magnetopause crossing on March 27, 1968, at 1818:58 UT. The mag- 
nitude \B\ of the field vector and its solar magnetospheric latitude 
angle 6 and longitude angle (p are shown {6 > 0 for fields with a 
northward component, ^> > 0 for fields with an eastward component). 
The sampling rate is 56 complete field vector measurements per sec- 
ond. The interval A-E consisting of 240 points comprises data segment 
5 in Table I, Approximate satellite location (GSM): X ^ 55000 km, Y 
^ -*47000 km, Z ^ 4000 km. 

mal vector N3 is denoted by AN3 it is seen that j ANa’Ni | = 
0.01 rad an^ | AN3*N2 | = 0.03 rad for data segment 5. Fur- 
thermore, Bs = 8.1 ± 0.4 7 (1 y = 1 nT) for that segment. 

We conclude that the minimum variance method in this case 
yields a very accurate determination of the magn^opause 
normal vector N3 and of the average field component ^3 along 
it. We also conclude that B 3 is significantly different from zero, 
a result also obtained by Aubry et ai [1971], using data from 
the UCLA magnetometer on board Ogo 5. The value of ^3 
obtained by the latter authors was ±12 7. The reasons for the 
difference between this result and the value ±87 obtained in 
our study are not known. However, for nearly half of the 
magnetopause crossings during the March 27, 1968, event the 
minimum variance analysis, applied to the zero level corrected 
GSFC data, yielded normal magnetic field components < 
Afia, a result that we think is not random and would be 
unlikely to occur if large zero level errors or spacecraft fields 
contaminated the data. The total zero level correction, ob- 
tained by comparison with the GSFC rubidium vapor magne- 
tometer, was about 2.5 7; the spacecraft field is unknown but 
was about 1 7 during preffight testing. On this basis we are 
confident that measurement errors contributed less than 2 7 to 
the total uncertainty in the values of ^3 determined from the 
GSFC data. We also emphasize that the differences between 
the UCLA and GSFC measurements are too small to have a 
significant influence on the results to be discussed in this paper. 

In order to compare the observed magnetopause structure 
to the predictions of the Levy-Petschek-Siscoe model in Figure 
2 we present the data in the principal axis system Nj, Nz, N3, as 


polar plots of Bi versus B 2 and Bi versus The result is 
shown in Figure 4. It is seen that the observations bear a 
striking similarity to the theoretical prediction for a crossing 
south of the reconnection region. Points A and E correspond 
to the magnetosheath and the magnetosphere, respectively. 
The segment A-B in the left-hand diagram in Figure 4 may 
represent the rotational discontinuity in which the tangential 
field vector rotates by approximately 180°, Note that the field 
magnitude does not remain constant during the rotation. The 
interval B-C, where the tangential field undergoes little direc- 
tion and magnitude change, may correspond to the traversal of 
the narrow wedge of uniform flow and field immediately on 
the magnetospheric side of the rotational discontinuity. The 
interval C-D may correspond to the slow expansion fan in 
which the field magnitude increases with little change in field 
direction. The segment D-E, in which the field magnitude 
remains more or less constant but the field direction adjusts to 
the final magnetospheric orientation at E, has no counterpart 
in the MHD model. 

Note that the polar plot of the tangential field in Figure 15 
of Aubry et al. [1971] i.hows only a portion of the field rotation 
(A-B), Mainly on the basis of the changing field magnitude in 
the rotation, these authors concluded that (A-B) was not a 
rotational discontinuity. 

The time durations of the various portions of the crossing 
are listed in Table 2. It is seen that the segment A-B occupies 
the largest portion of the total crossing time. On the basis of 
the normal vector orientations, Aubry et a/. [1971, Figure 16] 
have interpreted this crossing and one occurring about 2 min 
earlier Aubry et al,, 1971, Figure 13] as being the result of 
a broad indentation in the magnetopause, traveling tailward 
along the surface with a speed comparable with the magneto- 
sheath flow speed, Fq. In such a model the angle between the 
flow vector, taken to be tangential to the average magneto- 
pause, and N3 is about 55°. Thus the thickness 5 of any portion 
of the magnetic field record of duration Ar is 

5 = K)(cos 55°)A/ 

The last column in Table 2 lists the thicknesses of the various 
magnetopause segments calculated from this formula and with 
Vq = 200 km/s. On this basis, it is seen that the total thickness 
of the entire niagnetopause region is about 470 km, an esti- 
mate that exceeds the one given by Aubry et al. by a factor of 
2. This discrepancy derives from the use of different total 
crossing time A/. 

With a proton gyroradius of about 70 km, say, the entire 
crossing is no more than 3^4 gyrodiameters thick. The individ- 
ual features B-G, G-D, and D-E have thicknesses less than this 
gyroradius but much greater than the electron gyroradius. In 
such circumstances one might expect considerable deviations 



TABLE 1. Results of Minimum Variance Analysis of the GSFC Magnetic Data From the Ogo 5 
Magnetopause Crossing on March 27, 1968^ 1818:58 UT 


Segment 


Eigen Values, 7^ 

Normal Field 

Error in Ng 

B Vectors 

A, 

A2 

A3 

± A^3, 7 

ANg-N, 

ANg-Nz 

1 

48 

76 

1.5 

0.15 

3.4 ±3.4 

0.02 

0.13 

2 

96 

267 

23 

0.58 

10,0 ± 1.2 

0.02 

0.06 

3 

144 

463 

83 

1.3 

8.1 ±0.7 

0.02 

0.04 

4 

192 

570 

124 

1.6 

8.0 ±0.5 

0.02 

0.03 

5 

240 

722 

130 

1.9 

8,1 ±0.4 

0.01 

0.03 


The five nested data segments are centered at 65938239 nis UT 
Fifty-six vector samples per second. 
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Fig. 4. Polar plots of the observed magnetic field for the Ogo 5 
magnetopause crossing on March 27, 1968, at 1818:58 UT. Format is 
the same as that in Figure 2, and the curves should be compared with 
the predicted curves for a crossing south of the reconnection line. The 
dashed circular arc, of radius Rc and centered at represents 

the best fit of (10) to the data segment A-B. Data segment 5 in Table 1 
is plotted. The components of the vectors Ni, Ng, and N, on the GSM 
Cartesian axes (A^, T, Z) are (-0.3794, 0.1470, 0.9135), (-0.8467, 
-0,4532, -0.2787), and (0.3730, -0.8792, 0.2964), respectively. 

from simple MHD theory. The principal observed dis- 
crepancies are (1) the field magnitude does not remain con- 
stant during the ‘rotational discontinuity’ (A-B) and (2) there 
is a field direction change (D-E) at the inner edge of the ‘slow 
expansion fan’ (C-D). These effects are discussed in the next 
two sections. 

Structure of the Rotational Discontinuity 

By considering the conservation laws across a rotational 
discontinuity, Hudson [1970, 1971, 1973] has shown that in a 
plasma with nonisotropic pressure the magnetic field magni- 
tude and the plasma state need not remain unchanged across 
the discontinuity. On the other hand, the quantity (px + 

2po) does remain constant, from which we conclude that in the 
present case the pressure px perpendicular to the magnetic field 
must increase by Apx = 3.6 X 10“*° N/m^ as the magneto- 
sheath plasma flows across the rotational discontinuity, At a 
nominal density of 30 particles per cubic centimeter such a 
pressure increase would correspond to an increment ATx - 0.9 
X 10° °IC in the perpendicular temperature Tx, which does not 
seem unreasonable. In reality the density as well as the temper- 
ature would change across such a discontinuity. 

Unless particle reflections from the downstream side are 
important, an increase in the perpendicular pressure across a 
rotational discontinuity requires nonadiabatic particle behav- 
ior and may occur as a result of the scattering of particles off 
fl uctuating electromagnetic fields or as a result of steep magne- 
tostatic field gradients, Hudson’s analysis does not yield any 
information about the internal structure of a rotational dis- 
continuity, but neither of the above effects is likely to be 
important unless the thickness of the discontinuity is small as 
it is indeed observed to be in the present case. Unfortunately, 
no adequate self-consistent theory seems to exist for rotational 
discontinuities with thickness of the order of the ion gyro- 
diameter. However, even if such a theory did exist, it would 
probably not be directly applicable to the present case. Indeed, 
it is unlikely that even Hudson’s jump conditions should be 
used across the segment A-B. The reason is that the extreme 
proximity of the slow expansion fan (segment C-D in Figure 4) 
to the rotational discontinuity (A-B) is bound to influence the 
structure of both. Some ions will spend part of a gyro orbit in 
the rotational discontinuity and another part in the expansion 
fan, so that conditions in the segment B-C do not approach the 


gyrotropic equilibrium assumed in deriving Hudson’s jump 
conditions. 

Clearly then, it is not justifiable to treat the various segments 
of the crossing as separate entities, Nevertheless, it is useful to 
consider the following ‘pseudotheory’ v/hich probably con- 
tains the essential ingredients which influence the magnetic 
siTucture of a thin rotational discontinuity. 

Let us adopt a one-dimensional model of the rotational 
discontinuity in which 

B = (Bi(xs), B2(.Y3), Ba) 

E ^0 

and in which the plasma flows across the current layer from 
the magnetosheath into the magnetosphere. The frame of ref- 
erence chosen for this analysis is such that the flow is exactly 
antiparallfel to the magnetic field on both sides of the disconti- 
nuity. In the satellite frame this is not the case, and in that 
frame E ^ 0. However, nonrelativistically, the magnetic field is 
the same in both frames. Note that the axis points from the 
magnetosphere outward into the magnetosheath, so that the 
plasma flow has a component in the negative Xa direction. We 
also assume the electrons to be cold so that they move through 
the magnetopause by sliding exactly along the magnetic field 
lines. In other words, we assume that the electron length scales 
and drifts are unimportant in establishing the structure of the 
field rotation (A-B). It then follows that the electron current 
may be written 

je = Kf:;.B/|B| (1) 

where w, and B are the local electron density, electron 
speed, and magnetic field, respectively, and where e is the 
magnitude of the electron charge. Conservation of electrons 
yields 

Je, ^ neVeBs/lBl = const = /loe^o^a/lBol (2) 

when «o, 1^0 and Bo are the electron density, electron velocity, 
and magnetic field, respectively, in the magnetosheath. Note 
that Bs is constant, that the number density of protons is 
everywhere the same as that of electrons, and that the velocity 
i?o in the magnetosheath is common to protons and electrons, 
so that the net current in the jCg direction vanishes. Equation 
(2) may be solved for wi;e/|B| , and when the result is sub- 
stituted into (1) there results 

j. = «(^t>oB/|Bo| (3) 

The proton current is more complicated to deal with. For 
layers that are much thicker than the proton gyroradius the 
proton motion consists of the guiding centers sliding along the 
field lines, thus effectively canceling the electron current given 
above, but with first-order drifts across the field lines which 


TABLE 2. Time Duration A/ and Thickness 6 of Various Segments 
of the Magnetopause Crossing 


Segment 

A/, s 

6, km 

A-B 

3.12 

358 

B-G 

0.24 

28 

C-D 

0.56 

64 

D-E 

0.20 

23 

A-E - : : 

4.12 

473 


See Figures 3 and 4. 

The thickness is calculated from the formula 5 - >o(cos 55°) A/ 
with Ko = 200 km/s. 
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then account for the net current. Structures of this type have 
been analyzed by Su and Sonnerup [\96^]. For layers of thick- 
ness comparable with the ion gyroradius, one expects only a 
partial cancellation of the electron current plus a drift current 
parallel to the layer which is much larger than that described 
by first-order orbit theory and which exhibits much less de- 
tailed spatial structure across the layer. Denoting the latter 
current by jp and accounting for the partial cancellation of the 
electron current by the factor A (0 < ^ ^ 1), the two nontrivial 
components of Ampere’s law become 

dBo , , 

Aio Ji = = Mo(7pi+ \ fio | ) 

(4) 

Mo 72 = = Mo(7p 2 4- I Bol^ 

where Mo is the permeability of free space. These two equations 
may be combined into the complex form 

dBt 

iapoBt ^ -ipoipt (5) 

where / = {— 1V^^ a = /cwo^Uo/ 1 Bo | , and 
Bt Bt -f iB2 

jpt =7pi + iJp2 ^ ■ 

If we consider jpt as a known nonhomogeneous term in the 
linear equation (5) for Bt, we find the following solution to 
that equation: 


For «o ^ 30 cm’® the thickness 5 ab = 358 km given in Table 2 
is obtained with k = 0.37. 

It may be added that the greatest difficulty in developing 
self-consistent Vlasov theory or two-fluid theory of the rota- 
tional discontinuity is associated with obtaining aperiodic so- 
lutions. Exact solutions for infinite wave trains have been 
known for a number of years [Bell, 1965; Lutomirski and 
Sudan, 1966; Sonnerup and Su, 1967; Abraham-Shrauner and 
Feldman, \977]/\ bQ only known aperiodic solutions are those 
obtained by Eastwood [1972, 1974] by numerical integration of 
the Vlasov-Maxwell equations. They are restricted to the fire- 
hose limit so that no net flow of the plasma across the layer 
takes place. Further, the magnetic field component is identi- 
cally zero. An aperiodic adiabatic solution again with zero net 
flow but with Bz 0 has been given by Francfort and Pellat 
[1976]. Oblique cold plasma whistler soliton solutions are also 
known Kellogg, 1964], but these only allow field rota- 
tions of exactly 27 t. 

3. One may attempt to make further use of (7) by assum- 
ing the complex proton current to be approximately constant 
and equal io jpt in the rotational discontinuity. One then finds 

Bt ^ ~7pt/oi -h (Bm +7pi/«) exp {-iapoXi) (10) 

In the B 1 B 2 plane th^s formula is represented by a circle of 
radius Rc = (Bm +7pt/«) centered at B^ = ^p./cc, Bc^ = 
~Jp 2 /^^ The best fitting circle, shown by a dashed curve in 
Figure 4, is centered at B^^ = -8 7, 5^, = +5 7, so that with 
a = kn(^Vo/\Bo\ and |Bo| = 39 7 


Bt = (exp -/aMoXa) - fioi ' jpt(x)(c\p iapoxjdxl (7) 

where Smi is the complex field at X3 = 0. 

The following observations should be made concerning the 
formula (7): 

1. The complex exponential multiplying the right-hand 
side describes rotational behavior of the magnetic field vector 
in the B 1 B 2 plane with a nonconstant field magnitude given by 
the parentheses. The sense of the rotation of the field vector is 
the same as the one observed in the magnetopause crossing 
examined in this paper. It also agrees with earlier observations 
[Sonnerup and Cahill, \96S; Sonnerup and Ledley, 1974] and 
with the theory of and Sonnerup [i96S], Clearly, this polar- 
ization is the direct result of the electron current flowing along 
the magnetic field. 

2. The rotation of the field vector is periodic rather than 
aperiodic. Thus the solution obtained describes an infinite 
stationary wave rather than an isolated rotational discontinu- 
ity. However, by assuming that approximately one half period 
of this wave can be used to describe the structure of a rota- 
tional discontinuity in which the field rotates by 180®, one 
obtains a thickness 5 given by 

X ^ JL. = 1^0 1 /Q\ 

Mo« pt/toeVok Vo k 

where is the Alfven speed in the magnetosheath and \i is 
the ion inertial length, i.e., 

X/ = (9) 

Since apart from the correction factor for nonisotropic pres- 
sure, Uo ^ Uaqi one concludes that the lower limit for the 
thickness of the rotational discontinuity is of the order of ttx,, 
i.e., about 132 km for n^ = 30 cm"®. This occurs when k - I. 
As /: decreases below unity, the thickness gradually increases. 


= 0.21 n^evtik . 

”0.13 V ^ 

In the field given by (10) the curvature drift current derived 
from first-order orbit theory may be shown to be 

V'p.r,,. = Jpr + y'p. = — exp (-/«MoXa) (12) 

where Un is the ion velocity along B and the angle 6 = is 
shown in Figure 4. In calculating the average proton currentjp 
fromjhis expression by integration over jCa, it is seen thatjp^ 
and y’p2 are proportional to the averages of (cos i)/B* and 
(- sin 5)/5^ respectively. Thus it is seen th^ y'pj is negative 
regardles£ of the variation of with §, while jp^ - 0 when fic, 
== 0, and y’pj > 0 when < 0. We conclude that the directions 
of the proton current components given in (11) are compatible 
with the directions expected from the first-order curvature 
current. The magnitude of the first-order drift current is far 
less than that given by (1 1 ), which is not surprising. Computer 
simulations of particle orbits in thin rotational discontinuities 
show that protons may become semitrapped and execute 
meanderlike orbits with very large displacements parallel to 
the layer. 

The present model in which jpt is constant within the layer 
and then changes abruptly to zero at the edges is by necessity 
crude. Thus we do not wish to overemphasize the importance 
of the detailed agreement between the circle and the measured 
segment (A-B) in Figure 4. 

4. It would appear that the present model does not pre- 
clude the occurrence of an electric field £*3 along the X3 direc- 
tion within the current layer. Such a field may be produced by 
miniscule deviations from charge neutrality in the layer. It 
leads to ExB/£® drifts of both electrons and ions but not to 
any additional currents. and Sonnerup [1968] have shown 
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that in the approximation of first-order orbit theory, ^ 0 
must occur in a rotational discontinuity in a hot plasma. 

Edge Current 

Given the extreme thinness of the segment C-D (see Table 
2), it would be unreasonable to expect the MHD theory of 
slow expansion fans to represent accurately the structure of 
this portion of the magnetopause. In the interpretation offered 
in this paper the reason for the small thickness is that the 
crossing occurred very close to the reconnection line. In this 
low-latitude region of space, one might also expect finite gyro- 
radius effects at the magnetospheric edge of the expansion fan: 
the protons will penetrate approximately a proton gyroradius 
further into the magnetospheric field than the electrons. Since 
in the reconnection model the plasma in the expansion fan is 
moving away from the reconnection line nearly parallel to the 
magnetopause, a proton current due south would occur at the 
inner edge of the expansion fan C-D in Figure 4. This current 
is expected to induce a magnetic field component along the 
negative B 2 direction on the magnetospheric side, as is indeed 
observed in segment D-E. The magnitude of this field com- 
ponent is seen to be ^82 ~ 21 y. In the edge-current model 
described above we have 

(13) 

where 6 de is the thickness of the segment D-E and it has been 
assumed that the particle density n of the streaming plasma 
rises linearly from zero at the magnetospheric edge (point E in 
Figure 4) to the value «D at the point D. 

The small difference in field magnitude between points A (39 
7) and E (44 7) indicates that a substantial amount of magne- 
tospheric plasma is also present. However, this plasma is 
nearly stagnant and does not contribute substantially to the 
edge' current. With 5 de = 23 km and v - 900 km/s, the latter 
number being taken from the compressible magnetopause re- 
connection model of Yang and Sonnerup [1977], (13) yields «d 
== 10 cm“^ For a very thick expansion fan (C-D), one might 
expect «D - 0. However, in view of the thinness of the entire 
fan and edge current layer, 5 cd + — 87 km, a density of 10 

cm^® at point D does not seem unreasonable. 

Parker {\961 a, b) has discussed the generation of a trans- 
verse magnetic field component in the magnetopause when the 
plasma is streaming along the field. This effect, which has been 
analyzed quantitatively by Su and Sonnerup [197 1, Figure 4], is 
in all essential respects equivalent to the edge effect discussed 
here. 

It should also be noted that the above discussion is relevant 
only very near the reconnection line. At higher latitudes the 
magnetic field configuration is such that the magnetosheath 
electrons may penetrate deeper into the magnetosphere than 
the protons, simply by flowing along the magnetic field lines 
(see Figure 1 ). Such a situation would lead to an induced field 
AB 2 of the opposite sign. 

Conclusion 

As is evident from the previous discussion, it is not entirely 
appropriate to analyze the present crossing in tenns of sepa- 
rate MHD structures such as a rotational discontinuity and a 
slow expansion fan separated by a wedge of uniform flow and 
field. The total thickness of the crossing is too small for such a 
procedure to be acceptable, and the crossing should be consid- 
ered as a single unit. For example, the ‘expansion fan’ segment 
(C- D) appears to be a recovery from the weak magnetic field at 
the inner edge (B) of the rotational segment (A-B) rather than 


an actual increase in field magnitude above the magnetosheath 
level as is assumed in the Petschek model. Nevertheless, if one 
asks what magnetopause structure should be expected at a 
location south of, but very near, the reconnection line of a 
standard MHD reconnection model, then one answer, a rea- 
sonable one, would be in terms of a structure of the type 
discussed in this paper. In other words, the observed crossing 
may be thought of as an embryonic Petschek-type structure 
which at larger distances from the reconnection line evolves 
into a full-fledged MHD magnetopause reconnection configu- 
ration. 

The preceding arguments do not imply that one should 
think of reconnection as occurring in a quasi-steady fashion 
over the entire frontside of the magnetosphere. Rather, it is 
tempting to suggest that reconnection was confined to the 
tailward-convecting magnetopause indentation associated 
with this crossing and the one immediately preceding it (which, 
however, had a much smaller normal magnetic field com- 
ponent, Bs = “2.7 7) and that the indentation itself was the 
direct result of local erosion of the magnetosphere caused by 
reconnection. 

Finally, it is noted that the magnetopause structure ob- 
served for the crossing discussed here appears sufficiently lami- 
nar so that it would be worthwhile to attempt a self-consistent 
particle field simulation of this structure on the computer. 
Such a simulation would give a valuable check of the magne- 
topause thickness and it would indicate whether a time- 
independent equilibrium is in fact possible. 
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I . Introduction 

Magnetic-field reconnection has been proposed as a basic energy-conversion process 
which may occur in many parts of the universe. Its primary function in the cosmic 
scheme is to prevent the build up of excessive amounts of magnetic energy in asso- 
ciation with intense electric current sheets formed in highly conducting plasmas. 
Tlie reconnection process is thought to cause a relaxation of such configurations, 
cither partially or completely, and cither continuously or sporadically, toward their 
lowest energy (current-free) state. The magnetic energy released during reconnec- 
tion is converted into kinetic and internal energy of the plasma. The process causes 
the transfer of magnetic flux and plasma from topological cells with excessive flux 
to cells deficient in flux. This fact provides the basis for a precise definition of 
reconnection to be given in subsect. 3.4. Reconnection is also often referred to as 
magnetic field merging or magnetic field annihilation but, as will be seen, the three 
terms should not be used synonymously. 

Figures 1-5 show examples of cosmic current sheets where reconnection may 
occur. Figure 1 represents the field produced by two pliotosphcric dipoles which 
gradually move tow'ard each other [94]. In the absence of reconnection, a current 
sheet of increasing length forms between the dipoles in the highly conducting solar 



Fig. 1. Qualitative time sequence for two dipoles moving toward each other on the solar sur- 
face. A cunent sheet A— B develops during timc 0 < r < fo- Rapid reconnection sets in at f « to 
and relaxes the configuration toward a potential field in the short time c/q. 


RfPlBii 


» 




Fig. 2. Cunent^heet formation caused by stretcliing of magnetic loop on the sun (after 

Carmichael (14J). 
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Fig. 3. Sector structure of the interplanetary magnetic field in the ecliptic plane as observed by 
IMP'l in 1963. Positive and negative signs indicate the direction of the measured interplanetary 
magnetic Held away and toward the sun, respectively (Wilcox and Ness (123 1). The regions of 
inward and outward Helds are separated by a current sheet which intersects the ecliptic along 
the four spiral lines. (Copyright by American Geophysical Union, 1965.) 
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atmosphere above them. If reconnection suddenly sets in, the magnetic field may 
relax toward a potential one, as Indicated in the last picture of the sequence. This 
represents a possible, perhaps even plausible, meclianism for a solar flare [57,921. 
Figure 2 illustrates current-sheet fonnation caused by the stretching of magnetic 
loops on the sun during rapid plasma ejection [l4]. Fig.3showsmtcrplanetary mag- 
netic sectors with different polarity [123] separated by a current sheet. Figure 4 
shows the magnetopause current layer, formed as the solar wind presses the 
intcrplanctaiy magnetic field against the terrestrial field, as welt as the tail current 
sheet, resulting directly from tangential stresses exerted by the solar wind on the 
magnetic field in the two tail lobes. The topology sliown in the figure was first 
proposed by Dungey [32]. Figures shows the magnetic field configuration 
expected for a rapidly spinning planetary magnetosplierc such as that of Jupiter 
[48,67], 

All of the above examples, and many possible other ones, such as supernova 
remnants [70], accretion disks [76], and galactic dynamos [87] , illustrate cosmic 
situations in which magnetic field reconnection may occur. However, we do not 
know with certainty that the process does in fact take place in any or all of these 
geometries. And if it does take place, we still do not know much in detail about its 
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Fig. 5. Magnetic field configurations for a rapidly spinning magnetosphere containing low* 
energy plasma (Gleeson and Ay^ord [48]). (Copyright by American Geophysical Union, 1976) 


dynamics. Are both continuous and sporadic reconnection possible, and if so, what 
are the plasma parameters and geometries in which these two modes arc to be 
expected? What are the conditions for onset of reconnection? What is the energy 
conversion rate? In spite oftwenty years oftheoretical effort, recently summarized 
in a brilliant manner by Vas 3 diimas [1 18] , as well as several laboratory experiments 
[10,31,84,85,114] and computer experiments [4,42,43], no universal agreement 
exists concerning the answers to most of these basic questions. Even in the most 
recent literature, opinions about the cosmic occurrence of the process range from 
full acceptance [118] to outright rejection [2,3]. On the other hand, there is con- 
clusive evidence that reconnection occurs in tokamaks and other fusion devices as 
an end product of the resistive tearing-mode instability [12C -122]. 

One of the difficulties with the cosmic reconnection research effort to date is 
that to a large extent it has lacked the detailed integration of theoretical and 
experimental work essential to the effective advancement of our knowledge 
concerning the process. On the one hand, an extensive but rather abstract body of 
theoretical work exists [118], concerned primarily with the steady-state process 
and utilizing the fluid description. The latter is likely to be inadequate for the 
analysis of certain critical aspects of the process. On the other hand, laboratory 
experiments [10,11] indicate the importance of sporadic reconnection. However, 


the plasma parameters jn these experiments arc sufficiently different from those 
prevailing in mo.st cosntic applications so as to pose serious difficulties in the appli- 
cation of the laboratory results in cosmos. A wealth of indirect observational 
evidence in the terrestrial magnetosphere, both at tlic magnetopause and in the tail, 
suggests that if the process occurs, it is likely to do so sporadically rather than con- 
tinuously. In current observational magnetospheric work, the reconnection process 
is often invoked to account for a great variety of observations but with little effort 
to check theoretical predictions in detail- or to consider alternate interpretations. 
The result is that the observational case for the occurrence of the process in the 
magnetosphere is not as solid as it might be. For other astrophysical applications, 
the situation is even worse. 

On balance, our best opportunity for learning about reconnection as a viable 
cosmic energy conversion process is likely to be in the earth’s magnetosphere. It is 
difficult to account for the overall dynamic behavior of the magnetosphere without 
tavoking limc-dcpeiidcnt transfer of magnetic ilux from closed to open field lines 
and vice versa. And such transfer is one of the principal features of the reconnec- 
tion process. The magnetosphere offers the unique advantage of permitting in situ 
plasma and field observations with probes that are much smaller than relevant 
plasma length scales. Thus an intense magnetospheric observational program with a 
focus on reconnection, coupled with a theoretical effort aimed at the geometries 
and plasma parameters prevailing at the magnetopause and in the magnctotail 
would seem to have high potential for success. What is learned about reconnection 
in the magnetosphere may then be applied to other cosmic systems which do not 
permit in situ observations. It is seen that a research effort focused on magneto- 
spheric reconnection is likely to lead to significant advances in our understanding of 
many other astrophysical and cosmic problems. 

It is the purpose of this chapter to provide a concise qualitative summary of the 
present state of reconnection theory and observations, with special reference to the 
earth’s magnetosphere, and to brittg into focus a number of specific problems and 
questions concerning the reconnection process in its magnetospheric application 
which should be studied both theoretically and obscrvationally. The organization of 
the paper is as follows. First, a number of basic concepts are introduced via a quali- 
tative discussion of steady two-dimensional reconnection in sect. 2, and of possible 
nonsteady and/or three-dimensional configurations in sect. 3. With this background, 
the more detailed technical discussion in subsequent sections can be presented in a 
compact fashion. Specifically, sect. 4 deals with the external flow region, which is 
usually described in terms of the fluid appio.ximalion. Section 5 discusses one-fluid 
and two-fluid approaches to tlie plasma dynamics in the diffusion region, which is 
the site of the field reconnection process itself, and in which plasma microinstabili- 
ties arc likely to be important. Section 6 discusses possible mechanisms for the 
generation of finite resistivity in the diffusion region and for the onset of reconnec- 
tion. Section 7 contains a brief summary of present observational evidence for or 
against magnetospheric reconnection. Finally, sect. 8 provides a summary of out- 


¥ 


■V-l, 


■ rV- 








ffi^«?fnffi-f -apti^^ffa^y tQiiiamni 

-' f - e- ' ,/fe ji- : f -■:! ^ ■: I S" I 1 "~"i U: ' * 


52 BM,O.Sonnerupf Magnetic field reconnection 

standing problems along with certain recommendations concerning the organization 
of future reconnection studies. 

Three comments should be made about the scope of the cliapter. First , it does not 
attempt to provide a historical perspective. Rather it is organized to elucidate basic 
physical principles and recent significant approaches to the development of ade- 
quate theories of cosmic reconnection. Second, the paper docs not attempt to cover 
all direct and indirect evidence for or against reconnection in the magnetosphere, 
on the sun, or elsewhere in cosmos. Third, the paper docs not deal with applications 
in tokamaks and other Jaboratory devices where the physical boundary conditions 
are such that spatially periodic behavior results. It should be stressed, however, that 
vigorous interaction between fusion plasma physicists and cosmic physicists on the 
problem of reconnection is likely to be of substantial benefit to both groups. 

. 2. Plane steady-state leconnection: a qualitative picture 

In order to develop an understanding of certain basic features of magnetic field 
reconnection, it is desirable first to e.xamine the simplest possible qualitative model 
of the process. To this end, consider the two-dimensional time-independent 



Fig. 6. Bade plane reconnection configuration. Solid tines are magnetic field lines; dashed lines 
are streamlines. The shaded region at the center is the diffusidn region. 


B.U,0, Sontiemp, Magnetic field reconnection 53 

electromagnetic field configuration shown in fig. 6. The magnetic field ^ Is con- 
fined to the jcy plane and has a hyperbolic (X-type) null point at (he origin. An 
electric field of the form £ = £*0X15 present along the direction perpendicular to the 
plane of the figure. Since V X £=^ 0 in a steady state, and since partial derivatives 
with respect to 2 are assumed to be zeroi;it follows that£o is independent of x and 
y, i.e., the electric field is uniform. This electromagnetic field is imagined to be 
imbedded in an electrically conducting fluid or plasma. In the following subsections 
wc examine scvcrtil aspects of this physical model: flux transport, external plasma 
dynamics, nature of the region around the magnetic null point, and electromagnetic 
energy conversion. The treatment is qualitative. More detailed discussion of existing 
analyses is presented in later sections of the chapter. 

2.1. Flux transport 

It is well known [ 8 , 83 | that £ -£ = 0 is a sufficient condition for the flux transport 
velocity \)£ = to move points which are on a given magnetic field line at 

one instant in such a way that they remain linked by a field line at all later times. 
For example, points which at a certain instant are located on field lines C|C2and 
D1D2 in fig. 6 will move in such a way that at a later time they arc located on field 
lines CiC2 and DiP^t respectively. Thus, a set of points, originally located on a 
field line and subsequently moving with t>£>, may be thought of as representing a 
“moving field line”. This fact explains the use of the term flux transport velocity 
for Wjp. Note K at the reconnection process may be discussed entirely without 
reference to moving field lines and that indeed the latter concept might become 
invalid if substantial electric fields parallel to the magnetic field should develop. 
However, in tlic present simple model no such parallel fields occur except in the 
region very near the magnetic null. The use of the concept of moving field lines is 
then just another way of referring to the electric field £q. In this model, the use of 
the term “reconnection” to describe the process is best understood in terms of 
moving field lines. As the lines C1C2 and DjD2 move withUf^ toward each other 
througl] positions CiC^ and D1D2 they ultimately reach location C'/C2 and D'|1>2 
where the lines meet at the origin. The surfaces through these lines and perpendi- 
cular to the plane of the figure arc called scparatriccs, because they separate 
families of field lines of different topological origin. When the lines have reached 
this critical position, they appear to be cut and reconnected so that at still later 
times they are connected as CV'D|" and as shown in the figure. It is evident 

that the reconnection may he thought of as leading to a transport of magnetic flux 
from flux cclIsQ) and (D across the separatrices into cells (|) and 0 . 

2.2. External plasma dynamics 

Up to this point the description of the reconnection process has contained lio 
reference to plasma dynamics. Indeed, the process may well have been imagined to 
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occur in a vacuum. In such rn instance, or if the ncld configuration is imbedded in 
a weakly conducting plasma, few restrictions exist on the magniludc of 
the magniludc of uh- And the magnetic field will be equal to, or nearly equal to, a 
vacuum configuration with an angle ft of nearly between the intersecting 
separatrices at the origin. The coupling between the electromagnetic field and the 
plasma is weak or absent. But in virtually all cosmic applications of interest, the 
field configuration in fig. 6 would be imbedded in a plasma of high electrical 
conductivity. Indeed, in many cases Coulomb collisions may be considered entirely 
absent and the condretivity , if such a term is to be used, is associated with plasma 
turbulence and/or inertial and gyro effects, occurring near the magnetic null. Away 
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from that point, the coupling between the B field and the plasma is strong and the 
piasma dynamics of the process will have dramatic effects in determining the 
detailed magnetic field configuration and perhaps in limiting the magnitude of the 
electric field /Tq. We now outline some basic features of the plasma dynamics of the 
reconnection process* 

First, it is observed that in a collision-free plasma the guiding centers ofcliargcd 
particles move with some velocity Up under the influence of the electromagnetic 
field in fig. 6. In the drift approximation, which is expected to be valid, except In 
the immediate vicinity of the origin, the component of Up parallel to the Jty plane 
and perpendicular to^ is identical with the flux transport velocity Thus, in that 
plane, and as long as £ • ^ = 0, the magnetic field lines may be thought of as moving 
with the plasma or vice versa. We note that the simplified magnetohydrodynamic 
description also yields this result in the limit of an mfinite electrical conductivity. 
The region away from the magnetic null in which plasma and fields move together 
is referred to as the convection region. Qualitatively the plasma motion is the one 
shown by the velocity ^ rows in fig. 6. Plasma approaches the origin along the 
positive and negative x axes and leaves along the positive and negative y axes. The 
motion may be the result of an external electric field Eq applied between capacitor 
plates at 2 = ±//. Alternatively. £*0 ^ polarization field created by an 

impressed plasma flow, specified in terms of a prescribed inflow rate at large ixl 
values, say. The details of the overall flow and field configuration will depend on 
these and otiier boundary conditions in a manner discussed in subsect. 4.1. However, 
all MHD models are expected to have in common the occurrence of large-amplitude 
standing waves in which the plasma is accelerated into the exit flow along the ±y 
direction, as shown in fig. 7. In incompressible analyses, these waves arc Alfve'n 
waves: in compressible flow they arc slow shocks approaching the switch-off Hmil. 
The occurrence of these standing wave patterns is related to the fad that the 
propagation speed of these modes is very small in directions nearly perpendicular to 
the magnetic field. Thus, by arranging the angle between the wave norma! and the 
upstream magnetic field to be sufficiently near 90°, the wave front can remain 
stalionar>v even for very small inflow speeds along the ±x direction. The set of 
waves divides the flow field into two inflow regions and two outflow regions. These 
regions do not coincide exactly with llie four flux cells in fig. 6. Because the separa- 
trices arc located upstream of the standing waves, parts of cells (|) and © overlap 
the inflow regions. 

Tlic standing waves contain concentrated electric currents, directed along the z 
axis, as shown in fig. 7. The JXB force associated with these currents serves two 
purposes: it balances the difference in perpendicular momentum and in pressure of 
the plasma across the shock, and it accelerates the plasma in a direction tangential 
to the shock. It should be emphasized that currents are by no means confined to 
flowing only in the wave fronts. Distributed currents /, may occur througliout the 
flow field. In particular, as will be shown in subsect. 4.1, the mrent distribution in 
the inflow region may influence the reconnection process in a crucial way. 
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An approximate balance of the magnetic shear stress at the shock and the exit 
momentum flow • yidds 

PiUjUj -BxBJpq (2.1) 

where pj is the plasma density in the inflow. Further, Ui;^i and U 2 >^ 2 » are inflow 
and outflow speeds and magnetic fields, respectively. They are related via 

EQ-V1BX-V2B2. ( 2 . 2 ) 

If Uj is eliminated between eqs, (1) and (2), we find 

«2-^i/(MoPi)'^ = Wai (2.3) 

and 

A/A,=y,/UA, -B 2 IB 1 , _ (2.4) 

Thus it appears that, regardless of the inflow speed, the exit speed U 2 is always of 
the order ofthe Alfven speed uA|> based on inflow conditions. Also, for fixed i?i, 
the magnitude of the magnetic field in the exit flow increases with increasing 
Alfven number A/a , in the inflow. When A/aj = 0, B 2 -O and the configuration 
reduces to a current sheet. When A/aj 1 the two fields are approximately equal. 

In steady-state reconnection models, the inflow Alfven number A/aj is com- 
monly used as a measure of the reconnection rate. 

For very small values of A/aj, and in a collisionless plasma, the plasma ejection 
along the ±y axis, postulated in the model in fig. 6, may become gradually replaced 
by an ejection at r ~ and z ~^h, respectively, of electrons and positive ions 
meandering in the current layer, as suggested by Alfven [1] and discussed further 
by Cowley [23]. The charge separation effects in that case lead to an electric field 
which is a function of the coordinate z. This limit will not be dealt with in the 
present paper. 

23, Region near the magnetic null 

The preceding discussion has dealt with plasma motion away from the magnetic 
neutral point at the origin In fig. 6. Let us now briefly consider the region imme- 
diately adjacent to that point. As the ongin is approached, the flux transport velo- 
city ve tends to infinity. Thus it is evident that the plasma can no longer move with 

* In this calcutation it is assumed that theplasma has a negligible velocity component along the 
y direction as it enters the shock. Tins assumption is not always valid. See subsects. 4.1 and 4.2. 
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U/. ill !hc .vr plane. In fact, as the plasma approaches the origin from both sides it 
imisl be brought to res! for syimnctry reasons. In hydrodynamic terms, the 
magnetic neutral point is also a double stagnation point. The region in which the 
plasma velocity deviates significantly from w/? is referred to as the diffusion region; 
its dimensions are denoted by 2x* and 2y* as indicated in fig. 6. in this region finite 
conductivity effects of some type must come into play, allowing the current 
density to remain finite at the null point forFo#0. Three main possibilities exist. 

(i) In a coilisional plasma with large but finite cicctrical conductivity o, the half 

width Jc* of the diffusion region is expected to adjust itself in such a way that a 
balance is established between the rate of magnetic fiux convected into the 
diffusion region and the rate of diffusion of the flux through the semistagnant 
plasma in the diffusion region. The ratio of these two transport rates is measured by 
the magnetic Reynolds number Thus we expect R^ cs. I, i.e., Jt* is 

of the order of the resistive length : 

X* ^{Pqov\)~^, (2.5) 

We note thatx* decreases with increasing conductivity and increasing Uj. Since Ui = 
ve-EqIBi, Bi being the magnetic field at (jc-±jc*,y -0), increasing Uf corre- 
sponds to increasing assuming 17 1 to remain fixed. 

(ii) In a collision -free plasma one might expect the value of x* lo be determined 

instead by the scale of the particle orbits near the null point. Four such scales may 
be of relevance: the electron and ion gyroradil and the electron and ion inertial 
lengtl^s. Further discussion of these scales is presented in sect. 5. An equivalent 
eicclrical conductivity may be imagined in this case, such that (he effective 
residence time of a particle (an electron or an ion) replaces the usual collision time 
r in the expression o- ne^rlm (w = particle mass). This residence time Is found to 
be Inversely proportional to U| so that JC* becomes independent of U| 

and lienee ofFo for fixed Bi , For further discussion, see subsect. 6.1 . 

(iii) In each of the above two cases, the current density or the gradients in the 
diffusion region may become sufficiently large to cause plasma microinstabilitics. 
The resulting plasma M'rbulcnce will lead to a reduction in the effective conduct- 
ivity, as discussed in subsect. 6.2. 

Whether the plasma dynamics in tlie diffusion region is described in a continuum 
fashion, i.e., by use of an effective conductivity, or in terms of individual particle 
orbits near the magnetic null point, il ls easy to see that the net current / in the 
d iffusion region will be along the positive z axis so that £ */> 0. Thus the diffusion 
region, along with the entij^ slipck system, acts as a dissipator of electromagnetic 
energy. 

*We note that the overall conservation of mass in the diffusion region yields 

Pi«i>’*=P20jx*, (2.6) 
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which may be combined with eqs. (2.2) and (2.4) to yield 

^A, =y|/yA, -(j>2lP\){x*ly*), (2.7) 

Assuming the density ratio P2/P1 to vary only moderately with A/aj, wc sec that 
the diffusion region is very elongated along the >^ axis for small A/aj values. Addi- 
tionally, in a collisional plasma the thickness X* increases with decreasing A/a|, as 
shown by eq. (2.5): 

( 2 . 8 ) 

Combining eqs. (2.7) and (2.8) it appears thatx* "^A/aIiJ'* a collisional 
plasma [case (i)] while const., y* ^ AfJ^\ in a collision-free case dominated by 
inertia! resistivity [case (ii)] . 

Finally, we estimate the separatrix angle a in the outflow (see fig. 6). Near the 
magnetic null point we may write: 

Bx^ay; By-bx; (2.9) 

where a and b are positive constants, and the angle a - 2 tan*"* (tf/^i)*/^. Estimating 
and we find by use of eq. (2.7) 

« = 2 K*7>’‘XP2/Pi)'^ 1 

^2lan-J[(p,/pj)'/*yirX.3, (2.10) 

indicating that the range of Alfven numbers A/a | front zero to (p 2 /p|)*^^ corre- 
sponds to an a range of zero to ^ir. The latter value corresponds to b-a, i t. to a 

current-free state, because /j — fl)/Po* 

2.4. Energy conversion 

The reconnection model described in this section serves as a steady-state converter 
of electromagnetic energy into plasma kinetic and internal energy. For example^ the 
rate of electromagnetic energy flow into and out of the diffusion region may be 
estimated as follows: 

inflow = 8y*^£oi9| /pq; outflow = 8 x*A£'o^ 2 //^V 

where the diffusion region has oecn taken to be a rectangular box with sides 2y*, 
2x* and 2/r. Thus the n^'t rate of inflow of electromagnetic energy is 


=(87*A£'o5l//ioXl - X*i?2/.y*^l) 
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which upon use of eqs. (2.2) and (2.7) may be written < ; 

H'™=-16vV/(/1?/2Mo)«'a,Ma,(1-/W^,P.^ (2.M) 

It is evident from this approximate expression that the energy conversion rate has a 
maximum at some value of the reconnection rate A/a | intermediate between 0 and 
a niaxnmim value, which in the present approximate set of relations appears to be 
A/ai = (Pz/Pi)*^^ ll^3t Il'em “ 0 l^oth for A/aj = 0 and for A/aj = 

For the latter value of A/a j, the configuration near the null is current-free and sym- 
metric a = ^7r). In such circumstances one may expect (P 2 /Pi)*^^ — I- Thus 

^ I appears as a theoretical upper Hinii for the reconnection rale (based on 
plasma conditions at x=x*,y-0). It is, however, by no means assured that 
boundary conditions at large distances or plasma processes in the diffusion region 
will always permit this upper Umil to be reached. 

The net rate of increase of kinetic energy of the plasma may be expressed as 
follows 

Wkt-PiViSy’h(iV2-^v])~&y’h(Bil2ixo)vAiM\,0-MA\) (2-12) 

and conservation of energy requires the difference fi'em - h'ke io be the rate of 
increase of the internal energy of the plasma, This latter rate may include 
thermal as well as nonlhcrmal parts, for example in the form of run-away electrons. 

The analysis given above applies to the diffusion region. But usually only a 
minute part of the total energy conversion occurs there, the main part taking place 
in the shocks. In approximate terms, the formulae (2. 1 1) and (2. 1 2) may be modi- 
fied to be valid for the entire reconnection geometry by replacing by where 
2L is the heiglit of the total configuration, as shown in fig. 6. Also, all quantities 
bearing the subscript I (which are evaluated at x = x*,V— 0) should be replaced by 
quantities bearing the subscript <», i.c. they should be evaluated at = 0. 
Depending on (he nature of the boundary conditions, the inflow may be such that 
A/a« differs significantly from A/a I (see subsects. 4.1 and 4.2). 

The phrase magnetic field annihilation has been used to describe the reconnec- 
tion process. In the light of the preceding discussion, this term appears appropriate 
only in the Hmit of small A/aj values where the magnetic field ^2 exit flow is 
small (or absent as in Alfv(:ii’s model, mentioned earlier [l,23]). Henceforth, anni- 
hilation will refer to situations where A/a j is sufficiently small so that the diffusion 
region occupies the entire length of the current sheet, i.e.,jv*>i^. By combination 
of eqs. (2.7) and (2.8) this is seen to occur for 0 <Mai ^ [(P 2 /Pi)/(/^®^Ai^)]*^^- 
In reconnection, energy conversion occurs on a time scale comparable to the 
Alfven wave time Ta^ L iu (assuming the inflow regions to extend to lx| — i), 
while in annihilation the scale is (taI'd)*^^* being the time for purely resistive 
decay of a current sheet i.e., to ^ is enormous in most cosmic applica- 

tions, so that reconnection rather than annihilation is required to account for the 
rapid energy release in solar flares, geomagnetic substorms, etc. 
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3 . Flux transfer in time*dependent and three'dimensifinal confi^:. rations 

The two-dimensfona) steady reconnection model outlined in sect, 2 is useful as a 
vehicle for introducing certain basic aspects of reconnection. But it appears likely 
that in any real cosmic applications of the process, three-dimensional and temporal 
effects are important, perhaps even dominant. For this reason it is useful to con- 
sider briefly a few reconnection configurations which incorporate these effects. To 
date, the plasma dynamics associated with such geometries has not been dealt with 
in a substantial way, so that the discussion is confined mainly to the electro- 
magnetic field topology and flux transfer aspects of the process. In the following 
subsections we describe the two-dimensional but time-dependent double inverse 
pinch configuration, a simplified steady-state three-dimensional magnetopause 
topology and a possible three-dimensional time-dependent magnetotail configura- 
tion. Finally, in subsecl. 3.4, a general definition of reconnection is given. 

3 J. Plamtime‘d€pendent geometry 

A plane vacuum magnetic field geometry associated with the double inverse pinch 
laboratory experiments {10| is shown in fig. 8. The X type magnetic null point is 
again located at the origin. The magnetic field is maintained by the currents / in the 
two metal rods at the center of flux cells 0 and @, and a return current 2/, 
flowing in the plasma along an outer envelope, which coincides with the outermost 
field lines in flux cell (5). In the experiments, the current / increases with time so 



Fig. 8. Field configuratioi; in double inverse pinch experiment (after Bratenahl and Baum 
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that magnetic flux is generated comimiajly at the two rods, i.c., in cells© and ©. 
If wc assume for a moment that no plasma is present, the flux in cell ©i> .'cases 
proportionately so that magnetic flux may be thought of as being transpo' ' from 
the rods into cells © and © and from there across flic separatrix info cell©. It is 
of interest to calculate the electric field responsibic for this flux transport The 
vector potential for the vacuum configuration is given by 

where the rod separation is 2c, the minor diameter of the return-current envelope is 
2dy and the radii rj and Tj are measured from the two rods as shown in fig. 8. Note 
that/4 = 0 at the envelope. In the experiment, the current/ and the envelope diam- 
eter both increase with time: !n a more general case, the rod separation might be 
imagined to depend upon time also. But for our purposes it suffices to consider the 
time variation of the current / and the diameter d. Then, the electric field is 

.rir, J 

Since at each instant remains constant on a magnetic field line, the instantane- 
ous electric field has (he same value on a given field line but its value changes from 
one line to another. In particular, on the separatrix it has the value 



Thus for increasing current / and diameter d, is positive as required for flux 
transport into cell ®. 

In the presence of a plasma, the field configuration is modified as folIows. The 
electric field now drives plasma currents in the vicinity of the magnetic null line, 
causing a .field deformation of the type shown fay the dashed lines in fig. 8. The 
separatrix intersection angle a falls below its vacuum value of These effects 
imply an excess of magnetic flux in cells© and©, a deficiency in cell @, com- 
pared to the vacuum configuration, which is the lowest energy state. Thus, a certain 
amount of free magnetic energy is stored in the system. However, at the same time 
a considerable amount of flux transport into cell© takes place. That is, reconnec- 
tion occurs continuously *. The principal difference between the present case and 
the steady-state model in sect. 2 is the spatial nonuniformity of the instantaneous 
electric field. This effect occurs because in the nonsteady case some of the fiux 


* By contrast, ref. [15] analyzes a hyperbolic-field'' collapse, where cr decreases from ^ to 0, 
without any reconnection. 
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transported in the xj' plane is being deposited locally, causing a field magnitude 
increase at each point. Associated with this flux accumulation, a plasma compreS' 
sion must also occur. But this would appear to be a relatively minor effect so that 
the steady model in sect. 2 may provide an adequate instantaneous description of 
the flow away from the rods and the return envelope. Thus the essential qualitative 
features of the reconnection flow may be obtained by examination of a sequence of 
steady-state configurations. 

Impulsive flux transfer events are observed in the double inverse pinch experi- 
ments. It appears that^ as the magnetic field and associated plasma currents near the 
null point grow, anomalous resistivity associated with ion sound turbulence sets in 
abruptly with an associated rapid increase of electric field and decrease of currents 
at the null point. The net result is a much more rapid flux transfer into cell @ and 
an associated relaxation of the entire magnetic field configuration toward its 
potential form with the separatrix intersection angle a increasing toward 
Evidently the stored free magnetic energy described in the previous paragraph is 
being rapidly converted into plasma energy. These events occur on a time scale 
much shorter than that associated with I, Hence it is unlikely that they may be 
described, even approximately, by a sequence of steady-state configurations. But 
the conditions for onset of such an event may perhaps be identified by examination 
of such a sequence. 

3.2. Steady three-dimensional geometry 

A thrce-dimensionai magnetic-field configuration of interest for steady-state 
magnetopause reconnection is obtained by the superposition of a dipole and a 
uniform field of arbitrary direction . This topology, shown in one cross section in 
fig. 9, has been discussed extensively in the literature [24,32,127]. Two hyperbolic 
magnetic null points X| and X 2 are formed in the plane containing the dipole 
moment vector and the uniform field vector. A basic topological property of such a 
null point is that many field lines enter it forming a separatrix surface and two 
single field lines leave it along directions out of that surface, or vice versa. The 
separatrix surfaces associated with Xi andX^ intersect along a circular ring located 
in a plane through the two points and perpendicular to the plane of fig. 9. This ring 
is referred to alternatively as a singular line, a reconnection or merging line, a criti- 
cal line, an X line, or a separator line. At a chosen point on the ring the magnetic 
field does not vanish in general, but it is directed along the ring. Only at X| and X 2 
is the field intensity zero. If the uniform field is exactly anti parallel to the dipole 
field a degenerate situation anses in which the magnetic field vanishes at each point 
on the ring. 

A schematic picture of the two separatrix surfaces is shown in fig. 10, in a con- 
figuration that may be rppropriate for magnetopause reconnection. The upper part 
of the figure shows a view in the antisolar direction of field lines on the separatrix 
surface associated with the null point X 2 : the lower part shows the same view of 
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Fig. 9. Field lines in the plane of the neutral points Xi and X 2 for a uniform magnetic field and 
a dipole field (Cowley 124]). Dipole moment vector at right angles to the uniform field. (Copy- 
right by American Geophysical Union, 1973) 

the X| separatrix. The total picture is an overlay of the two diagrams with the 
reconnection line connecting Xj aiid Xj. Bart of the solar-wind electric field Esw *s 
impressed across the configuration and must be sustained along the reconnection 
line. Thus, in the vicinity of that line a strong electric field component is present 
along the magnetic field. Unless special circumstances exist, such parallel electric 
fields do not arise in highly conducting plasmas. However, it is believed that the 
field lines on the separatrix and its immediate vicinity bend to become nearly 
parallel to the reconnection line extremely close to that line, as shown in fig. 10. 
Thus parallel electric fields occur only within the diffusion region which surrounds 
the reconnection line and in which finite resistivity effects permit their presence. 
Fig. 10 suggests that it may be possible to study reconnection in this geometry by 
use of a locally two-dimensional model which is then applied to each short segment 
of the reconnection line. Such a rriodel will be similar to that discussed in sect. 2, 
but with an added magnetic field component Thus the reconnection of 

fields that arc not antiparallel is obtained. Further discussion of such geometries is 
given in subsect. 4.4. The dynamics of the motion near the points X] and X 2 has 
not been studied to date. It may well be that these points mark the end points of t 
reconnection line segment on the front lobe of the magnetopause surface. 
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Fig. 10. Schematic of separatrix surfaces for magnetopause reconnection. Lower figure shows 
separatrix of the null point X|; upper Hgure that of X 2 . The two figures are to be super- 
imposed. 


Referring to fig. 4, which represents a cut through the earth’s magnetosphere in 
the noon— midnight meridional plane, it is seen that reconnection at the magneto- 
pause, as described above, serves to transport magnetic flux from the interplanetary 
cell(D and from the front-lobe magnetospheric cell (2) into the polar cap cells (|) 
and 0. 

3,3. Time-dependent three-dimensional geometry 

As a final example of reconnection geometries of cosmic interest, consider the 
magnetic-field topology associated with the formation of a reconnection bubble in 
the geomagnetic tail. The evolution of the field geometry in the noon— midnight 
meridional plane is shown in fig. 11. Note the formation of an X type and an O 



Fig. 11 . Fonnation of reconnection bubble in the geomagnetic tail. Schematic of field 
configuration in the noon-midnight meridional plane. 


type neutral point. The bubble originally has a very Small longitudinal dimension. 
As it grows in size in the noon-midnight plane, it also occupies an increasing 
longitude sector. The actual three-dimensional magnetic field topology of such a 
bubble is not known, but it may be represented sdiematically by an X type and an 
0-typc null line as in fig. 12. The points A, X, B and O in that figure all emerge at 
the same place at the time of onset of reconnection. Subsequently they move apart 
as the reconnection process continues and the bubble grows. An electric field exists 
along the reconnection line AXB but none, or almost none, along the O line AOB. 



Fig. 12. Three-dimensional sketch of reconnection bubble with the reconnection line along 
AXB and an O-type magnetic nuUline along AOB. 
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This field presumably has an inductive and an clcctrostaiic part which tend to 
cancel along AOB while adding along AXB. 

3.4, Definitions 

On the basis of the preceding discussion we now formalize the definition of several 
terms, used in the magnetic-field reconnection literature: 

(i) A separatrix is a surface in space which separates magnetic field lines belonging 
to different topological families. By necessity the separatrix is everywhere tangen- 
tial to the magnetic field. The field linc^constitu ting the surface originate at a 
hyperbolic neutral point in the field. 

(ii) A separator is the line of intersection between two separatrices or the line of 
intersection of one separatrix with itself. The separator is also called reconnection 
line, merging line, or X line| The terms neutral line, singular line, or critical line 
should be avoided, since they may refer to the 0-lypc topology as well. 

(iii) The diffusion region is a plasnia channel, surrounding; the separator, in which 
resistive diffusion, cdused by coHisional processes, turbulence, or intertial effects, is 
important. In : a highly conducting plasma, the diffusion region is irnbedded in a 
much larger convection region, in which magnetized plasma moves toward aitd 
away from the separator, in the inflow; and outflow regions, respectively, and In 
which dissipative effects arc confined to shocks. 

(iv) Magnetic-field reconnection is said to occur when an clcctric-ficId component 
£■() (induced or electrostatic) is present along a separator or a macroscopic portion 
thereof. It is proposed that the term magnetic-field annihilation be reserved for the 
case where the separator has degenerated (for dynamic puposes *) to a surface (c.g., 
the ^surface separating two half spaces containing antiparallcl uni-directional fields). 
The term magnetic-field merging may be taken to encompass both reconnection 
and annihilation. 

(v) The local instantaneous- reconnection rate at a chosen point on a separator is 
measured by the instantaneous magnitude of the electric-field component Eq along 
that line. It is desirable to express this rate in a nondimensional form by dividing 
the electric field by the product of a characteristic velocity and a characteristic 
magnetic field. The latter two quantities may be taken to be the Alfvdn speed 
and magnetic fielder at a chosen reference point, denoted by the subscript r, in the 
inflow, such as (Xr-x\ y^-O) or (Xr-Lr yr=0). Since f’o/^r represents a 
characteristic flow speed, the dimensionless reconnection rate takes the form of an 
Alfven number: 


In steady, two-dimensional (^* = 0, d/dz == 0), models the electric field Eq is 
* See comments in subsects. 2.4 and 6.1. 


constant tliroughout the jcy plane so that =^ £’o- With the reference point in the 
convection region, and on the x axis where = EqIBj is then thfc plasma flow 
speed toward the separator at the reference point and Mq is the local Alfvdn 
number, <^/o = il/Ar- In nonstcady flow, the electric field at the reference point, £*p, 

In general differs from £o» and A/o ^A/Ar- 

Vasyliunas [1 18j has defined magnetic merging as “the process whereby plasma 
flows across a surface that separates regions containing topologically different 
magnetic field lines”; he takes the magnitude of that flow as a measure of the 
merging rate. For reconnection in a highly conducting plasma, such that Rm = 
fiQOvL»l, the two definitions are essentially equivalent. However, the one 
adopted here, in terms of an electric field component along the separator works 
also for flows at arbitrar>^ ^m- If corresponds to the occurrence of flux rather than 
plasma transport across the separatrix, because flux transport is but an alternate 
way of referring to the electric field *. Note also that for the degenerate case of 
magnetic field annihilation there is no plasma flow across a separatrix. There is, 
however, an electric field -nd a corresponding magnetic flux transport. 

4. The convection region 

The plasma dynamics in the regions away from the immediate neighborhood of the 
reconnection line usually is described by use of continuum equations. Nonstcady 
solutions liavc not been found to date, which describe rapid configuration changes 
such as might be associated with impulsive flux transfer events in the double inverse 
pinch experiment (for a circuit modelj see Bratenahl and Baum, [1 1]). Three- 
dimensional solutions also have not been obtained, Hence the discussion in the_ 
present section is confined to steady-state plane reconnection. 

The incompressible assumption corresponds to the limit where is the 

ratio of plasma pressure to magnetic pressure. While this limit is invalid in most 
cosmic applications, it has the advantage of yielding simple analysis. Thus it 
provides an opportunity to study certain basic features of the reconnection flow 
without undue mathematical compl leal ions. We first describe two incompressible 
reconnection flows with fundamentally different behavior. Certain compressibility 
effects arc considered in the second subsection. The third subsection discusses 
asymmetric reconnection configurations, perhaps applicable to the magnetopause. 
The fourth subsection deals with the reconnection of magnetic fields that are not 
antiparallel, a common situation at the magnetopause. Finally, a partial sin^e- 
particle model is discussed briefly. 


* This equivalence is seen most clearly [119] by casting Faraday's law into the form of a 
conservation equation, viz., in subscript notation, aHj/ar 9/dx/ « 0, where is the 

antisymmetric (Levi-Civita) unit tensor. 


JWb = (£'o/fir)K. 
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4 A, Two incompressible symmetric flow models 

Figure 13, reproduced from Vasyliunas [118], shows a field and flow map for a 
reconnection model initially analyzed by Petsclick [90] and subsequently consider- 
ably refined and improved by Vasyliunas. The mode! contains a set of four Alfvcn 
discontinuities which compressible flow may be identified as slow-mode shocks 
and across which the piasma :s accelerated into the exit flow regions. Note that the 
plasma flow converges toward the x axis in the inflow and that the magnetic field 
intensity decreases on that axis for decreasing l.vl values. As pointed out by 
Vasyliunas, this behavior is characteristic of fast-mode expansion of the plasma as it 



Fig. 13. Pctschck’s reconnection model with A/a* ® O-l- Magnetic field lines (solid lines) and 
streamlines (broken lines) are shown. (Vasyliunas [118].) Fast-mode expansion in the entry 
flow (Copyright by American Geophysical Union, 1975.) 
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approaches the reconnection line. Because the fast-mode propagation speed is 
infinite in the incompressible limit, such expansion is by necessity an elliptic effect, 
that is, no standing expansion wavelets are possible. The maximum reconnection 
rate in this model corresponds to an Alfven number A/a| of about one in the inflow 
just adjacent to the diffusion region. But because of the increase in flow speed and 
decrease in magnetic field associated with the fast-mode expansion, the Alfvcn 
number, A/aoo? ^*1 large distances upstream is considerably less than unity. Values in 
the range 0.05 < Af a«. maximum rate are obtained (see ref. [118], 

fig. 12). Recently, Soward and Priest [MOJ have reexamined Pelschek’s reconnec- 
tion geometry by use of an asymptotic approach, valid away from the reconnection 
line. Their analysis in all essential respects supports the conclusions summarized 
above. 

Fig. 14, also taken from ref. [1 18], shows a flow and field map for a different 



Fig. 14 . Slow-mode reconnection model with A/a| “ (J-5. (Vasyliunas 1118].) Slow-mode 
expansion in the inflow is concentrated to waves from four external comers. (Copyright by 
American Geophysical Union, 1975) 
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model which is the sole noiisingular member of the similarity solutions 

derived by Yeh and Axford [131]. This model contains a second set of Alfvcn 
discontinuities located upstream of the slow shocks and originating at external 
comers in the flow, as shown in the figure. These discontinuities represent the 
incompressible limit of slow-mode expansion fans centered at the external corners. 
They cause a large iefleclion of the plasma flow away from the jc axis and a 
substantial increase i*T iflc|c! magnitude. It is now generally agreed that these discon- 
tinuities will not occur in any real situation. Rather tliey represent a suitable 
mathematical lumping of slow-mode expansion effects in the inflow. Tlic maximum 
reconnection rate in this model * is Mai=0 On the x axis this value 

remains constant, independent of jx|. However, this is a result of the lumping of 
the slow-modet effects. In a model where these effects arc spread over the inflow 
region the value ofA/^ ori the jf axis will decrease with decreasing |jc| in association 
with a decrease in plasma velocity and an increase in magnetic field. Thus, in reality 
it is unlikely that the inflow into the diffusion region can occur at as high as 
(1 +A^); more likely that value corresponds to the maximum A/a» at large |x| 
values. Further discussion of this point Is given in subsects. 4.2 and 5.1 . 

The two models discussed above represent two extreme sets of conditions in (he 
inflow: pure fast-mode and pure slow-mode expansion. In any real application both 
effects may be present. Vasyliunas [1 18J has pointed out tliat from a mathematical 
viewpoint the difference between the two models is related to the boundary condi- 
tions at large distances from the reconnection line. Ear Upstream, the fast-mode 
mode! is essentially current free and has a nearly uniform flow and magnetic field, 
while the slow-mode model contains substantial currents which bend the magnetic 
field lines and cause a deflection of the flow away from thc\x axis. Vasyliunas has 
further suggested that the former state of affairs may obtain when a demand for 
magnetic fliix originates at the current sheet itself (the plane) or in the exit flow, 
as may be the case in the geomagnetic tail, while the latter set of conditions may 
correspond to externally forced inflow such as at the magnetopause. In this 
context, it is worth noting that slow-mode expansion effects have been argued 
[132] to be present outside the subsolar magnetopause regardless of whether or not 
reconnection occurs there. 


4.2. Compressible symmetric models 


A detailed compressible analysis of the external region of Petschek's reconnection 
geometry is not available at present. On the other hand, the slow-modc expansion 
model has been iextended to include compressibility effects. An Imthermal analysis 
was given by Yeh and Dryer [130]. But the isothermal assumption leads to 

• The estimates given in sect. 2, viz., U 2 ~ waj (wi)max ~vaj» assumed a negligible flow 
component along the y axis as the plasma enters the shock. Such a component is present in this 
model, the result being that the exit flow speed »2 and the maximum inflow speed (uDmax 
both exceed uA| by a factor (1 ^ y/2). 
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Fig. IS. One quarter of a symmetric compressible slow-mode magnetic-field reirannection model 
for Af^i =0.7 and 0i^2uoPi/^l =5 ^Yang and Sonnerup 1125)). (Reprinted by permission 
from C.K. Yang and H.U.O. Sonnerup arid the University of Chicago Press. Copyright by the 
American Astonomical Society, 1976.) 


unacceptable entropy variations with decreasing entropy acrosT the shocks and 
increasing entropy across tlic expansion waves. More recently, an analysis has been 
performed by Yang and Sonnerup (125), which assumes isentropic flow in the 
inflow and uses tlie ordinary jump relations for slow shocks. It is found that the 
expansion-wave discontinuities in the incompressible solution do indeed dissolve 
into slow expansion far- centered at the external corners in the flow (see fig. 15). 
It might be tliought tliat the reflection of these fans in the jc axis, and the subse- 
quent interaction of the reflected waves with the shocks, shown schematically in 
fig. 16, may be treated exactly by the method of characteristics. However, it is 
found that the flow from region ® in the figure, across the last expansion wavelet 
and the innermost portion of the shock, cannot be deajt with without the inclusion 
of elliptic (fast-mode) effects. This is cxtrememly difficult to do. Thus, in the main 
part of their workp Yang and Sonnerup, after calculating the isentropic plasma and 
field changes across the fans, considered them to be lumped into a single discon- 
tinuity, i.e., they ignored (he reflection altogether. While such a {procedure is 
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Fig. 16. Schematic showing “reflection” of slow-modc expansion fan in the j: axis (Yang and 
Sonnerup (125j). [Copyright by American Astronomical Society (liinjvcrsity of Chicago Press, 

1970) 

perhaps justified in a first attempt to study compressibility effects in the external 
flow, it nevertheless seriously limits the usefulness of the resulting solutions. The 
width of the slow expansion fans in the inflow increases dramatically with increas- 
ing compressibility, he., with decreasing values of pj ^ 2poPi/B], so that for /3| ^ 1 
the lumping of the fan Into a single discontinuity is difficult to defend. Further- 
more, except perhaps for very large values, conditions immediately outside the 
diffusion region arc not adequately represented so that the solution may not be 
used to provide external boundary conditions for compressible matched diffusion- 
region analyses. However, the analysis is valid at large distances from the origin, and 
it is of interest to examine Its predictions concerning flow and plasma conditions in 
the exit regions, A^hen conditions typical of geomagnetic tail reconnection are 



dashed curves to the Soward -Priest ( 1 10{ analysis of the fast-mode expansion model. 


subsfiluicd, flow speeds in the range of 1000 kin sec“* arc calculated, in rough 
agreement * with observed proton speeds in the tail during energy-release events 
|39,63J. The analysis also predicts exit flow speeds considerably greater than the 
fast-mode propagation speed so that standing transverse fast shocks may be present 
in the two exit flow regions, causing a decrease in flow speed and an associated 
increase in plasma density, temperature, and in the exit magnetic field. 

Yang and Sonnerup [1 25] also calculated the change in plasma and flow proper- 
ties along the X axis in fig. 16, caused by the reflection of the slow expansion fan, 
but ignoring the elliptic effects mentioned earlier. The solid curve in fig. 17 shows 
the resulting relationship between the Alfven numbers and Af^^, in regions® 
and © of fig. 16, respectively. For comparison, the corresponding relationship for 
the fast-mode model, developed by Soward and Priest [MOJ, is shown by the 
dashed curves. It is evident that the different distant boundary conditions for (he 
fast-mode and the slow-mode models may lead to profoundly different inflow 
conditions into the diffusion region for the two models. 

* The agreement is however not lurncicntly detailed to support this particular reconnection 
configuration over others. 
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4.3 i A symmetric models 


A qualitative reconnection model for the asymmetric How and field conditions 
existing at the magnetopause was first dcscrihed by Levy ct al. [72,91]. In tliis 
model, shown in fig. 18, the magnet oshcath plasma is assumed to carry with it a 
compressed interplanetary magnetic field which is due soutli so that a field reversal 
results at the magnetopause (see fig. 4), In impinging on the earth’s field, the 
plasma encounters a wave system consisting of an intermediate wave (rotational 
discontinuity; large amplitude Alfve'n wave) followed by a narrow slow expansion 



pz 


Fig. 18. Upper hair of compressible stow-niodc model of magnetopause reconnection for ^ 
0.2 and p| “ 2 moPi/Bj = 2. The intermediate wave (IW) marking the magnetopause is shown as 
a dashed line. The slow mode expansion fan (SEF) is shaded. (Yang and Sonnerup [126].) 
(Copyright by American Geophysical Union, 1977.) 
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fan. The intermediate wave, which marks the magnetopause, accomplishes the field 1 

direction reversal and an associated plasma acceleration parallel to the magneto- 
pause and away from the reconnection line. The niagnclic-field magnitude remains 

constant across this wave but it then increases to its higher magnetospheric value in J 

the slow expansion fan across which the plasma pressure also is reduced to match g 

the pressure in the magnetosphere, taken to be equal to zero in the model. The ^ 

quanlitalive details of this model have not been worked out for fast-mode expan- ^ 

sion in the in How. But for the incomprc.ssibIc slow-modc expansion model, various 

types of asymmetries in the flow and fieldhavc been analyzed [26,82,106]. Certain 

compressible couterparts of these geometries have been studied by Yang [124] by 

use of the procedure of lumping slow-modc effects in the inflow, discussed in 

subsect. 4.2. In particular, the case of vacuum conditions in the magnetosphere has ! 

been reported in detail [126]. A typical resulting field and flow configuration is as 

shown in fig. 18. The model predicts the following features of the magnetic field: ■ 

(i) a small magnetic field component normal to the magnetopause; (ii) rotational . 

behavior of the magnetic-field component tangential to the magnetopause; (tii) a n 

gradual increase in magnetic field intensity inside the magnetopause. In terms of 

plasma flow, the model predicts: (i) How of magnctoshcath plasma normal to and 

across the magnetopause with speed equal to the Alfven speed based on the normal 

magnetic field component; (ii) tangential acceleration of the plasma to speeds of 





G9 03 27 18:16:56:34 UT 

Fig. 19. Polar plots of magnetic field at the magnetopause. Left hand figure shows the field 
components and Bi tangential to the magnetopause during an OGO-S crossing; right hand 
figure shows the nearly constant magnetic-field component 83 normal to the magnetopause. 
The field is given in units of y (ly = Int). Intermediate wave or rotational discontinuity is 
segment A1-A2 of the left-hand trace: the slow expansion fan is segment A2-A3. The segment 
A3--A4 may be caused by a finite gyroradius effect not contained in the MHD model. 
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the order of 500—750 km sec^* •• as it crosses the magnetopause; (iii) no change in 
plasma density or temperature as it crosses the magnetopause; (iv) an isentropic 
decrease in density and temperature and an associated velocity increase as the 
plasma expands in the slow expansion fan inside the magnetopause. To date, the 
predicted plasma behavior has not been observed. At various times some but usually 
not all of the predicted magnetic signatures have been seen [I08J. An example is 
shown in fig. 19. 

4A, Reconnecting fieldt forming an arbitrary an^le 


In his original paper on reconnection, Pc tsclick [90] observed that in incompres> 
sible twO'dimcnsional reconnection flow, a constant magnetic ficid component 
could be added to any solution without changing the flow or magnetic ficid config- 
uration in the xy plane (refer to fig. 6). Thus, it is possible to generate solutions 
that describe the reconnection of fields that form an arbitrary angle. This provides 
a link between the traditional cosmic reconnection models and the problems of 
reconnection in nearly force^free field configurations, such as the tokamak [120- 
122]. This procedure has provided the basis of a number of attempts to describe 
the dependence of the cross-magnelospheric electric potential difference on 
magnetosheath field direction, assuming the Tonner to be caused by magnetopause 
reconnection [49,58,59,107]. The result of these geometrical analyses is that the 
potential should have a principal angle dependence given by the function * 



J (gp/gj - cos 
I+(ao^i)^-2(^o/^t)cos^ 



(4.1) 


where 0 denotes the angle bctwceii the tangential field B^ in the magnetosphere and 
Bo in the magnetosheath. For cos 0 > no reconnection occurs and^^(O) sQ. 
Observations [13] indicate a low energy injection rate from the solar wind into the 
magnctospheric ring current system when 0 <\n, a result which appears com* 
patible with eq. (4.1). Observations [12] of the magnetosperic cusp location as a 
function of 9 also arc in qualitative agreement with this equation. 

Recently, Cowley [25,28] has pointed out that incompressible solutions exist in 
which ^2 is a function of jc and j?. Thus the assumption underlying eq; (4.1), of one 
and the same value of on the magnetospheric and the magnetosheath side of a 
typical magnetopause reconnection model, may be invalid, ^t is noted that this 
assumption corresponds to a situation where thei net current in the magnetopause 
(and in the diffusion region) is parallel to the separator. 

In the incompressible MHD approximation the equations describing the flow and 


field in the jo’ plane arc completcy uncoupled from the differential equations for 
the velocity component and for B^. However, as pointed out by Cowley [25], aii 
indirect coupling exists via the requirement that the line integral * dl = 0 for any 
closed loop which encircles the diffusion region. This requirement poses an addi- 
tional constraint on the shape of the diffusion region, a constraint that does not 
arise when B^ = 0, or B^ = const., and that does not appear to be automatically 
satisfied by the diffusion region of plane reconnection geometries. Thus it is not 
clear at the present time whether Cowley's criticism of eq. (4.1) has a firm founda- 
tion in incompressible MHD theory. But even if it doesn't, the use of a constant Bg 
in the real compressible magnetopause flow situation to construct reconnection 
geometries for arbitrary $ values remains hypothetical. Further theoretical study of 
this problem requires detailed analysis of compressible external and diffusion-region 
flows and appropriate matching of these flows at the edge of the diffusion region. 



• In ref. (49 J the functional dependence is [/(d)l*/^. *^ifr 20. Hill'* (591 coUisionlcss reconnection model. The magnetic-field change across the slow 

•• If the magnetosbeath neld magnitude exceeds the magnetosphere field, the subscripts o ind I >*><«* becomes weak fot fi a 2«ioP/S* -► 0 with the principal field reversal occurring in a current 

ate to be interchanged. layer at 
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4.5. Ofllisiofiless model 

No complete or nearly complete collisionless models for the external reconnection 
region have been developed to date. But certain partial results have been obtained 
by Hill [59]. He suggests that for small j3 values the principal field reversal is 
accomplished by a current sheet located on they axis, as shown in fig. 20, with only 
a small amount of field-line bending at the slow shocks. Hill does not treat the flow 
and field configurations in the inflow or in these weak shocks. Rather lie assumes 
that, away from the magnetic null point at the origin, the field lines fortii an angle x 
with the current sheet. He then proceeds to discuss the properties of the sheet. One- 
dimensional self-consistent Vlasov equilibria of such sheets have been obtained 
numerically by Eastwood [36,37]; an approximate analytic theory using the 
adiabatic invariant [34,65,104] associated with the meandering particle orbits in 
the sheet has also been given [38j. However, the result primarily used by Hill is 
obtained directly from the overall stress balance at the sheet, without reference to 
the sheet structure: in a frame of reference sliding along they axis (see fig. 20) with 
a speed such tliat the reconnection electric field Eq Vanishes, the usual firehose limit 
must apply. Assuming the rnagnetic moment of individual particles to be preserved, 
Hill shows that for inflow along the jc axis * this condition may.be expressed as a 
local reconnection rate 

£*o/ua^~( 1 -tt)*^sinxcosx, (4.2) 

where Ua — and.ff are the Alfven speed and magnetic field, respectively, 
in the region adjacent to the current sheet. Further, the anisotropy factor a of the 
incident particles is defined by 

(4.3) 

Note that or = 0 corresponds to isotropic pressure, I to firehose conditions 
(taking account of the fact that the total plasma density p is twice the density, pi„, 
of the incident particles). 

The preceding result is incomplete in that the angle x must be a function of the 
reconneetjoh rate. Also, the rate should refer to conditions on the x axis in the 
inflow region. Equation (4.2) is nevertheless interesting because it suggests that 
pressure anisotropy in the incoming plasma may be an important factor. In particu- 
lar, it appears that reconnection may cease altogetiier for a = I . 

Using the same approach, Hill has also considered the case of reconnection of 
fields of unequal magnitude and with a constant ^^2 component present. He obtains 
the formula (4.1) for the angular dependence of the reconnection potential. 

* Hill’s analysis also includes an unspecified velocity component uy of the incoming particles, 
which we have set equal to zero. 


The particle energization in a model of HilEs type is seen to be the direct result 
of inertia and gradient drifts in the current sheet, moving positive ions in the 
direction of the reconnection electric field, electrons in the opposite direction. It is 
also important to note that the energized plasma will be streaming out nearly 
parallel to the y axis, i.c., for small angles x, nearly parallel to the magnetic field on 
both sides of the current sheet. By contrast, the symmetric fluid dynamical models 
yield an exit plasma flow that is perpendicular to tlie weak magnetic field in the 
two exit flow regions (sec figs. 13 and J4 below). 

S. Fluid description of the diffusion region 

A complete theoretical treatment of the reconnection problem requires not only a 
self-consistent solution for the external flow, but also an internal, or diffusion- 
region solution which describes the essential dissipative pliysical processes operat- 
ing in that region, and which joins smoothly to the external solution. In the present 
section we review attempts to describe the diffusion region in terms of continuum 
equations which incorporate effects of plasma microinstabilities, if any, by means 
of an effective conductivity a. 

A brief discussion of onc-fluid theories is given in subsect. 5.1. In magneto- 
spheric applications of reconnection, the collisional resistive length (/ioWj)“* is 
much smaller than relevant inner plasma scales, such as the electron inertiaUength. 
Thus onc-nuid theory is directly applicable only if turbulent processes generate an 
effective resistive length which exceeds these inner scales. Bui even if that condition 
is not met, one-fluid theory serves the important purpose of providing an oppor- 
tunity for a careful mathematical treatment in one region of plasma-parameter 
space. 

Tlie two-fluid description of the diffusion region is dealt with in subsects. 5.2 and 
5.3 . The former discusses the importance of the electron inertial length in deter- 
mining the width 2x* of the diffusion region when no collisional or turbulent 
resistivity is present. In the latter section, the importance of Hall currents and of 
the ion-inertial length and gyroradius are discussed. 

5.L Ouc'fluid models 

Detailed studies of the flow and field configuration in the diffusion region, using 
onc-fluid magnetohydrodynamics, have utilized two approaches: scries expansion 
around the magnetic null point, and development of integral or “lumped” equa- 
tions for the entire diffusion region. Additionally, certain exact solutions exist. 

Series expansions have been given by Priest and Cowley [93] and by Cowley 
[27]. They found that in incompressible flow, and assuming analytic behavior at 
the null point, a plasma velocity of the form - —kiX, Vy = kiy to lowest order 

yields a magnetic field behavior of the form = ^2 By-k^ Jf, i-c., the field 

lines touch at the neutral point, rather than cross at an angle d#0. The latter 
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type of configuration may however be obtained by assuming a third-order, rather 
than a first-order zero in Ujc(x) at x = 0. Furthermore, Yeh [128] has shown that 
the flow and field behavior near the null point need not be analytic. Logarithmic 
terms are possible In the expansion. Whether or not such terms arc in fact present 
can be determined only by matching of the series expansion solution loan appro- 
priate external solution, which has not yet been done. It also appears that the 
inclusion of compresf’bility effects will invalidate the above-mentioned result of 
Priest and Cowley. Finally, Cowley [27] has shown that series expansions yielding a 
non-constant field Bjfxyy) are possible. But the question of whether such solutions 
may be matched to corresponding external solutions with non-constant (see ref. 
[25]) has not been resolved. 

The first lumped analysis of the diffusion region was performed by Parker [86] 
with application to Sweet’s resistive current-sheet model [113] of a solar flare. The 
analysis yielded the following expression for the reconnection rate in incompres- 
sible flow: 

-Va, (5.1) 

This result is obtained directly from eqs. (2,5) and (2.7) wim pj = P i - Originally, 
the formula was used with y* replaced byX, the half-length of the current sheet. It 
then describes field annihilation (see subsect. 2.4) and yields a value of that is 
far too small to account for an energy-release time of the order of 10^ sec in a solar 
flare, or for the observed potential difference of 20-100 kV across the terrestrial 
magnetosphere. Later Petschek [90] used the formula (5.1), now with 2y* 
represenfing the height of the diffusion region and with y*«X, to describe the 
diffusion-region flow in his model. The basic point made by Petschek is that cq. 
(5.1) determines, not the reconnection rate, but the height y* of the diffusion 
re^on. In agreement with the discussion in subsect. 2.3, y* is then seen to be propor- 
tional toA/xf- 

A more detailed lumped analysis was performed by Sonnerup [103] in an 
attempt to develop a diffusion re^'on solution for the slov/-modc reconnection 
geometry in fig. 14. The treatment is incomplete because it docs not take account 
of the momentum balance. Additional criticism has been offered by Vasytiunas 
[118] on the basis that the model implicitly assumes an abrupt switch from finiie 
to infinite electrical conductivity at the ou‘er edge of the diffusion region. Con- 
sidering the extreme simplification of the external flow in this model, with slow 
expansion effects lumped into a single discontinuity (see subsect. 4.1 and fig, 14), 
this latter criticism is probably not of serious consequence. The following qualita- 
tive results of the analysis are likely to remain valid for the slow-mode reconnection 
model in fig. 14: 

(i) A relationship exists between 3^d the magnetic Reynolds number 
which is of the form given by Parker, eq. (5.1), for small values of 
and which yields ^ Q for AfA, " (I +\/2).Thus, the dimensions of the diffusion 
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region shrink toward zero as A/a j approaches its maximum value. 

(ii) When the slow-mode expansion in the inflow is concentrated into discontin- 
uities as in fig, 14, the increase in magnetic field and decrease in flow speed, 
described in subsects. 4.1 and 4.2 must occur in the outer portions of the diffusion 
region. Thus, as the origin is approached along the ±x axis, the magnetic field first 
increases tlicn decreases to zero while |Ujc( decreases montonically to zero. Thus, in 
the outer portion of the diffusion region the field behaves almost in a frozen-in 
manner. This may account for the large diffusion region width x* found in the 
analysis (103,1 18J. 

Two exact solutions exist which describe flow near the magnetic null point. Yeh 
[129], obtained shock-free similarity solutions by assuming resistivity and viscosity 
to increase linearly with distance from the origin. It is not clear how such assump- 
tions can be reconciled with an exterior solution in which dissipative effects are 
confined to shocks. A different type of exact solution has been discussed by Priest 
and Sonnerup [95,109]. It describes an incompressible two or three-dimensional 
M HD resistive stagnation-point flow at a current sheet. The field lines are straight 
and parallel to the current sheet. Thus, purely resistive magnetic field annihilation 
without reconnection occurs, as illustrated in fig. 21. These solutions represent a 


y 



Ffgl 21. Magneitc field lines and streamlines for stagnation point flow, v = (-ktX,k2y,k2T), at 
a cunent sheet. The diffusion-dominated region is shaded (Sonnerup and Ftiest [ 1 C:/ ]). 
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Fig. 22. Nondimcnsional magnetic-field profiles /J (uQo/k | )“* * for the configuration 

in fig. 21 in the resistive limit (ot - ^ 2/^1 v ^3 - ^^^nc stagnation point flow for a ~ 1 , 
axisymmetric flow for a - The frozen field profile fora - 1 is shown by dashed line (Sonne- 
rup and Priest [109J). 

generalization of a case studied by Parker [88]. The resulting magnetic field profiles 
are shown in fig. 22. Three features are of interest. First, tlie cliaractcristic width of 
the current layer is of the order of the resistive length as expected. Second, the 
frozen field condition applies at large |.x| values and leads to a gradual increase in 
the field magnitude \By\ as [x| and jUj^j decrease. As resistivity effects become 
increasingly important reaches a maximum and then decreases to zero atx - 0. 
This is precisely the behavior described in (ii) above. Third, a nenconstant value of 
82 is possible. 

5.2, TwO'fluid effects: electron scale lengths 

In the two-fluid description of an electron— proton plasma, the ordinary Ohm’s law 
is replaced by a generalized form (see, e.g,, Rossi and Olbert [96]). 


£-+«X/f=-^+gf|^+V(/«+i,/)l-— V-Pe+— (5.2) 
a ne^Ldt J ne ne ' ' 



As pointed out by Vasyliunas, the terms on the right-hand side may introduce a 
variety of plasma scale lengths into the problem. The first term yields the resistive 
length \ = (pqOVi)^'; the second set of terms yields the electron inertial length 
the off-diagonal parts of the electron stress tensor term Pe 
yield the electron gyroradius; the last term, describing the Hall effect, yields (v/J 
Ui)X/ where = is the ion inertial length. The importance of the 

diagonal terms in has not been studied: with isotropic pressure and isentropic 
flow, these terms are cancelled identically by an electrostatic field. 

To illustrate the effects of the electron inertia terms we now generalize the 
stagnation-point flows discussed by Sonnerup and Priest [109] to include two-fluid 
effects. Assuming incompressible flow and diagonal stress tensors for ions and 
electrons, the flow and fields arc of the form: 

o = -xk |X + + zkjz; B = + zB^ix ) ; (5.3) 

where the quantities A:j, X '2 and are constants such that lfc| =1^2 +A: 3 . These 
assumptions lead to a Bernoulli-type pressure integral! 

P=Po +* 32 ^) -(2Mor'(5>+5?) 

of the momentum equation, and to the following components of the Induction 
equation [the curl of cq. (5.2)] : 

X’A-, .V (A-jX? - I/a/oo) b; - - My = 

X|*, a B;” + (A:3Xj - J/moo)B;'-A:,.yB:-*3Bj = 0. (5.4). 

It is seen that only the elect ron-iinertial length and the resistive length appear in eqs. 
(5.4). The Hall current term in cq. (5.2) is curl-free and is cancelled exactly by a 
Hall electric field ^^^(x). Thus the ion inertial length docs not appear. The solutions 
of eqs. (5.4) in the resistive limit are illustrated in fig. 22. In the inerlial limit, the 
odd solutions are shown in fig. 23 for various values of a = A' 2 //:i- As expected, the 
width of the magnetic-field reversal region is now of the order of the electron 
inertial length regardless of the flow rate. Further, it is observed that foro = I,i.e., 
for plane flow, the current density is logarithmically infinite at x = 0, a conclusion 
also drawn by Coroniti and Eviatar [22j. Thus some form of plasma microinstabil- 
ity or other effect by necessity must be present to reduce the current density to a 
finite value. 

The off-diagonal terms in the electron pressure tensor will provide a finite 
electron gyroradius correction to the preceding results. When the electron gyro- 
radius greatly exceeds Ac, it will replace the electron inertial length as the minimum 
width of the layer [118]. However, a detailed calculation of these effects is. difficult 
because the appropriate form of the off-diagonal pressure tensor terms is not known 
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Fig. 23. Nondimensional magnetic-field profiles for the configuration in fig. 21 in the clcctron- 
inenial limit [a = ^2/^1 » ^3 ~^1 ^^2^ 7 Plane siagnation-point flow fora = 

1 ; axisym metric flow for a = The frozen field profile for or “ 1 is shown by dashed line. 


in a field reversal region of width comparable to the orbit scale. 

Vasuliunas [J 18] pioneered the study of electron inertial effects in the diffusion 
region. In an approximate lumped analysis which neglected compressibility (later 
included by Coroniti and Eviatar [22]), off-diagonal stress tensor elements, and 
Hall-current effects, he showed that the diffusion region in the inertial limit (i.e., 
neglecting resistivity) is hyperbolic in shape with width 

x’(y)-iM\^y^ (5.5) 

The exact analysis given here in all essential respects confirms Vasyliunas* results 
for small values of It also provides magnetic field profiles By(x) whereas 

Vasyliunas neglected By within the diffusion region. Note that the formula (5,5) 
yields xXO) ^Xc;y* =^ Ae/AfAj, in agreement with the behavior quoted in subsect. 
2.3. [between eqs. (2.8) and (2.9j]. 

53. Two fluM effects: ion scale lengths 
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the ion inertial length and the ion gyroradius did not appear in the previous 
analysis. The former js introduced via the Hall current term J XBjne in cq.(5,2). lt 
may be cancelled by a Hall electric field only when it is curl free, which was the 
case for the stagnation-point in subsect. 5.2. The ion gyroradius as introduced via 
the off-diagonal terms in the ion pressure tensor. We now demonstrate that in plane 
reconnection flow these effects imply the presence of Hall-current components/^ 
and jy as well as a macroscopic flow v^x^y) and a field component y). 
Assuming d/3z = 0 and omitting electron inertia and pressure terms for brevity, the 
z components of the momentum and Induction equations arc, respectively: 
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where the particle density n and the conductivity o have been assumed constant. In 
the second equation, the Hall term (the last term on the right) is seen to be of the 
form (B ' i.e., it vanishes only wlicn the current remains constant along a 

field line. Otherwise it becomes a driving term in the second equation forcing values 
of Uj and B^ different from zero. Similarly, in the first equation it appears that the 
stress terms will usually force values of and B^ different from zero. These terms 
are expected to introduce the ion gyroradius as a characteristic length into the 
problem. However, it is difficult to discuss these effects in detail, because the form 
of the stress terms in a thin field-reversal region is not known. Thus, we confine 
attention to the Hall current term, (^ * V/^)/«e. While this term vanishes in the 
stagnation*point flow discussed in the previous section (and indeed along the x axis 
of any configuration), it cannot vanish throughout the diffusion region. We may 
estimate the magnitude of by approximately equating the third term on the left 
and the second term on the right in eq. (5.7): 

Vjc^BJ^x'^iMne^Bxhjzl^x. 

With Ux a^d 3/dx ^ 1/x* wc then find 

BjBi -hlx\ 


It is not clear how a diffusion region of the small physical dimensions implied by 
eq. (5.5) can be joined to an external solution with slow shocks of thickness [20] 
comparable to, or greater than, the ion inertial length. Thus we are led to ask why 


where Xj = (Wi/wie^)*/^ is the ion inertial length. Thus it appears that values of x* 
much less than Xi would give rise to unacceptably high values of Bg. A similar 
companson between the terms ^x^U 2 /dx and (nc)"* Bj^JJdx yields '^ X|/x*. 
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Again, x* « Xj leads to an unacceptshle result. Vasyliunas |1 19) has pointed out 
that off-diagonal electron pressure tensor terms, not shown in eq. (5.7), may 
possibly cancel the Hall term. However, ilicrc seems to be no obvious physical 
reason to expect such a cancellation. And problems with the off-diagonal stress ten- 
sor terms and P,,, i/ eq. (5.6) wouW 

Detailed analysis ot the effects described above is not available at present. 
However, a nonvanishmg field component would imply the presence of 

Hall currents /x ^Pq * dBjdy and iy^ dB^/dx in the diffusion region. The 
expected current flow and field pattern is shown schematically in fig. 24. The 
behavior of the By and B^ components indicated in the figure should be easy to 
identify in magnetjc-field vector measurements from a satellite which crosses the 
diffusion region at the magnetopause or in the geomagnetic tail. 

The reason for the appearance of the Hall current component/^ with the direc- 
tion shown in fig. 24 may be understood by noting that for o the generalized 
Ohm’s law [eq., (5.2)] implies that apart ifrom electron-inertia and gyroradius 
effects the magnetic field is frozen into the electron component of the plasma. 
Thus the electrons fiowing toward x == 0 are brought to rest over a distance of the 
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Fig. 24, Schematic piciure of Hall cunent loops in the diffusion region. Also shown is the trans- 
verse magnetic field induced by these currents. 
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order of the electron inertial length or the electron gyroradius. If the ions are 
similarly brought to rest over a distance coinparable to the ion inertia) length or the 
ion gyroradius, a relative motion of electrons and ions results, leading to currents 
/x in the direction sliown in the figure. Charge conservation then implies the 
presence of /^, as shown. 

The 2 component of the force / XB associated with the Hall currents also leads 
to an acceleration of the plasma in the ±z direction. This effect is caused by drift 
and meandering motion of the ions in the current sheet. It may be interpreted as an 
ion current in the layer. Indeed, ifx* ^Ai, principal current component /j is 
earned by the ions; ifx* it is carried by the electrons. 

It is evident that in the absence of plasma resistivity (classical or turbulent) the 
electron Icngtli scales must play a significant role in the diffusion region structure. 
But from the preceding discussion it appears possible that tfic ion length scales 
may dcterniinc ilic overall width 2v* of the diffusion region while the electron 
length scales give the size of the detailed structures of/x,/y, and /g ncarx * 0. From 
the preceding discussion, it is concluded that, even without plasma turbulence, the 
electromagnetic structure of the diffusion region may be far more complicated than 
previously assumed. 

6. Non-fluid effects in the diffusion region 

In niagnctosphcric and interplanetary applications of the reconnection process, 
collisional resistive effects in the diffusion region are negligible. Thus an effective 
resistivity In that region must dcriyc citlicr from inertia effects or from plasma 
turbulence. The former effects were dealt with in siibsccts. 5.2 and 5.3 from a fluid 
point of view. However, to develop a physical understanding of inertial phenomena 
in the diffusion region, it is useful to obtain an appro.ximate expression for the 
effective inertial conductivity. This is done in subsect. 6,1. Subsection 6.2 examines 
plasma instabilities which may generate steady-state turbulence in the diffusion 
region, but with details provided only for the ion-acoustic instability. Subsection 
6.3 discusses several threshold effects that might be of importance for the onset of 
reconnection and for the identification of situations in which reconnection may not 
occur. Particle acceleration in nonstcady reconnection is discussed in subsect. 6.4. 

It will become quickly apparent that most of the material in this section is 
speculative in nature. Different processes may occur in different applications. It 
appears that no systematic effort has been made to sort out which mechanisms 
dominate in different parts of plasma parameter space. 

In reading the present section, it will be useful to refer to the typical values of 
several physical parameters giyen in table 1. 

6J. Inertial resistivity 


The concept of inertial resistivity was first discussed, in the context of magnetic 
field reconnection, by Speiser [lOS]. The basic idea is that a particle spends only a 



Table 1 

Physical paramatcrs for various reconnection sites 
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finite amount of time in the diffusion region and thus can pick up only a finite 
amount of energy from the electric field £joi. The inertial conductivity is written as 

OmcTt = fie^rlrn, ( 6 . 1 ) 

where n is the average particle density in the diffusion region andT is the effective 
time available for acceleration in the electric field. The formula (6.1) may be used 
with cither the elcetron mass m = m^ or the ion mass /« = Wj depending on whether 
the diffusion region current is principally an electron current or an ion current. 
Both cases may be treated the same way so that the particle mass /ii will be left 
without a subscript. 

The average displacement, Az, along the electric field £*o, of a particle in the 
diffusion region may be obtained from a simple mass balance over a box of 
dimensions iif* X X Ar: ~ 

2niVi2y*^-nv22x*2y* = 2n2V22x*^ ( 6 . 2 ) 

where is the average particle speed along the electric field -Mlt • 

As before, the subscripts ! and 2 refer to conditiosn at the points (jc*, 0) and (0, 
y% respectively. From the first equality in cq. (6.2) we thus obtain 

f—nx*Jri2Vi. - (6J) 

When this expression for f is substituted into cq. (6.1) there results 

o\nt,t=n^e^x'lntVim . (6.4) 

Thus the inertial conductivity is very high for low reconnection speeds, Ui, i.e., 
when the configuration approaches a current sheet. This is the behavior referred to 
in subsect. 2,3. 

Expressing the basic balance of field convection and diffusion as ^ajnertt^i^* = 
I , which is valid for small reconnection rates *, we find 

x\ =0/1^) X, (6.5) 

where X| is the inertial length X| =(m//io«ie^)*^^. For m —iWe this result is in 
agreement with Vasyliunas’ formula (eq. (5.5)]. See also Coroniti and Eviatar (22]. 
But the calculation gives no clue as to whether the electron or the ion inertial 
length is to be used. 

It should be realized that the value of x* given by eq. (6.5) represents a lower 

* With inertial rcristivity, pure field annihilation i$ found to occur for OO/ai 
( compare subsect. 2.4). ”S mall reconnection rates” implies \\jL< Mji ^ « 1 . 
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Jimit; The calculation assumed the time T available for free acceleration in the 
electric field to be equal to the residence time of a particle within the control box. 
In reality T must always be less than the residence tinie because the niagnciie field 
docs not vanish within the entire box. Thus the inertial condiictivily is less than 
that given by cq. (6.4) and x* is correspondingly larger than in eq. (6.5). When 
diamagnetic currents become important, i e,, for 0^ -IpoPilB] > I , it inay be 
shown that x:* gradually approaches a magnitude of the order of the gyroradius 
instead. 

The previous estimate of -X* applies only for small values of the reconnection 
rate. To understand this fact, we note that the expression PqOViX* ~ 1, which was 
used in obtaining cq. (6.5), derives directly from Ohm’s law in the approximate 
form f^-oBg-oviBi with j^-PQ^dBylbxf^ large reconnection 

rates it becomes important to incorporate tlie term oX^ in the electric field, as 
well as the curvature term dBJby in the expression for/-. The latter effect leads to 
a muitiplicative factor (I —A/a Jpj/P 2 ) the right-hand side of eq.(6.5) (see cq. 

(6.7) below] so that the size of the diffusion region decreases tow'ard zero as the 
reconnection rate approaches its maximum value. Thus, for large reconnection 
rates, the required width of the diffusion region may be substantially less than the 
relevant plasma scale (the inertial length or gyroradius, depending on Pi). It is diffi- 
cult to reconcile such a situation with the nature of the particle orbits in that 
region. Therefore, it is conceivable that steady-state reconnection with inertial 
resistivity as the dominant effect in the diffusion region is not possible for large 
reconnection rates. 

The average electrostatic particle energization in the diffusion region may be 


obtained directly from the first equality in eq.(6.2); 

f = eAzEo =^nev2X* EolniVi . (6.6) 

' . -1 

But we also have ncU^=/ {bBylbx — bBJby). Approximating by 

Bilx\bBJby by B2ly*\n^ = W|/u 2 we find 

u ^nevz “ ( 6 . 7 ) 

and from eq. (6.6) 

] , ( 6 . 8 ) 


where the relations Eq = Ui^i and vx jvj ~ ^i/wai have been used. This 

formula agrees with eq. (2. J I). Maximum acceleration occurs for small reconnec- 
tion rates and densities: for inagnctospheric conditions 6 l~10kcV. It is 
emphasized that eq. (6 8) represents the average energy gain. A small number of 
particles moving nearly along the reconnection line may gain larger amounts of 
energy in the electric field £q. 
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6,2. Plasma turbulence 

A variety of plasma instabilities may .serve to generate plasma turbulence in the 
diffusion region and an associated turbulent conductivity Ofurb- now discuss 
such effects in an assumed quasi-sicady stale of reconnection. Onset effects are 
dealt with in subsect. 6.3 

The tearing instability, cither in its collisional resistive version [45,97,1 17] or in 
the collision -free cicciron-incrtial version [17,30,44,61,69] has been studied 
intensely in the context of reconnection. It generates a pattern of alternating X 
type and O type magnetic neutral lines in a current sheet. But most analyses of 
this instability pertain to current sheets with a vanishing magnetic-field component 
Bx perpendicular to the sheet. In other words, in the present application either 
tlic reconnection rate 3/^^ is very small or the magnetic field lines at tlie null 
point arc “touching” rather than “crossing”, as discussed in subsect. 5.1. Schindler 
[99] has pointed out that for Bx^Q the collisionlcss tearing instability may still 
proceed as long as the gyro-period Tg InmleBx of a particle in the field Bx exceeds 
tlic instability growili time Tq — CvVuch) where Ujh is the thermal speed 

and R/^ the gyroradius. This condition may be applied cither to electrons (electron 
tearing) or ions (ion tearing). In rough terms,, non-gyrotropic behavior of the 
particles is required for these instabilities to be possible. While the nature or 
existence of steady-state tearing turbulence docs not appear to have been esta- 
blished, one cannot exclude the possibility that such turbulence could be of 
importance in the diffusion region [19,21,47]. 

Parker [89] has suggested that interchange instability may serve to enhance the 
How rates in Sweet’s [113] current sheet model. In the geomagnetic tail, the 
instability would be driven, not in the diffusion region itself, but rather by the pres- 
sure gradient and field curvature in the near-earth section of the tail plasma sheet 
(sec fig. 4). A detailed analysis, including the impeding effects of the ionosphere, 
has been given recently by Kan and Chao [66]. It indicates growth times of the 
order of a few hours with ionospheric coupling, a few minutes without such 
coupling. The situation relative to the level of steady-state turbulence is not clear. 

Huba ct al. [64] have proposed tliat the lower-hybrid-drift instability may 
provide anomalous resistivity in the diffusion region. It appears that the threshold 
for this instability is sufficiently low to permit the diffusion region width to be of 
the order of the ion inertial length. 

Hacrcndel [50] has discussed the possibility that the electron-cyclotron drift 
instability, which has a current threslioid somewhat less than that of the ion- 
acoustic instability, may generate turbulence in a diffusion region of width 2x* 
equal to a few electron-inertial lengths. However, its importance has been ques- 
tioned by Coroniti and Eviatar [22] on the basis that the gyrocoherence required 
by the instability may not be available in the diffusion region. They also note that 
when the electron drift speed exceeds the threshold for the ion— acoustic instability 
the electron— cyclotron drift mode goes over nonlinearly to the ion— acoustic one. 
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The ion-acoustic instability has been proposed [9,22,41,101] as a likely agent 
for the generation of turbulence in the diffusion region. It will be dealt with in 
some detail in the remainder of this subsection. This is done for illustrative pur- 
poses and not as an indicator of a universal preference for this particular 
mechanism. On the other hand, the ion-acoustic instability appears in fact to occur 
in laboratory reconnection experiments [9,85] . But it is probably not relevant to 
magnetospheric reconnection. 

For a current-driven instability such as the ion-acoustic one to occur, the 
current density in the diffusion region must exceed a certain minimum value, 
corresponding to a critical current velocity i.e., / > neoj.. If Hall /cur rents, ate 
present, as discussed in subsect. 5.3, the total current must be considered. Here we 
shall confine attention to the component/,. According to cq. (6.7) we then find 

=-^fj 1 >«euc (6-9) 

A'o’f L P2 ‘J 

where, for the ion-acoustic inrtability 

Vc=(kTilm,y/^IATjT{). (6.10) 

Combination of eqs. (6.9) and (6.10) yields 

x’Ae<(2//3i)'/" n -(Pi/P2)Afi,I/(7V7y, (6.11) 

where ft = The; function f{TjT\) is shown in fig. 25. It is seen that 

is of the order unity for Te ^ Ti so that for small A/aj, and for ft of order 
unity or less, the critical diffusion region width is of the order of the electron iner- 
tial length. For large values of ft,x* must be considerably less than suggesting 
that only a sub-portion of the diffusion region may contain ion— acoustic turbulence 
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Fig. 25. Function /(7«/7i) * (*Ti/mc)*/^/uc where is the critical current velocity for onset 
of ion-acoustic instability (after Fredricks [40]). 
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(compare the logarithmic singularity in /, at x = 0, discussed in subsect, 5.2). 
Coroniti and Eviatar [22] indicate ^< 5 as a condition for their analysis to remam 
valid. For greater /3 values, the size of the turbulent region approaches the wave- 
length of the ion-acoustic turbulence. For A/^j of order unity the diffusion region 
must also be very small for ion-acoustic turbulence to occur. 

For high temperature ratios TJT\,{\\q instability may occur for x* considerably 
larger than Ag but probably not as large as \\ (see fig. 25). A large temperature ratio 
may pcriiaps be generated temporarily by electron run-away in a current sheet at 
the onset of reconnection (sec next subsection). For example, in the double inverse 
pinch experiment the collisional resistive length considerably exceeds (sec 
tabic I), so that run-away must occur to initiate the ion-acoustic instability. But it 
appears unlikely lliat a large temperature ratio TjTi could be sustained on a steady 
basis in a diffusion region of width much greater than Ag since in most of such a 
region the run-away would have to occur transverse to a strong magnetic field. We 
conclude that steady-state ion— acoustic turbulence, driven by the current com- 
ponent /„ is unlikely to be important unless the diffusion region width, 2x*, is of 
the order of the electron inertial length. At the same time it is observed that the 
Hall current component jy discussed in subscct. 5.3 (fig. zA) may be sufficiently 
intense to drive the instability in parts of a diffusion region of total width compar- 
able to the ion inertial length. 

Coroniti and Eyiatar [22] have examined the question of tlic turbulent satura- 
tion of the ion-acoustic instability in detail. They conclude that the current veloc- 
ity will remain close lo the threshold value given by cq. (6,10). The resulting weak 
steady-state turbulence is adjusted to give the value of turbulent conductivity 
required to satisfy — ^ with jc* given by the equality in cq. (6.1 1). On 
the other hand, common estimates of the effective electrical c^onductivity 
associated with the ion— acoustic instability, in a state of turbulent saturation, such 
as(sce,e.g., refs. [40,115]) 


’ We Lwpe ihte \Ttl J * 


where iOpe ~ {it e^/coWe)*^^, give a much too low value of the conductivity, even at 
the critical current velocity = other words, with reasonable reconnection 
speeds and with x* satisfying cq. (6.11), one finds pQOf^it^UiX* « \, vfhlch is 
impossible in a steady slate. In table 1, this fact is manifested by the inequality 
^turb ^e» where A|ujb the turbulent resistive length. 

6.3, Onset of rapid reconnection 


There is ample observational evidence relating to solar flares, to the earth’s magne- 
totail, and to the double inverse pinch experiment, to indicate that occasionally 
rapid reconnection is switched on in an abrupt, almost explosive manner. At the 
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earth’s magnetopause, if reconnection actually occurs thcrcViiie switch-on appears 
more gentle and may be a direct consequence of the interplanetary field turning 
southward so that the angle between the reconnecting fields exceeds some criticid 
vaJue (compare siibscct. 4.4). ft is natural to assume that the niore explosive events 
might be associated with a plasma instability and/or an abrupt decrease in the 
effective conductivity in a current sheet or in the diffusion region of a slowly 
reconnecting configuration. Five such possibilities, all speculative at present, arc 
mentioned below: 

Thermal instability, h has been proposed [16,19,56] that the fiash phase of a 
solar flare may be associated with a thermal instability. For example, explosive 
solutions of the electron energy equation, i.e., solutions which yield an infinite 
temperature in a finite time, are known to occur when collisional Joule dissipation 
dominates the equation. This instability is not relevant for magnctosphcric applica- 
tions or for the upper solaratmosphcre (case II, in table 1). 

in) Beta threshold. It may be hypotliesized [106] that* in a collision-free piasina, 
reconnection is suppressed for liigh |3j values but may occur for small p|. Thus, any 
current sheet in which /3j decreases gradually from some initially large value may be 
converted to a rapidly reconnecting configuration when a critical /?i value is 
reached. In the geomagnetic tail, an abrupt decrease in /?| value occurs if the plasma 
sheet in which the tail current sheet is imbedded shrinks to a ihickiiess equal to the 
current sheet thickness. On the other hand, at the subsolar magnetopause, Lees 
[71] and Zwan and Wolf [132] have described a magnetosheath plasma depletion 
mechanism (by escape along the magnetic field lines) which would tend to maintain 
a value of /Ij of order unity or less. Tlie /?i threshold is not relevant to the double- 
inverse pinch laboratory experiment, and probably not to solar flares because is 
small in these applications (table 1). 

(iii) Current threshold. Assume tiiat a current sheet with little or no reconnection 
gradually thins from an initial width of an ion gyroradius or more toward the 
electron inertia! length, in response to an increased external total pressure, p| + 

2iUO' In this process the current density in the sheet increases gradually. When the 
threshold for onset of current-driven plasma insiabilitics is reached, c.g., for the 
ion-acoustic instability, when cq. (6.1 1) is, satisfied, a reduction in effective 
electrical conductivity takes place in the layer. If this reduction occurs sufficiently 
rapidly, the inductance of the system wiH allow us to consider the current density 
initially to remain essentially unchanged. Instead an inductive electric field ^-(x, f) 
is developed within the sheet to maintain the current density. The magnitude of 
this electric field is larger, the larger the conductivity reduction. 

For a turbulent conductivity of the size Used by Coroniti and Evlatar (22J,£‘^ is 
of the size usually esfimated for steady-state reconnection. This electric field, which 
is initially confined «to the current sheet, is subsequently spread by fast-mode 
expansion waves propagating outward from the sheet as the configuration converts 
itself to one of steady or quasi-steady reconnection. 


If the conductivity is reduced to the level given by cq. (6.l2),£’j. may be one or 
two orders of magnitude larger than typical steady-state values. Smith [102] has 
pointed out that the width 2x* of the current sheet must then increase. As pointed 
out in subsect. 6.2, for reasonable flow rates we find //oOU|.y* « 1 when x* ~ 
Ae and with a given by eq, (6.12). Thus an increase in x* occurs in order to bring 
noovix* toward unity as required in a steady state. The rate -9/?/9r associated with 
(lie increase in x* is the principal source o f E^. But the main resiilLdf the increase in 
X* is that condition (6.1 1) ujtiniatcly is violated so that the ion— acoustic instability 
is quenched. Smith [102] proposes that the process may then repeat itself. A state 
of pulsating reconnection is established. See also Bratcnahl and Baum [II]. While 
the above arguments were given in terms of the ion -acoustic instability, other 
mechanisms may produce similar effects. 

(iv) Tearing threshold. In a collision-frcc current sheet with a vanishing normal 
magnetic field component, electron tearing should be normally present, unless it is 
su pressed by some agent such as pressure anisotropy [18] or velocity shear [60]. 
With a nonvanishiiig normal magnetic field component Bx.m threshold for the onset 
of collision-frcc tearing docs exist, as mentioned in subsect. 6.2. If is originally 
large, no tearing occurs. But as \Bx\ gradually decreases it may set in when the 
gyroperiod in B^ exceeds the growth time. Schindler [99,100] has noted that this 
threshold may be exceeded for ions (but not electrons) in the geomagnetic tail 
current sheet during the thinning of that sheet which occurs in the expansive pliase 
of the geomagnetic substoriti. Since the tail at; this lime has free energy available for 
dissipation, the ultimate result of the onset of ion tearing should be large-scale 
relaxation (via reconnection) of the tail towards a state of minimum free energy, 
rather than merely the generation of tearing turbulence in the sheet. Further 
development of tlie ion-tearing instability theory has been given by Galeev and 
Zclcny [46,47]. 

(v) Interchange instability. An abrupt onset of interchange turbulence in the 
geomagnetic tail [66,89] may occur if the ionosphere becomes decoupled from the 
tail plasnia sheet by the development of electric fields parallel to the magnetic lines 
of force. 

6.4. Particle acceleration 

One of the most important, and at the same time most poorly understood, aspects 
of magnetic-field reconnection is its presumed ability to accelerate particles to high 
energies. Observations in the magnetospheric tail indicate the occurrence of ener- 
getic electron and proton bursts [6,68,116] during times when reconnection may 
be going on. And it should be remembered that curability to observe reconnection 
on the sun and in the far reaches of cosmos depends critically on the generation of 
energetic particles and on the electromagnetic radiation they subsequently produce. 

Particle acceleration may occur either in turbulent small-scale electric fields or in 
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the large-scale reconnection field E^, Both types of acceleration are 

expected to be operative principally in high eurrent diffusion region 

and the shocks. To discuss turbulent acceleration one must understand the nature 
of the dominant micro-proccssc*! in these regions. Since no such understanding is at 
hand, the discussion in this section is confined to particle acceleration in the large- 
scale reconnection electric field. 

In many, but not all, cosmic applications, the total potential difference asso- 
ciated with a steady reconnection electric field is. sufficiently high to account in 
principle for observed particle energies. However, it is only in the small diffusion 
region that particles have the opportunity to move along tlic electric field for any 
considerable distance, And even there, most particles have short residence times and 
undergo a correspondingly small energization, as sitown by cq. (6.8), Thus, steady- 
state reconnection docs not appear to be an effective mechanism foj- tlic accelera- 
tion of particles to very higfi energies f 1 05]. Additionally, in applications such as 
the geomagnetic tail it is necessary to account for particle energies which exceed 
the steady-state cross-tail voltage by an order of magnitude or more. One is there- 
fore led to consider the possibility of particle acceleration during nonstcady recon- 
nection (1 16] . T\yo possible advantages are gained. First, tH inductive electric 
fields may, in principle at least, become ntuch stronger than the quasistatic ones 
during steady reconnection. Second, the nonconservative nature of £* permits 
acceleration within more localised regions of space. For example, betatron accelera- 
tion to high energy may occur in a small region of space where the particles 
experience a large increase in magnetic field intensity.. By contrast cIcctroMatic 
acceleration requires particles to move large distances along the separator. 

The lack of nonsteady reconnection models prevents a detailed analysis of 
particle acceleration. But the simple model given below may serve as an illustration 
of how electron energization might occur in the diffusion region. A resistively 
decaying one-dimensional current sheet, perhaps generated as described in subsect. 
6.3, may be crudely described by 


\yBilx\lx 


\x\<x* 

\x\>x* 


where the sheet width x* is an increasing function of time and Bi is the constant 
field outside the sheet. Assuming no inflow irilo the sheet, the associated electric 
field is 

£ = z^i(\ - (6.14) 

The direction of this field is such that it drives the particles toward the center of 
the current sheet (or = 0). A particle accelerating freely atx = 0 in this electric field 
may be shown to gain an amoun t of energy given by 

A<5 = mc^{[l +(e^iAx*/2mc)^]*^ — 1}, 


(6.15) 
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provided it doesn't leave the system (at z = ±/i) during the time it takes the current 
sheet to widen by Ax*. In lhc geomagnetic tail /?o — 20nt, and for Av* = 500 kiii 
cq. (6.1 5) predicts a possible energy gain of 1 .07 MeV for electrons (protons would 
gain a similar amount of energy only if Ax* 15000 km).Sincc an electron in tins 
energy range traverses the entire tail in less than a second, it would appear that an 
unreasonably large value of dv*/d/ is needed. But this is not neccs.sarily so. If the 
widening current sheet is located along the separator AXB of t!;c rcconncctioti 
bubble in fig. 12, electrons may be accelerated as they move from A to B along the 
separator AXB. They may then return from B to A by gradient drift in the vicinity 
of the O type neutral line BOA where the electric field vanishes. Subsequently they 
reenter the acceleration region at A. By cycling electrons through this loop many 
limes the energy gains predicted by eq. (6.15) may be achieved even for small 
values of dx*/dr. 

The illustrative example discussed above emphasizes that it may be necessary to 
consider three-dimensional time-dependent configurations in order to account for 
particle acceleration in the reconnection process. For further illustrative calcula- 
tions, see refs. [56,1 16]. — 

7. Magnetosplieric evidence 

Much of the observational evidence concerning the possible occurrence of recon- 
nection in the magiictosphcre has sbecn summarized by Burch [12]. Relevant 
references may be found in his paper and are not, for the most part, repeated here. 

A large amount of evidence exists indicating a rclarionship of various magneto- 
spheric activity indices to the southward component of the interplanetary magnetic 
field. Also, spatial asymmetries in a variety of polar-cap processes appear to be 
correlated with the orientation of the interplanetary magnetic field. Such evidence 
is compatible with, but docs not prove, the occurrence of reconnection at the 
magnetopause. This body of observations will not be discussed here. Instead, we 
focus, in subsect. 7.1, on observations relating directly to the transfer of magnetic 
flux from closed to open field lines, and vice versa, in the magnetosphere. If such 
transfer in fact occurs, reconnection of some form must take place. If not, there is 
no need for it. Subsect. 7.7 contains a brief discussion of direct measurements of 
magnetic field and plasma in the vicinity of what may have been reconnection sites. 

7.1. Flux transfer evidence 


The case for the occurrence of flux transfer in the magnetosphere from closed to 
open field lines is based on four sets of observations, discussed below: 

(i) Existence of open field lines in the tail Anderson and Lin [5] have studied the 
shadowing effects on solar electrons (<^ > 20 keV), produced by the moon when It 
is located in the geomagnetic tail. They provide persuasive evidence that a substan- 
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tial amount of magnetic flux iivlhc two tail lobes occurs on open field lines, i.e., on 
lines that intersect the earth’s surface in one place only. But tlic observations do 
not establish how large a fraction of tlic tali flux is on open lines at a given distance 
from the earth. Thus, it is not known for a fact how large a fraction of the earth’s 
polar-cap field lines, i.e., lines emerging at latitudes above the auroral oval, that arc 
open. A popularly held view is that all are open. But llcikkila [53], questioning the 
soundness of _l_his view, has drawn attention to observations niy McDiarniid ct al. 
[77,78] which indicate the common occurrence of trapped particle pitch-angle 
distributions in the day-side cusp region as well as poleward of discrete auroral arcs. 

(ii) Flux erosion from the front-lohe magnetosphere, 'Ihe magnetopause is observed 
to move closer to the earth when the interplanetary field develops a soutliward 
component. At the same time, the dayside polar cusp moves to lower latitudes. 
These effects cannot be accounted for by simple compression of the magneto- 
sphere, Maezawa [79] has estimated that flux on closed field lines is removed from 
the magnetosphere front lobe in an amount estimated at about 10* weberduring a 
typical event. Either this flux is transferred to open field lines in the polar cap by 
dayside magnetopause reconnection or it is moved into tlic tail wliilc remaining on 
dosed field lines. In the latter case, the fiux might be added to the lohe of closed 
field lines in the tail or it might possibly be placed on open field lines by reconnec- 
tion at the tail magnetopause. The popularly held view is that the flux is transferred 
to open flux by reconnection somewhere pn tlic dayside magnetopause. 

Flux addition in the open /a// A siibslanlial body of evidence indicates 
that the magnetic field intensity ia the tail starts to increase shortly after the onset 
of a southward component of the interplanetary magnetic field while at the same 
time the asymptotic tail cross-section increases [80]. The observed concurrent 
gradual thinning of the tail plasma sheet (which Is believed to contain the closed tail 
field lines) argues against these effects being caused by an increase of flux on closed 
field lines in the tail. Rather they indicate an increase of flux on open field lines in 
the two tail lobes by an amount estimated at 1—2.5 X 10* weber. If the auroral oval 
(and the dayside cleft) is associated with the scparatrix between closed and open 
field lines, the motion of this oval to lower latitudes following the southward turn- 
ing of the interplanetary field [62] supports this interpretation. But the evidence, 
while strong, is not conclusive. If dosed field lines occur in the tail outside (i.e., 
above and below) tlie plasma sheet, the flux on open field lines could conceivably 
remain unchanged. 

(iv) Polar cap electric fields and convection, measurements [55| in the 

polar ionosphere indicate an average Voltage difference across the polar cap of the 
order of 65* kV, corresponding to a magnetic flux transport across the cap from 
dayside to nightside a« a rate of about 2 X 10* weber/h [80], Ion flow measure- 
ments [52] over the polar cap show fiow patterns that carry particles and, unless 
£ magnetic flux poleward across the dayside cleft in a narrow longitude 
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sector. Most but not all of the flmv in the cap region occurs near the cquatorward 
edge of the cap adjacciit to an abrupt flow reversal, below which the return flow to 
the dayside occurs. While the exact location of the scparatrix between closed and 
open field is not known, it is difficult to locate it in such a way that these results do 
not imply a transfer of flux from closed to open field lines. 

in spite of the ambiguities in the interpretation of the observations listed above, 
their mutual consistency in terms of flux transfer rales is impressive and lends 
credence to the idea that flux transfer from closed to open field lines does occur in 
the magnetosphere. However, a far greater body of simultaneous observations by 
satellites at different locations in the magnetosphere needs to be examined in order 
to establish the validity of the idea in a conclusive manner. It is noted that on the 
average, any flux transfer from closed to open field lines must be balanced by a 
reverse transfer from open to closed lines. Tail reconnection, occurring sporadically 
in connection witli the expansive phase of magnetic substorms, is thought to 
accomplish this latter transfer but conclusive evidence is not available (see subsect. 
7.2). 

7,2.AIcasurimicnts near reconnection sites 

In a strict sense, direct evidence for reconnection Consists of in situ observations of 
the hyperbolic magnetic field configuration associated with a separator and an 
electric field along that tine. The electric field observation may convincingly be 
replaced by tlie observation of plasma energized in the reconnection process 
(compare subsects. 2.4 and 4.2). 

Hones et al. [63] have reported observations of proton fluxes and of magnetic 
fields in the geomagnetic tail at geocentric distances in the range of 25-32 earth 
radii. They have found substorm events in which tailward proton flows at speeds up 
to lOOOkni/sec and an associated southward component of the magnetic field 
occurred during the storm expansion phase, followed by earthward flow and a 
nortliward field component during recovery. Such observations arc consistent with 
a separator moving tailward past tlie satellite. However, evidence concerning the 
magnctic-ficId component perpendicular to the tail current sheet is not entirely 
convincing unless tlic field is measured near the center of the sheet, which w'as not 
tlic case. And recently Lui et al. (74,75] have challcugcd observations purporting to 
show the formation of a near earth reconnection line during substorms. Observa- 
tions of proton jctiing in the tail [39], of energetic particle bursts [6,68,98], and 
of lunar shadow patterns of electron fluxes [73], while generally compatible with 
tail reconnection, nevertheless cannot be claimed to provide unambiguous proof of 
the occurrence of the process. 

At the dayside magnetopause, magnetic field components perpendicular to the 
magnetopause have been observed [108], although not as a permanent feature not 
even when the magnetosheath field opposes the terrestrial one. The narrow jets of 
energized plasma, predicted by magnetopause reconnection models (e.g., fig. 19) 
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and flowing nearly tangential to ihc magnetopause, have no been seen *, even 
thougli satellites such as HEOS 2 Iiave had tlie riglit position and attitude to observe 
them [50,5 1 1. These facts along with recent observations of a plasma boundary 
layer inside the daysidc magnetopause [29,35] suggest that magnetopause recon- 
nection, if it occurs, may be more sporadic and more localized than originally 
expected. Furthermore, the possibility of reconnection in the cusps and elsewhere 
on the magnetopause surface, rather tlian near the sub-solar point, needs to be 
examined [50]. 

The absence of observations of plasma energized by daysidc reconnection has led 
Heikkila [54] to suggest that no such reconnection occurs, i.e,, that the magneto- 
pause is an electrostatic cquipotential. This suggestion is difficult to reconcile with 
the presence of magnetic field components perpendicular to the magnetopause, 
unless one is willing to accept potential differences of the order of 50 kV along 
field lines extending from the magnetopause into the solar wind; or unless one 
argues that such perpendicular components are never present over any substantial 
part of the daysidc magnetopause. 

8. Summary and recommendations 

In this paper we have given a reasonably detailed description of the present status 
of our understanding of reconnection. The picture that emerges is of a process, 
siniple in concept but extremely complicated and multifaceted in detail. Nonlinear 
magnetohydrodynamic processes in the external flow region, governed by distant 
boundary conditions, are coupled to non-linear microscopic plasma processes in the 
diffusion region in a manner not clearly understood. And It appears that reconnec- 
tion may operate in entirely different ways for different plasma parameters and for 
different external boundary conditions. Steady reconnection may be allowed in 
some cases, forbidden in others, with intermediate situations involving impulsive or 
pulsative events. 

On the whole, our theoretical and empirical knowledge of reconnection is poor. 
Yet the process plays a key role in solar-flare theory as well as in our present con- 
cept of the dynamic magnetosphere. And it appears as an unwanted feature in 
tokamaks and other fusion configurations. These facts, along with the potential 
importance of reconnection in other parts of the cosmos, amply justify vigorous 
research efforts related to reconnection in the following five areas: solar-flare and 
astrophysical observations, magnetospheric observations, laboratory experiments, 
computer simulation, and analytical model building. The first area, while extremely 
important, is too broad to be commented upon here. In the remaining areas the 
following recornmendations are made: 

* However, a layer of energetic electrons of unknown origin has been discovered outside the tail 
magnetopause [7,81]. 
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{\) Magnetospheric observatiom and experiments s\\ovAA \nc\\xdQ: 

(1) A coordinated program to establish (or deny) the occurrence of flux transfer 
across separatrix surfaces, and to study other global consequences of reconnection; 

(2) Direct observations of magnetic field, plasma, energetic particles, and fluctuat- 
ing as well as steady electric fields, near magnetospheric reconnection sites. Multi- 
salellilc missions arc needed to separate spatial and temporal effects; 

(3) Perhaps active experiments, such as the release of barium clouds near reconnec- 
tion sites. 

(ii) Laboratory experiments. The observation of impulsive flux transfer svenU aild 
of ion-acoustic turbulence in the double inverse pinch experiment illustrates the 
importance of such experiments in shaping our understanding of reconnection. Yet, 
(excluding fusion devices) the double inverse pinch appears to be the only operatirtg 
reconnection cxpcrimcnl in the US today. A substantially expanded laboratory 
program is needed with four principal goals: 

(1) Simulation of solar-flare reconnection; 

(2) Simulation of niagnci ^spheric reconnection; 

(3) Study of basic plasma processes of importhnee in rcconneciion. such as slow 
shocks and anomalous resistivity; 

(4) Exploration of reconnection in plasma heating devices. 

(iii) Computer simulation provides a potentially very powerful tool for the study of 
reconnection. Magnctohydrodynamic codes, and ultimately self-consistent particle- 
fields codes should be developed. It is particularly important to build into such 
simulations the effects of Inertial and anomalous resistivity in the diffusion region. 

(jv) Analytical models of reconnection should emphasize the following interrelated 
problems: 

(1) Nonsteady and three-dimensional effects; 

(2) Plasma processes in the diffusion region; 

(3) Particle acceleration; 

(4) Reconnection of fields that are not antiparallel. 

It Is throUgl) vigorous activitiesun the aforementioned areas, and effective inter- 
action between scientists involved in them, that our understanding of the reconnec- 
tion process may be most rapidly advanced. To bring about such a state of affairs, 
two proposals arc made: 

(A) That a special working group be assembled with the charge of promoting 
effective research on all aspects of the reconnection problem and with membersliip 
drawn from the five research areas discussed above. 

(B) That NASA and other funding agencies develop coordinated programs of 
support for reconnection research. 

The importance of the reconnection concept is such that we can ill afford the 
present somewhat haphazard approach to its study. 
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Note added in proof 

Several developments relating to reconnection have occurred after the comple- 
tion of this paper. New computer models have been reported f 133,137,139,144], 
Laboratory experiments [142] simulating the flow of magnetized plasmrx past the 
magnetosphere have revealed magnetic field configurations suggestive of reconnec- 
tion in the tail as well as at the dayside magnetopause (sometimes at the sub-solar 
point, sometimes in tlie cusps). Space experiments [1 40] have revealed the presence 
of a substantial electric fieldj tangential to the dayside magnetopause. An analytical 
model of resistive current layer decay has bfeen developed [138], Mathematical 
studies of plasma motion near magnetic null points [135,143] have been reported. 
And an exact compressible MHD mpdel of the convection region flow, based on fig. 
14, has been developed [I36J, It contains intermediate, rather than fast or slow 
waves in the inflow regions. Finally, the ion tearing mode hasbeen studied further 
[134] and its existence has been called into question [141], 
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ABSTRACT 

Mass, momentum, and energy transferred from the solar wind and 
magneto sheath into the magnetosphere must cross the boundary of the 
earth's magnetic field, the magnetopause. This paper reviews the pos- 
sibilities for convective and diffusive transport across the boundary 
and discusses the dependence of such processes on the orientation of 
the interplanetary magnetic field. It is first shown that, in the 
absence of transport processes, the magnetosheath magnetic field just 
outside the magnetopause always builds up to a dynamically significant 
level. The possibility of macroscopic plasma flow (convection) across 
the magnetopause is then discussed both for the case of a vanishing 
and a nonvanishing magnetic field component normal to the boundary. 

In the former case, convective flow does not appear possible, unless 
substantial electric fields occur parallel to the magnetic field. The 
latter case includes exactly field-aligned flows as well as magnetic- 
field reconnection. Particular attention is given to the problem of 
particle energization during reconnection. Finally, a brief discussion 
is given of diffusion processes and the constraints placed upon them 
by existing observations. It is shown that the efficiency of diffusive 
as well as convective transfer should be expected to have a strong 
dependence on the angle 6 between the magnetospheric and the magneto- 
sheath magnetic fields with a maximum for 0 = IT, a minimum for 0 = 0. 

I. INTRODUCTION 

All mass, momentum, and energy transferred from the solar wind to 
the magnetosphere must cross the magnetopause. The processes whereby 
such transfer takes place are poorly understood at present, and they 
constitute one of the roost important unsolved problems in magnetospheric 
physics. It is the purpose of this paper to provide a brief qualita- 
tive review of possible transfer mechanisms and the extent to which 
they are able to explain plasma and field observations near the day- 
side magnetopause. The nature of these observations has been discussed 
recently by Eastman and Hones (1979), Paschmann et al. (1979), and by 
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Russell and Elphic (1979). A study of the magnetic structure of the 
magnetopause may be found in Sonnerup (1976). In the present paper, 
particular attention will be paid to the question of how the transfer 
rates may be modulated by the orientation of the interplanetary mag- 
netic field. 


Traditionally ^ transfer processes are classified in three groups; j 

(i) Transport by convection, i.e., in the present application, 

by direct macroscopic flow of .plasma across the magnetopause. ? { 

. ■ r‘l 

. i 

(ii) Transport by diffusion, i.e., by processes that are con- 
sidered microscopic on the temporal and spatial scales of interest. fj 

■ . • . ■■ ' Li 

(iii) Transport.by radiation, i.e., by waves propagating across 

the magnetopause. ■ ^ 

The last of these categories is rather special because it involves 
no mass transfer. Apart from the electromagnetic radiation, there may 
also be momentum and energy transfer across the magnetopause via com- .^ | 
pressive MHD waves, as pointed out by Axford (1964). This topic and - * 
the related one of the Kelvin-Helmholtz instability, while important, 
are beyond the scope of the present paper. - | 


The distinction between convection and diffusion, depends on one's 
definition of the terms macroscopic and microscopic. In other words, r| 

it depends on the space and time resolution available (or desired). -f r 

Turbulent transport may be considered a diffusive process, described 
by an effective "eddy" diffusivity, if one is concerned only with time 
and space scales long compared to the correlation times and correlation j | 
lengths of the turbulence. With higher temporal and spatial resolution,- - 
it should be properly considered a convective process. For the pur- 
poses of this paper, time scales shorter than a few proton gyroperiods ; 

and length scales shorter than a few proton gyroradii must be considered. { 
microscopic. (The gyroperiod and gyro radius of a IkV proton in a 50nt 
field are 1.3 sec and 92 km, respectively.) Hi' 

i \ \ 

5 - :-S ^ ; 

The paper is organized as follows. First, a discussion is giyen 
of the "ideal" or nontransfer state of the magnetosphere with focus j 

on the magnetic field near the magnetopause. This is an appropriate ■: } 1 

starting point, for the transport across the magnetopause usually ap- " i ' • 

pears to be sufficiently modest so as to be considered a small pertur- 
bation on the nontransfer state. Furthermore, well established cor- fl | 

relations between various geomagnetic effects and the interplanetary Ll i 

magnetic field (e.g., Burch, 1974), indicate the importance of under- 
standing the nature of the magnetosheath magnetic field immediately Tj 
outside the magnetopause. - j 
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Second, a presentation is given of existing convective transfer 
mechanisms by which plasma flows directly across the magnetopause, 
usually under the influence of an electric field parallel to that 
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surface. The cases where the magnetic field component normal to the 
magnetopause vanishes and does not vanish are considered separately. 

, Third, a brief suimnary is given of diffusive processes that may 
bd operative at the magnetopause and the constraints placed on them by 
existing observations of the plasma boundary layer immediately inside 
the magnetopause. 
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2, THE NONTRANSFER STATE 

' i ■ 

j The problem of ideal flow past the magnetosphere was dealt with 
lii great detail in the sixties (e.g., Spreiter and Alksne, 1969). In 
general, the approach was to calculate the magnetopause shape and the 

• flow configuration under the assumption that the interplanetary mag- 
netic field plays no dynamic role other than that of rendering the 
pljasma a continuum. Tlie magnetic field configuration outside the mag- 
netopause was calculated afterwards (Alksne and Webster, 1970) by use 
of the frozen field condition. The point to be made here is that, no 
matter how weak the interplanetary field, such a procedure always leads 
to a violation of the basic assumption that the Maxwell stresses may be 
ignored. We illustrate this point by considering ideal steady axisym- 
metric supersonic flow past a blunt-nosed impenetrable and perfectly 
diamagnetic object, as shown in Figure 1. In each roeridional plane, 
the streamlines in such a flow are given by \i»(R,0) = const., where 
V>(R,0) is Stokes’ stream function and the notation in the figure is 
used.^ 



Figure 1. Ideal IIHD flow past an impenetrable diamagnetic object L 
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The velocity v Is given by 

V *= (Vijjx^)/pRsln0 , (1) 

\ f • 

p being the density and ^ a unit vector in the azimuthal direction. 

The frozen magnetic field condition 

E + V X JB = 0 (2) 

with 

E = - V $ (3) 

ijmplieci that not only the magnetic field lines but also the streamlines 
are equipotentials. Since each streamline lies along the intersection 
b|etween a meridional plane, (J) = const., and an axisymmetric stream 
surface, = const., the electrostatic potential may be written 
$ = . The functional form of this latter expression may be eval- 
uated in the region upstream of the bow shock (subscript: where the 

flow and field are uniform. With ^ say, we find 

. . ^CO = = -v^B^Rsiia0sin(j) (4) 

But the stream function for uniform flow with speed is 

'L = -^Pc„v„R''sin^0 . (5) 

so that ■> 

^oo = (6) 

One concludes from Equation (6) that the potentisil distribution in the 
entire flow field must be given by 

4>(R,0,<J>) = -By-2v^\|j/p„ sintl) ^ ^ (7) 

By use of the frozen field condition, the component of the magnetic 
field perpendicular to the flow velocity vector may now be written as 

= JB - v(v*g)/v^ - V$ X (8) 

■ I , ^ , , - ^ , ■ 

Tire component of JB along v is not of interest here, but it may be ob- 
tained from ^ ^ = 0 . 

i The stream function vanishes on the axis (see Equation (5) for 
0j= 0) and since the streamline along the axis splits at the forward 
stagnation point and then covers the entire surface of the object, we 
conclude that \l> vanishes on that, surface. In calculating W, 
given by Equation 7, a factor- will be generated. Thus when W is 
substituted into Equation .(8) , ^^becomes infinite at “ 0 , i.e^ 
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thie entire surface of the blunt object (Alksne and Webster, 1970). The 
only exception occurs if p = 0 on the surface; this possibility will be 
discussed below. 

I As an example, the magnetic field intensity on the axis (0 = 0) 
has been calculated as a function of position using Lighthill’s well- 
known constant-density solution for hypersonic flow past a sphere (see 
Hayes and Probstein, 1959) for which the stream function in the region 
between the shock and the sphere is given by 

'P “ [PcoV„(R3sin0)V(3Oe')][3(l-e)^R/Rg)'* - 

5(l-4e)(R/Rj^ + 2(l-e)(l-6E)(R /R)] (9) 

, s s 

Here R is the shock radius (the radius of the sphere, Rj^j, is obtained 
by solving the equation ip - 0) and 

£ = (y-i)/(y+i) (10) 

Y being the ratio of specific heats at constant pressure and constant 
volume. The resulting Bj_ field distribution along the stagnation 
streamline is shown in Figure 2. 



Figure 2. Magnetic field as a function of position R/Rj^j on the stag- 
nation streamline (the x axis) for Ligh thill Vs solution with Y = 2. 
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One may understand the infinite value of the surface magnetic 
field in the following manner. Consider a plasma torus of major and 
minor radius a and respectively, moving towards the sphere as shown 
in Figure 1. As it approaches the stagnation point it must expand its 
circumference 2ira enormously in order to allow the object to pass 
through. In this process the mass of the torus is preserved so that 

p(2TTa)(Tf2)*) - cons.t. (11) 

i , 

At the same time, the magnetic flux trapped within the cross sectional 
area n£»^ is also conserved: 

= const. (12) 

The ratio of these two expressions is 

Bj_/(pa) = const. (13) 

Unless p decreases proportionately to the increase in a as the torus 
expands past the object, one concludes that must increase. In the 

limit of an infinitely small torus threading the stagnation streamline 
the fractional increase in a and hence in becomes infinite as the 
torus expands and moves along the surface of the object. 

One concludes from the preceding result that in ideal MHD flow the 
original assumption of negligible Maxwell stresses can never be true in 
;the magnetosheath flow near the magnetopause. The net effect of the 
action of these stresses must be either to decrease the density to zero 
at the surface of the object, as suggested by Lees (1964) and further 
developed by Zwan and Wolf (1976), or to change the flow topology in 
such a way that Bj_ may remain finite. A substantial decrease in den- 
sity immediately outside the magnetopause has in fact been seen 
(Crooker and Siscoe, 1975; Paschmann et al., 1979). Nevertheless, the 
density does not go to zero so that the dilemma remains. 

An example of a topological change is shown in Figure 3. The flow 
ceases to be axisymmetric and a stagnation line, parallel to the local 
magnetic field, forms on the nose of the object. At the two ends of 
the stagnation line the magnetic field has null points from which mag- 
netic null lines emerge. 

The formation of stagnation lines aligned with JB and magnetic null 
lines aligned with v may be predicted directly from the frozen magnetic 
field condition in the form 

(B/p) • W = V • V(B/p) (14) 

If \r = 0 and jB # 0, p 0, at a certain point on the surface of the dia- 
magnetic object then B • Vv = 0 there. Thus remains zero if one 
moves along the surface in the direction of B. Similarly, if ^ - 0 
and V 0, p 0, at some point, then v • V(^/p) = 0 there. Thus ^ 




Figure 3- Magnetosheath field and flow on the subsolar magnetopause. 
Streamlines are dotted, magnetic field lines solid. The line segment 
Bj-C is a stagnation line at which each incoming streamline splits. The 
line segments A~B and C-D are magnetic null lines at which each in-* 
coming magnetic field line splits. The points B and C have JB = 0 and 
V - d. 


remains zero if one moves along v. It may be predicted qualitatively 
that the length of the stagnation line increases with decreasing 3^ 
value, 3^ being the ratio of plasma to magnetic pressure in the solur 
wind. 


The flow near the stagnation region in Figure 3 is organized to 
be predominantly perpendicular to the magnetic field. The main result 
of this geometry is that the torus-shaped plasma element in the pre- 
ceding discussion must be replaced by a chain-link shaped one which 
threads all of the magnetosheath stream lines joining the stagnation 
line. The fractional increase in circumference of such an element is 
finite as it expands to allow passage of the object through it. Thus 
the cross section of the chain-link will remain finite and, with a 
constant magnetic flux trapped in that cross section, the magnetic 
field will also remain finite. 

' The implications of ideal MHD flow past a diamagnetic object are 
described above. In the real case of solar-wind flow past the mag- 
netosphere, deviations from the frozen-field condition at the magneto- 
pause may obscure the effects shown in Figure 3 either partially or 
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cbmpletely. If in reality there exist conditions such that the front 
side of the magnetosphere remains essentially impenetrable to the mag- 
netosheath field and plasma (the "ground-state" magnetosphere) then 
some of the features shown in the figure should be releA^ant. On the 
other hand, if field line interconnection across the magnetopause is 
important, the field topology will be quite different. 


1 The basic prediction of the MHD theory for flow past an impene- 
trable object is that the magnetic field immediately outside the mag- 
netopause builds up until it becomes dynamically significant. Poten- 
tially it therefore becomes sufficiently strong to exercise a control- 
ling influence on the transfer processes at the magnetopause. Without 
siich processes a 3 value of the plasma of order unity or less should 
be expected in the stagnation region just outside the magnetopause. 

Ori the other hand, observations of large 3 values would be a clear in- 
dication of a substantial violation of the frozen magnetic field con- 
dition, presumably in association with strong transfer across the mag- 
netopause. Both small and large 3 values have in fact been reported 
(see Paschmann et al., 1979). It is of interest to note that the 
smallest value occurred when the magnetosheath field was due north, 
suggesting that the transfer may be a minimum for this field direction. 


3., CONVECTIVE TRANSFER 


In order to discuss the possibility of convective transfer across 
the magnetopause we adopt a simple one-dimensional time- independent 
model of the local magnetopause structure. In such a model. Maxwell's 
equations require the normal component of the magnetic field and the 
tangential component of the electric field to remain constant across 
the layer. .Denoting the magnetosheath and magnetosphere sides by sub- 
scripts 1 and 2, respectively, we find 


B = B • n = B • n = B = B 
ni ~i ~~2 ~ nz n 

E = E. = E. 

— ti — 12 —t 


(15) 

(16) 


where the subscripts n and t refer to components normal and tangential 
to the magnetopause surface. The unit vector ^ is perpendicular to 
that surface and points away from the earth. 

We shall explore the consequences of the assumption that the elec^ 
trie field has no component along the magnetic field, i.e.. 


E 

• B 

= E, 


+ 

E B =0 


“"1 

— t 

-tl 


ni n 


• B 

= E^ 

* B .. 

+ 

E B = 0 

—2 

*~2 

— t 

— 12 


nz n 


(17) 


For given E and the normal flow velocity on either side of the mag- 
netopause may be written 
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(E^XB)/B^ + n • V)| (18) 

wherein is the velocity along 

■ i 

The model described above does not require the actual magnetopause 
to be stationary or smooth ^ It may be applied locally to a wavy and 
rapidly moving and deforming current layer as long as the frozen mag- 
netic field condition holds on the two sides. (This condition does 
not have to be, and generally is not, valid within the magnetopause 
structure itself.) However, when used to interpret satellite magnetic 
and electric field observations during magnetopause crossings, the pos- 
sibility of time aliasing must be considered. 

Observations of the magnetopause (Sonnerup, 1976) indicate the 
importance of examining both the case = 0 and the case ^ 0. 


4. 


THE CASE B = 0 
n 


We now examine the possibility of plasma convection across. the 
magnetopause when the normal magnetic field component is vanishingly 
small. It then follows from Equation (17) that E must be perpendicu- 
lar to » ss well as to B^’ = . For t 0 this is clearly 

possible ^only^if ^ and are^eithir parallel or antiparallel, i.e. 
if where R^is a positive or negative factor. Thus we have 

the following cases: 


(i) B M 


= 0, in which case v 


Either E^ 

‘1 — 2 — t 

that no convection occurs, or there is an electric 


= V 


fiSld 


so 


- 0 

component par- 
allel to the magnetic field on at least one side of the magnetopause, 
in general, this parallel electric field would have a magnitude compa- 
rable to the field It is questionable whether such a field can be 

sustained in a quasisteady state and on spatial scales sufficiently 
large to be classified as macroscopic. 


(ii) B = kB ; k < 0. The magnetic fields on the two sides of 

the magnetopause are antiparallel. The case k = -1 has been discussed 
by Alfven (1968) and by Cowley (1973). The prediction is that a vol- 
tage difference = B^/(y Ne) can be sustained along the magnetopause 
surface (and perpendicular ^ to the magnetic field). For B = 50 nt, 

N = lO^m we find A$ = 1.24 kV. This result should be compared to 
the typical value of 50 kV for the observed potential difference across 
the polar cap. It would appear that Alfven’ s mechanism cannot explain 
the observations. ' 

(iii) B = kB ; 0 < k < 1. If ^ 0 then the magnetopause must 

be identified with a fast MHD shock. Since this is certain not to be 
the case, we conclude that - 0 so that plasina convection across the 
magnetopause does not occur. 


- 10 - 


' (iv) B = kB ; k > 1. The magnetopause would then be a reverse 
fast shock,* which does not exist. Again, no plasma CGnYectlon across 
the magnetopause is possible. 

(v) Bj = B 2 . In this case, no restriction occurs on the tangen- 
tial electric field except that it be perpendicular to the magnetic 
field. From a local viewpoint, an arbitrarily large plasma flow may 
occur across the magnetopause. This is the case assumed by Alfven 
(e.g., 1958) in his theory of the geomagnetic storm. In this situation 
the magnetopause cannot be identified locally as a current sheet. 
Rather, it is defined in a global sense by the fact that the field 
lines on its two sides have different topological origins , On the 
magnetopause surface itself, the field lines converge to the two mag- 
netically neutral lines, shown in Figure 3. Unless E * ^ ^ 0, this 
topological feature will lead to a short circuiting along ^ of any 
transverse electric field in the equatorial plane. Since magnetde sub- 
storm activity correlates positively with the southward, not the north- 
ward, component of the interplanetary magnetic field, it appears that 
the short circuiting along the field lines may be remarkably efficient, 
a fact that needs to be taken into account in local models of plasma 
flow across the magnetopause (e. g. , Cole, 1974; Formisano et al., 1978; 
Lemaire and Roth, 1978). 

In summary, it appears that, unless the condition ^ ^ “ 0 is 

violated, the case B - 0 leads to the classic signatures of a tangen- 
tial discontinuity and does not offer any hope for convective plasma 
entry across the magnetopause. And abandonment of ^ ^ not a 

Step that can be taken lightly. To be sure, the electric field along 
jB is not identically equal to zero. Voltage differences comparable to 
the average particle energy may be easily sustained along the magnetic 
field lines. But the proton energy is of the order of IkV or less, 
which is insignificant compared to a 50kV potential difference across 
the polar cap. Various reconnection geometries, to be discussed in 
the next section, bypass this difficulty by admitting lar^e potential 
drops only along singular field lines where the X-type field topology 
permits of such voltage drops without the generation of enormous cur- 
rents. 


5. THE CASE B 0 
. n 

When the normal magnetic field B is finite. Equations (17) do 
not place a constraint on the value of the tangential electric field 
Rather, for given the two equations may be used to solye for 
B and E , the normal electric field components on the two sides of 
magnetopause. The plasma convection speed across the magnetopause 
is given by Equation (18) which shows that the parallel velocity may 
now yield a convective flow across the magnetopause even when - 0. 

In such a case we have interconnection of the field lines across the 
magnetopause but no reconnection. When 0 reconnection pccurs. 

Both cases are discussed below. 
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Equatlons (17) and (18) describe the kinematics of the flow across 
the magnetopause. The dynamics consists of the balancing of the Max- 
well stresses (the ^ x ^ force, being the Chapman-Ferraro current per 
unit width) against the change in particle momentum as the plasma 
crosses the current layer. It is dealt with most conveniently by 
transforming to a local coordinate system S' which moves along the 
magnetopause at a velocity v such that the electric fields ^ and ^ 
are transformed away. The relationship between v^ and is; 

E^ + V X = 0 . (19) 

— t — 0 — n 


For the case of antiparallel fields of equal magnitude, B^, (^^ = ) 

see Figure 4a, the momentum balance becomes ^ ^ 


2B B^/y = -2pv v' 
n t 0 n t 


( 20 ) 


where ± v' is the tangential flow speed in the frame S'. By use of the 
trigonometric relations 


B = Bsin6 
,n 

B^ = Bcos6 


1 V = -v' sin6 

n II 1 

J * ^|l j 


( 21 , 22 ) 


which follow directly from the geometry in Figure 4a and in which v*| 
is the flow speed along ^ in the frame S', one finds 


II 




(23) 


In other words, the flow speed along ^ is equal to the Alfven speed 
and the current layer itself should be identified as a large amplitude 
Alfven wave or rotational discontinuity. A more general analysis which 
permits of arbitrary orientations and magnitudes of and , and 
associated nonisotropic and unequal pressure tensors aAd unequal den- 
sities on the two sides of the magnetopause, may be found in Hudson 
(1970). The result is that the flow speed on either side is equal to 
the modified Alfven speed v^ ^ For simplicity 

we pursue only the case of isotropic pressure and antiparallel fields 
here. 


We now return to the rnagnetospheric frame of reference by means 

of the transformation velocity v which we first assume to lie in the 

plane of Figure 4a so that the ^electric field is perpendicular to 

that plane. Expressing the transformation speed as a fraction of the 

Alfven speed, v E rv. , we then find 
0 A 

= rv.B (24) 

t An * 

For small values of the angle 6 in Figure 4a the downward flow veloci- 
ties are 



( 27 ) 


fAS = (2rB^/U N)(Nv.B /B) = (2B/u )E, = • I 

0 An 0 t “t — 

In other words, the power supplied to the particles per unit area of 
the magnetopause is exactly equal to the rate of dissipation of elec- 
tromagnetic energy per unit area. Thus, the energization is exactly 
the one discussed by Heikkila (1975,1978) by use of Poynting's theorem.' 

We examine the meaning of Equations (24)-(27) for different r 
values : . 

(i) r < 0. In this case, • I < 0 and the particles lose 
energy in crossing the magnetopause. This case has been occasionally 
invoked for the front lobe of the magnetosphere with the reconnection 
site in one of the polar cusps rather than in the equatorial plane. 

The electric field is then opposite to the interplanetary electric 
field. *; 

(ii) r = 0. In this case = 0 and no particle energization 
occurs’. The plasma crosses the magnetopause by flowing exabtly along'; 
the magnetic lines of force with the Alfven speed. It is the case ofi 
Interconnection of field lines without reconnection. Such flows could 
conceivably occur on parts of the magnetospheric front lobe. Flow ve- 
locity reversals across the magnetopause have been seen (Paschmann et 
al., 1979) but it has not been established whether cases occur where 
the flow is exactly field aligned and Alfvenic. 

(iii) 0 < r < 1. The electric field is now parallel to the 
Chapman-Ferraro current Particles are energized and reconnection 
occurs. A reversal of the tangential component of the flow still occurs 
at,' the magnetopause. 

(iv) r = 1. In this case the plasma inflow' is perpendicular to 
the magnetopause. For N = lO^m B - 50nt we find A£ = 2,5 kV, while 
for N = 3 X 10^m~^, B = 35 nt, A£ = 400 V. The energization per parti- 
cle is substantial and it may be argued (Heikkila, 1975) that it should 
be easily observable. 

• 

(v) r >1. In this case no flow reversal occurs across the mag- 
netopause. Tlie particle energization is greater than in (iv) by a fac- - 
tor r. 


. It is an interesting fact that the energization given by Equation 
(26) is independent of the magnitude of the electric field E. , even 
though the mechanism of energization must be the displacement of the 
particles along the electric fields in the magnetopause. The explana- 
tion is as follows. Equation (24) shows that small values of E cor- 
respond to small values of B . But for small B values the drift dis- 
placement of individual particles along E is large. As E^ increases, 
B must also increase proportionately. The result is a smaller drift 
displacement in a larger electric field such that the total voltage 
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change experienced by a particle is independent of both and 
What does occur as E and B increase is that the inflow rate, f, of 
plasma increases with increasing E^ so that more particles are energized 
per unit time and unit area. 

Tlie derivation , of Equation (26) emphasizes that the particles ener- 
gized are the particles responsible for niaintiiining the momentum bal- 
ance in the current layer. ' For example, if this balance were main- 
tained by a small number of escaping radiation belt particles then the 
density N in Equation (26) should be the sntall density of these part i- 
fcles. The energization per particle would be correspondingly large. 
However, in MHD theories of reconnection it is the magnetosheath ions 
which maintain the momentum balance and are energised'. The energy gain 
of electrons is smaller by the mass ratio m^/m^. it is here that the 
reconnection concept, applied at the magnetopause, faces its most for- 
midable observational test. Three years ago, Helkkila (1975) poinced 
out that no observations available at that time indicated the presence 
of such energized protons just inside the magnetopause. More .recent 
Studies of the plasma boundary layer (Bserendel et al. , 1978; Eastmap 
and Hones, 1979; Paschniatui et al. , 1979) also haVe failed to reveal 
convincing evidence of proton energization. On the other hand, Mozer 
Ct al., (1979) have reported direct measurements of substantial tan- 
gential fields JEj.. It is also noted that an energetic electron layer 
bas been found outside the magnetopause (Meng and Anderson, 1970,1975; 
Baker and Stone, 1977,1978). However, it is difficult to construct a 
model in which those particles : are the ones responsible for the iiKjmen- 
tum balance in the magnetbpause. 

Helkkila (1978) has examined a large number of ways of accounting 
for the discrepancy between MHD reconnection theory and observations. 

He found them all wanting and concluded that recoiinection on the day- 
side magnetopause docs not occur. Yet, there are compelling. reasons, 
in particular the observed flux erosion of the magnetospheric front 
lobe when the interplanetary field has a southward component (e. g. , 

Aubry et al., 1970; Holzer and Slavin, 1978,1979), xdiich make it dif- 
ficult to discard reconnection. Below we examine several ways of 
avoiding Helkkila ’s conclusion; 

(a) Equation (26) represents an upper limit on the energization 
in two ways. First, it represents a maximum, because was chosen to 
be parallel to the net current _I in the magiietopause so that the elec- 
tromagnetic dissipation rate _E^ • ^ is a maximum. Although no conclu- 
sive proof is available, there'*arc strong theoretical arguments (Cowley, 
1974,1976) indicating that other directions of may be possible. For 
the case = ~”t 2 ’ ~t even be perpendicular to I so that the dis- 
sipation rate is zero. The transfotmatlon velocity v ii; Figure 4a is 

then perpendicular to, the plane of the figure and the particle velocity 
is equal to (v^ + v^)’^ on both sides of the magnetbpause. A second 

point is that the energization is less when the two fields J4. and 1)^.^ 
ate not antiparallel. Figure 4b shows a conib ination of the E>io effects. 
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It also suggests the range of, directions of the electric ‘ field for 
which reconnedtion-like flows occur, i.e., flows in which perpendicular 
magnetic flux is transported away from the reconnection line (which 
has-been tacitly assumed to be parallel to E ) . The calculation of the 
particle energization for an arbitrary configuration is a straightfor- 
ward geometrical task. 


Except for the extreme caseEj. = mentioned above, there is 

no possibility ■ of removing all of tne particle acceleration associated 
with dayside reconnection in this manner. But the two effects dis- 
cussed here may help reduce the energization to the point where it 
could have been overlooked. 


j (b) If the magnetopause is partly reflective,'' energized parti- 
clies will appear on both sides of the current layer. Thus the full 
energization is not discovered by comparing particle energies inside 
and outside the magnetopause. Rather it becomes necessary to separate 
the particles on the magnetosheath side ‘into incident and reflected 
populations before the enex-gization ca:i be established. It should be 
added that such a model of the magnetopaose would contain an abx-upt 
density change so that the internal boundary layer has significantly 
lower density than the magnetosheath. A completely reflective magneto- 
pause may be contemplated too (Sonnerup, 1976) but momentum conserva- 
tion then requires the two fields and B^^ to be exactly antiparal- 
lel (or possibly exactly parallel). 


Again, reflection effects do not remove the energization of parti- 
cles. But they do provide a set of circumstances in which it may have 
been overlooked. 


(c) The energy gained by the ions in the reconnection electric 
field may be transferred to the electrons, or possibly to plasnua waves, 
by some internal iiiagnetopause process. The nature of this process is 
unknown. Hov^^ever, as a simple example, assume that the magnetopause 
structure, viexjed in the frame S', contains a normal electric field 
E' . Equation (20) assumes no such field to be present. If E* is 
dxrected outward fx'om the earth, then, in the frame S', the ions will 
be decelerated, the electrons accelerated, as they cross the magneto- 
pause, while the energies of reflected particles wil]. x'cmain unchanged. 
Since the potential increase across the magnetopause cannot exceed the 
incident energy of the transmitted ions it is easy to show that their 
exit energy in the magnetospheric frame can be less than that given by 
Equation (26) by a factor of at most 4r/ (2x*-l) v However , if part of 
tjie incident ion population a.s refleGted, a situation may now arise 
where the average ion energy outside the magnetopause is greater than 
that inside. Energized electrons will appear in a layei' inside the 
miignetopause. It is not clear in what Circumstances an electric field 
E^^ could occur self-consistently in the magnetopause. One requirement 
Would be that in the absence of E* the magnetopause appeared substan- 
tially more reflective to electrons than to ions. 
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Figure 5. Magnetic field reconnection without plasma energization ^ 
Magnetopause is the wedge-shaped region between the two rotational dis- 
continuities. Quantities in this region are denoted by subscript 3. 
Note that |b |=|b |=|b I and jv l=|v i=|v I. 

3 1 2 3 12 

(d) Reconnection geometries do exist in which the reconnection 
component of the electric field is directed along the net current 
but where no particle energization occurs. Figure 5 shows an example 
of Such a geometry. The magnetopause is wedge-shaped and contains a 
strong magnetic field perpendicular to the plane of the figure. At the 
t^o edges of the magnetopause, rotational discontinuities deflect the 
flows incident from both sides without change of the particle speed. 

In the inflow regions the flow is predominantly perpendicular to the 
plane of the figure while in the outflow (i.e., in the magnetopause) 
the flow vector lies in that plane. Prom the point of view of the 
Poynting-vector flux, the inflow of electromagnetic energy from both 
sides is exactly balanced by an outflow in the wedge shaped magneto- 
pause regions so that no conversion of electromagnetic to particle en- 
ergy occurs. Mathematical details of such models may be extracted from 
the work of Cowley (1974) . 

! There is no observational information that would exclude the pos- 
sibility of a wedge shaped magnetopause. And the magnetopause often 
contains a transverse field of magnitude comparable to B. and B. . 

"'ll j ^ 1 2 

But the model does not predict a plasma boundary layer inside the mag-* 
netopause. Rather, the magnetopause and the boundary layer cpincide, 
and plasma is fed into this layer from both sides • 

(e) In the light of existing observations of the magnetopause 
structure it seems likely that reconnection, if it occurs at the 
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i^gnetopause* Is a patchy and highly time-dependent process. Assume 
that, on the average, reconnection signatures may be seen over an area 
equivalent to a strip of width A = 1R„ on the front lobe of the magne- 
tosphere. With a normal magnetic field component B = 5-lOnt such a 
dtrip can reconnect a sufficient amount of flux to replenish th^ tail 
flux in a few hours. The probability P of a satellite encountering a 
reconnection patch is equal to the area of the strip divided by; the 
drea of the entire front lobe. With a magnetospheric cross sectional 
ijadius of R = ISR^ one then finds P = TrR^A/2fTR^ = 1/30. I'/hen ins'trU- 
njental limitations of the plasma experiments are folded into the pic- 
tjure it seems possible that reconnection could occur even though a 
cjlear signature of it has not yet been seen^* in t.he plasma data. 

I It should be remembered that various magnetic signatures compati- 
ble with reconnection, such as the presence of nonvanishing value of 
and rotational behavior of the tangential magnetic field in the mag- 
netopause, have been seen (e.g., Sonnerup, 1976; Somierup and Ledley, 
1974,1979). Sonnerup and Ledley (1979) suggest that patchy reconnec- 
tion should produce indentations in the magnetopause. Such indenta- 
tions, when swept along the magnetopause, lead to rapid multiple boun- 
dary crossings. Thus it may be that the search for the plasma signa- 
tures of reconnection should focus on, or at least include, multiple 
crossings. Such crossings were excluded in a recent study of IMP6 data 
by Eastman and Hones (1979^ • 

■ i . ■ ■■ ■ ■ ;•••■ ■■■• ■■ ■ 

I In summary, it seems that a combination of some of the effects 

(a)-“(e) may be able to account for the discrepancy between classical 
reconnection theory and presently available observations. Better data, 
supplied by the ISEE and other experiments, and more detailed inter- 
comparison of magnetic field, electric field, and plasma data are 
needed before one can assess the impoi*tance of plasma convection across 
the magnetopause (with or without reconnection) in a reliable way. 


i Recent laboratory experiments by Podgorny et al., (1978) and 
Dubinin et al., (1978) seem to establish the possibility of two topo- 
logically distinct niagnetospheres, depending on the ratio of flow speed 
to! Alfven speed in the incoming flow. For low values of that ratio and 
a southward interplanetary field, a Dungey type configuration (Dungey, 
19J61) is established with magnetopause reconnection in the equatorial 
plane. For high values an entirely different geometry is observed with 
a ’Visor**, containing closed field line loops, covering the front lobe 
ofj the magnetosphere. It is suggested here that in the real magneto- 
sphere this visor has a thickness comparable to the ion gyroradius (as 
it, did in the experiment) and that it represents the dayside magneto- 
pause of a closed magnetosphere. If this interpretation is reasonable, 
th^n the magnetosphere with a southward interplanetary field would have 
two states: fully open for low values of the ratio of flow speed to 

Alfven speed, closed for high values of that ratio, the latter being 
the normal situation. A suggestion to this effect has been made by 
Sonnerup (1965) . It is also in accord with the work of Cassen and 
Szabo (1970) , who observed that viscous boundary layer solutions in 


field-aligned flow cease to exist below a certain critical ratio of 
flow speed to Alfven speed. They suggested a "disrupted" magnetopause 
in such circumstances. 

The patchy reconnection discussed above would appear as a pertur- 
baltion on the closed state. The question then remains whether the 
fully open configuration ever occurs. It would seem possible that ic 
does during magnetic storms that involve large southward interplanetary 
fields. Sonnerup (1971) has reported one possible direct observation 
of the reconnection line in the equatorial plane during a magnetic 
storm. 


6. ENERGY DISSIPATION ‘ 

It is of interest to consider briefly the expression for the elec- 
tromagnetic energy dissipation rate per unit area during reconnection. 
Using the notation in Figure 4b we find 

. 2 = EtI cosa (28) 

where a is the angle between and _I and the magnetopause current per 
unit length is 

I = (B. 7y )A+(B^ /B^ )^-2(B^ /B^ )cose 

tl 0 t2 ti t2 ti 

(2B. /y )sin(e/2) , (29) 

tl 0 

Here 0 is the angle between B and- B^ . For B^ = B^ , i.e., when the 

— tl — 12 t2 tl 

field magnitudes are equal, the last member of the equality results. 
Below, the last member in each formula for (P ref ers to this case. At 
present, the electric field magnitude and angle a cannot be pre- 
dicted from theory. However, it is of interest to explore the conse- 
quences of two reasonable sets of assumptions concerning these quanti- 
ties . : 

Assume that we take _E to be proportional to the unperturbed in- 
terpl anetary electric field vector so that E^ = K v^B^ (where 

+ b| in solar magnetospheric coordinates ^andv^ is the solar 
wind speed) and a = 90 - (The angle ^ is defined in Figure 4b.) 

In that case 

(P ■ 

(2K V B^/y )(B^ /B )sin^0/2 (TT/2<e<lT) (30) 

I I OO CX3 Q I; J OO 

The coefficient Kj may be a function of 0 as well. It is not clear 

precisely how B /B or 6 should be related to conditions upstream in 
t 1 

the solar wind. But it is evident that the power input will be a strong 

function of 0/ 2 with a maximum when B^ is antiparallel to B. (0=7i; 

tl. . : ‘ ^ . t2 
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ihe interplanetary field is due south). This model also gives (p- 0 
for O<0<Tr/2. In other words, it gives a "rectifier" type of behavior 
(Burton et al., 1975). 

Another possibUity is to insist that the electric field JE^ is 
along the current ^ and thct its magnitude is determined by a Pctschek 
type upper limit oh reconnection (see Sonnerup, 1974); 

= K V, B sin? (K =0.1-0. 2) (31) 

t 2 Aj_ t l 2 

where v. = B. sinaA^p p 

Ajl tl 0 1 

We then find ' 

K Bf /B/ cos0)^ 

^ 2 tl t?. tl 

^ 4Tp“ /1+(B^ /B^ )''-2(B^ /B^ )cos0 ' 

0 '^0 1 ta tl t2 L) 

[2K B^ /(p ,^m]sitvHe/2) * (32) 

2 tl 0 0 1 

In this case too, a very strong dependence on the angle 0/2 is 
predicted. 

From empirical considerations Perreault and Akasofu (1978) propose 
that the total energy input into the magnetosphere is^ 

*' (^/2) (33) 

OO CO ' ^ 

where 0 = ir/2 - tan * (B^/B^) . 

It is unlikely- that any I of the preceding expressions, which dcs- 
c^ribe only the magnetopause Dissipation, can be cast into this form. 
I|ut, assuming 0 and 0 to be proportional, the reconnection process 
deems to agree with the observational result that the power input de- 
piends strongly on the angle 0 with maximum input at 0 = IT and zero in 
put at 0 = 0 (or possibly for 0 less than some minimum value). 

7. DIFFUSION 

The reconneGtion model is not the only one to yield a strong de- 
pendence of the energy input on the angle 0 between the fields on the 
two sides of the magnetopause. In a purely resistive magnetopause the 
electromagnetic power dissipation rate per unit area is 

I •= n[i]i^ (34) 

where r|[l3 is the resistivity, assumed to be a function of the current 
I. With the expression (29) for the current, and for B = B , we 

find “ 

(P-^ (2B. /V )^sin^(0/2)n[(2B /)J )sin(0/2)] (35) 

■ ■ ■ ■ ■ ■ , t'l," 0. ■ ■ ■ , ' 0" ■■ 


~ 20 - 


If ri Increases rapidly with increasing current as it may when the re- 
sistive effects are produced by current or gradient driven microinsta- 
bilities, we find the electromagnetic dissipation rate to be a very 
strong function of 0 with a minimum at 6 = 0 and a maximum at 6 = tt. 

In a diffusion dominated situation, most of the energy input into 
thie magnetosphere is in. the form of mechanical and thermal energy of 
the particles which diffuse across the magnetopause. The electromag- 
nejtic portion given above is small. However, presumably this latter 
energy goes into the establishraent. of the microturbulence which in turn 
allows particles to diffuse across the magnetopause. Thus (P becomes 
the controlling factor for the diffusion so that its 6 dependence gov- 
erns the total mass and energy flow across the magnetopause. 

The detailed nature of the microinstabilities operative in the 
magnetopause is not known. Eviatar and Wolf (1968) have suggested 
two-stream ion cyclotron instability , Huba et al,, (1977) have proposed 
the lower-hybrid drift instability, and Hasegawa and Mima (1978) have 
calculated diffusion coefficients in the presence of kinetic Alfven 
wave turbulence. On the basis of these results it appears reasonable 
to assume that effective kinematic viscosities V in the required range 
of 10^-5 X 10^ m^/s (Axford, 1964) can be generated via microinstabili- 
ties in the magnetopause. 

The hydrodynamic Reynolds number based on a flow speed v - 100 
km/s, a magnetopause thickness h = 200 km, and a kinematic viscosity 
y = lO^.m^/s, is Re = vh/v == 20. A sheet jet in ordinary hydrodynamics 
becomes turbulent at Re = 30. I'Jhen account is taken of the stabilizing 
effects of the sheared magnetic field at the magnetopause, and of the 
fact that the flow is strongly accelerated, it appears unlikely that 
hydrodynamic turbulence and the associated eddy diffusivity is impor- 
tant near the subsblar point. (In the tail magnetopause, it may well 
be dominant.) 

Dayside plasma boundary layer observations reveal a bewildering 
array of possible structures (Eastman and Hones, 1979; Paschmann et 
al., 1978) including situations where little or no boundary layer is 
seen. The statement that such a layer is present implies that magneto- 
sheath-like plasma is seen inside the inner edge of the magnetopause. 
Often this plasma layer is much thicker than the current layer. If 
such a situation is to be explained by diffusive processes then the 
diffusion coefficient for mass must be much greater than that for cur- 
rent. Such differences are not uncommon for molecular diffusion pheno- 
mena (e.g., momentum diffusivity greatly exceeds heat diffusivity in 
laminar flow of oil); they may occur in a plasma as a result of micro- 
instabilities as well. Hasegawa and Mima (1978) find that kinetic 
Alfven waves yield a diffusion coefficient for electrons that is much 
greater than that for the current. Papadopolous (private communication, 
1978) has suggested that the unequal thicknesses of the current and the 
plasma layer may be accounted for by a turbulent wave field generated 
by the intense currents in the magnetopause and then spreading to the 
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adjacent boundary-layer region allowing effective mass diffusion there. 

I - 

i ' - 

On the other hand, in turbulent hydrodynamic boundary layers the 
dif fusivities for mass, momentum, and energy are approximately equal. 
Similarly, in the magnetopause application, if eddy diffusivity asso- 
ciated with gross hydrodynamic turbulence provided the dominant trans- 
port mechanism, then the magnetopause current layer would be as thick 
as the boundary layer and the two layers would occupy the same region 
in space. 

On the basis of the discussion given here, we conclude that dif^ 
fusion produced by plasma micro turbulence may not be ruled out as an 
Important, occasionally even dominant, transport process at the day- 
side magnetopause. Large-scale hydrodynamic turbulence appears to be 
a less likely agent, at least near the subsolar point. It must be 
pointed out, however, that the boundary layer is sometimes observed to 
^ave features which are difficult to explain on the basis of diffusion. 
For example, extremely steep density gradients have been seen at the 
magnetopause, followed by a nearly constant density in the boundary 
layer with another sharp drop-off in density at the inner edge of that 
layer (Paschmann et al. , 1979). 
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NOTES 

1. On leave from Dartmouth College, Hanover, N.H. 03755, U.S.A. 

2. The polar angle 0 is not to be confused with the angle 6 between 
the magnetospheric and the magnetosheath field. 

Another possibility, leading to similar results, is that particles 
of inagnetospheric origin flow outwards across the magnetopause. 

4.i After completion of this paper, a few ISEE magnetopause crossings 
have been found with substantially enhanced plasma flow speeds near the 
taagnetopause (Paschmann, private communication, 1979) . 

5. Rossberg (this conference) has shown improved correlation when 
B^ = B| + b| rather than = B^ + B^ + 
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ABSTRACT 


After a review of the properties and predictions of 
the closed and open models of the magnetopause, 

C^GO-5 magnetometer data are used to illustrate var- 
ious observed signatures of the magnetopause current 
layer and the adjacent plasma boundary layer • Among 
the topics touched upon are: fluctuations, diamag- 

netic effects, and field-aligned currents in the 
boundary layer; one- dimensionality of the magneto- 
pause; presence and absence of a magnetic field 
component perpendicular to the magnetopause; finite 
ion gyroradius effects. A brief summary is given 
of existing Vlasov theory for the description of 
tangential, rotational and contact discontinuities. 
Special attention ie paid to the tangential momentum 
balance and the jump conditions at a rotational dis- 
continuity. Finally, a discussion is given of low 
frequency fluctuations with emphasis on the' signa- 
tpres of the tearing mode. 


Keywords: Magnetopause, Boundary Layer, Tangential 

Discontinuity^ Rotational Discontinuity, Tearing 
Mode 


1 . INTRODUCTION 

The magnetopause is the thi.i layer in space which 
marks the outer* boundary of the earth Vs magnetic 
field. Immediately outside the magnetopause one 
finds the streaming solar-wind plasma of the mag- 
nejtosheath and its imbedded interplanetary magnetic 
field of time-variable direction, draped around the 
magnetosphere. Immediately inside , the magnetic- 
field direction is less variable and is in large 
measure controlled by the earth. The magnetopause 
Itself is an electric current layer which adjusts 
the magnetic field direction and magnitude from the 
magnetosheath to the magnetospheric state, in ac- 
cordance with Ampere’s law. Inside the magnetopause 
one finds a highly time -variable layer of streaming 
plasma of magnetosheath origin, at the inner edge 
of j which the magnetic field may undergo further, “ 
usually much smaller, changes in magnitude and 
direction. This layer is referred to as the bound- 
ary layer. 

The present paper attempts to briefly summarize what 
Is [known, and not known j both observationally and 
theoretically, about the electromagnetic structure 
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of the magnetopause and boundary layer regions . 
'Over the years, a substantial body of observationa 
facts has accumulated concerning the magnetic 
structure (Ref. 1-14) . Direct measurements of the 
DC electric fields near the magnetopause are far 
more difficult and results have started to emerge 
only very rece;ntly (Ref. 15; also A. Pedersen, thi 
conference) . Fluctuating electric ani magnetic 
fields have been studied only to a limited extent 
(Ref. 9, 16^17) , The existence of the plasma 
boundary layer has been established relatively 
recently (Ref. 18-26). The properties of this 
layer have been observed in fair detail (see revie’ 
by G. Paschmann, this volume) but its generation 
and relationship to the adjoining magnetopause 
remain obscure. From a theoretical viewpoint, the 
electromag net. Ic and plasma structures of the mag- 
netopause and boundary layer are of course inex- 
tricably intertwined, Thus, the theoretical dis- 
cussion to be given will, by necessity, deal with 
both aspects. 


There are three major scientific reasons for study 
of jthe magnetopause and boundary layer: (1) From 

the point of view of solar-planetary physics, it i; 
essential to understand the transfer of mass, mo- 
mentum, and energy from the solar wind, across the 
magnetopause, into the magnetosphere. Included 
among possible processes are diffusion due to 
micro- and macroturbulence as well as convective 
entry in association with magnetic-field reconnec- 
tion. (ii) Current sheets appear to play an 
important role in many cosmic systems. Any physi- 
cal understanding of the structure of, and dynami- 
cal processes in, the thin collision-free magneto- 
pause current layer may be directly used in other 
costnic contexts. In particular, if we can learn 
to understand the reasons for the occurrence, or 
absence of occurrence, of recohnection at the mag- 
netopause, and in the former case, if we can estab- 
lish the efficiency and physical ^nifes tatxons of 
that process, then our level of understanding of 
the conversion of magnetically stored energy in 
cosmos into plasma kinetic and thermal energy will 
be significantly advanced. (iii) The magnetopausi 
provides an excellent opportunity for plasma phys- 
icists to develop and test theories for linear and 
nonlinear micro- and raacroprocesses in current- 
carrying collision-free plasmas. 


The magnetometer data presented in the paper were 
obtained with the Goddard Space Flight Center 
fluxgate magnetometer ( J . Heppner, principal in-’ 
vestigator) onboard the satellite OGO-5. 


2 . GLOBAL CONS I DERATIONS 


For given interplanetary and iiiagnctospheric condi- 
tjlons the time- average magnetopause and boundary- 
layer structures presumably vary in an organised 
in|anncr with location. For example, it appears that 
the boundary layer thickness increases systematic- 
ally with increasing distance from the subsolar 
point (Ref, 21). However, because of unknown, 
usually rapid, radial motion of the local magneto- 
pause, the absolute thicknesses of the current; 

Ijayer and boundary layer have been difficult to 
aissess precisely (for a summary, see Ref. 27), at 
Ibast before the IS EE mission. Furthermorei, struc- 
tural changes associated with changing interplane- 
tary and magaetospheric conditions often obscure 
effects associated with different locations on the 
nagnetospheric surface. Since no systematic obser- 
vational studies exist which separate all of these 
effects^ we resort to a simpler approach: only 

structures on the front lobe of the magnetosphere 
will be considered, during conditions where there 
is a large and relatively abrupt change in field 
direction at the magnetopause, making it easily 
Identifiable in the magnetic data. It must be 
remembered that there are many recorded instances 

where these conditions are not met. 

( 

The closed and open models of the magnetosphere 
have played, and continue to play, an important 
role in the interpretation of observed magneto- 
pause structures. The former has its origin in 
the work of Chapman and Ferraro (Ref. 28). Note 
that in its modern form, this model does not imply 
ajtall of finite length but merely a vanishingly 
si^ll magnetic field component normal.! to the 
magnetopause. The latter model j in a simple form, 
was introduced by Dungey (Ref. 29)^ For our pur- 
poses, its most important property Is that ^ 0. 

Recently, these two conceptual models have been put 
on a firmer physical basis by the laboratory simu- 
lation experiments of Podgorny and coworkers (Ref. 
30). They have shown that for antiparallel fields 
on the two sides of the magne topause . the open con- 
figuration occurs for flow speeds comparahle to, or 
less than, the Alfven speed while a configuration 
similar to the closed model (but with a dispropor- 
tionately thick magnetopause, presumably resulting 
from the large ion gyroradius in the experiment) 
occurs at high flow speeds. In the real case, the 
solar-wind flow speed is large compared to the 
Alfvdn speed so that one might normally expect the 
closed model, perturbed by those micro- and macro- 
processes, including small-scale ■■patchy" recon- 
nection, which may operate in the real magnetopause 
but perhaps not in the simulated one. Because of 
Imperfect scaling, the conditions for which a 
transition between the two models occurs may be 
substantially different in the laboratory and in 
the real case. Thus, one cannot exclude the pos- 
sibility that the Dungey model occurs in nature. 

In the closed model the plasma boundary layer inside 
the magnetopause must be formed either by local 
diffusion of magnetosheath plasma across the mag- 
netopause or by nonlocal entry, e.g., at the cusps. 
In the open model, the plasma flows directly across 
a, possibly rather limited, longitude sector of 
the dayside magnetopause, as a result of the pres- 
ence of a nonvanishing normal magnetic field conipo- 
nent* In crossing the magnetopause the plasma is 
accelerated by the i X ^ force into two polevs^ard- 
directed jets, one in each hemisphere. Here ^ is 
the magnetopause, or Chapman-Ferraro, current. 


TltcsG jets supply the boundary layers of the open 
model . . I 

1 

Iii the closed model the plasma flow Is tangential " ^ 
and the associated electric field Is normal to the mag- 
ne,topause. Additionally, the magnetopause struc- : \ 
tiire itself is likely to contain a normal electric: ^ 
field. In the open imodel, normal electric fields 
are also present. But in addition, there is a 
tangential component at the magnetopause. This 
field leads to the drift of magneto sheath plasma 
into the magnetopause. It is associated with the 
occurrence of magnetic-field reconnection in or 
near the equatorial plane (reconnection in the 
cusps also remains a viable possibility), and its 
direction is such that * i > 0* 
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Figure 1. OGO-5 magnetometer record of crossing 

of the magnetopause (MP) and boundary , j 
layer (BL) , showing the field magnitude 
(ly = Int, dashed line), GSM latitude ^ , 
(6, solid line) and longitude (V^ dottec i 
line) . Satellite location (GSM): ■ J 

X « 48130, Y = - 202040, Z = 3330 km. 

I 3 . OBSERVED TIME SIGNATURES H 

iFigure 1 shows an OGO-5 magnetometer record of a 
traversal from the magnetosheath ^ across the mag- v . 
netopause and boundary layer, and into the mag- 
netosphere. The magnetic field is represented in ; ’ 
terms of its magnitude, GSM latitude 9, and longi- 
tude The latitude angle is initially positive ♦ ; 
(4*25°), indicating that the magnetosheath field ^ j 
points somewhat north of the equatorial plane. ^ ^ 
The angle turns negative as the magnetopause is ■ 
entered, and then positive again until it reaches - i 
the raagnetospheric value (+85°) at the inner edge : j 
of the magnetopause. The field magnitude has a 
broad minimum during the 80 sec period comprising 
the magnetopause but it then remains low for ' 

another 60 sec. This latter period corresponds 
to the traversal of the boundary layer. The dia- 
magnetic effect associated with the boundary-layer 
plasma is evident from the abrupt increase in field i 
magnitude at the inner edge of the layer . It cor- j | 
responds to a density change of 112 protons/cm^ ^ 
at an assumed temperature T - 2 x 10^°K. This 
unusually high density and the observed high field i 
magnitude are associated with a magnetopause loca- \ 
tion at 8.2 Rg . Note also that the field magnitude 
has a substantial maximum just inside the inner 
edge of the boundary layer, Neugebauer e t al . (Ref , ^ 

9 ) have ^^ssociated this latter effecc with the loss j 
of magnetospheric trapped particles within a gyro- ^ 
diameter of the magnetopause. 

■ - if 

Large fluctuations of both field magnitude and I 

direction are present, not only in the magnetopause. 



but In the boundary layer as well. Tills is a 
characteristic feature of these re^»ions (Ref. 21). 
Inside the boundary layer the field magnitude and 
latitude angle exhibit only slow variations. The 
longitude angle continues to fluctuate but since 
0 Is near 90® the corresponding direction changes 
are relatively small. 

It is noted thaty^changes in a systematic way at 
the inner edge of the boundary layer. This change 
may correspond to a field-aligned current sheet 
flowing from the boundary layer into the ionosphere 
In the manner described by ETastman et al. (Ref. 21). 
An alternative explanation, which is less likely on 
account of the high fluctuation level in the fv.id, 
is the finite ion-gyroradius effect described oy 
Parker (Ref. 31, 32) and discussed further in Sec- 
tion 5. 

If the magnetopause were stationary, the time dura- 
tion of the magnetopause would correspond to a • 
thickness of about 200 km. However, thickness 
estimates based on the analysis of multiple cross- 
ings and, more recently, on ISEE data (Ref. 14) 
Indicate typical magnetopause thicknesses in the 
.range of 500-1000 km, i.e., several proton gyro- 
diameters. In the present example it appears that 
the boundary layer was somewhat thinner than the 
magnetopause. In fact, unless one takes care to 
define the latter as the region where the main 
field direction change occurs, one may be tempted 
to incorporate the boundary layer into the magneto- 
pause and thus fail to identify it as a separate 
region. 
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Figure 2. 0 GO- 5 magnetometer record of boundary 

layer (BL) and magnetopause (>IP) cros- 
sing. (For further information concern- 
ing the March 27, 1968, erosion event, 
ste Ref. 7). (X = 54050, Y = - 45540, 

Z rs 3590 km.) 

Figure 2 shows another OGO -5 record, this time of 
an outbound crossing, which displays a clear mag- 
netic signature of the boundary layer, namely a 
depressed field magnitude and large fluctuations. 

It should be noted, however, that in many cases the 
field depression in the boundary layer is weak or 
absent I presumably because the pressure in the hot 
tenuous ring-current plasma is about the same as 
the pressure in the cooler but denser boundary layer 
plasma. Tire deep field magnitude minimum in the 
magnetopause itself is a very common feature. 

Til e field magnitude levels in Figure 2 are more 
representative of typical conditions at the magnet- 
opause, Reference 11 gives average values of 28 nt 
outside the magnetopause and 42 nt inside the bound- 


ary layer, corresponding to a density change of 
14 protons/ cm^ at 2 x 10 ^ “K. 

Theoretical aspects of the electromagnetic struc- 
ture of the boundary layer are discussed in a con- 
tributed paper (Ref. 33). In the remainder of the 
present report we shall focus on the magnetopause 
Itself, 


4. OBSERVED MAGNETOPAUSE STRUCTURES 


It is of interest to ask whether the magnetopause 
structure is approximately one dimensional or 
intrinsically two or three diraensionai. The former 
case implies that variations along the layer occur 
only on length scales much greater than its thick- 
ness. It then follows from ^ ^ = 0 that the 

magnetic field component normal to the magneto- 
pause, remains very nearly constant across the 
layer at any location and instant. The latter case 
implies that variations along the sheet occur on 
length scales comparable to the thickness. In such 
circumstances usually does not remain constant. 
If the magnetic vector data set taken during a 
magnetopause crossing has the property that a 
unique directibn ’ can be found along which all the 
vectors have one and the same component, then 

a strong case can ibe made to the effect that this 
direction is perpendicular to a one-dimensional 
sheet. Alternate possibilities are (i) that the 
two tangential derivatives in £ • _3 = 0 are large 
but equal and opposite; (ii) that time variations 
or jvariations in the attitude of the magnetopause 
exactly compensate .for existing space variations 
In Bn. These alternatives require very special, 
and therefore highly unlikely, circumstances. A 
fourth and more likely alternative will be discus- 
sed at a later point. On the other hand, if one 
cannot find a direction haying constant or nearly 
constant Bn from the data set, ^ this does not 
necessarily imply that the local magnetopause 
structure was two or three dimensional. Time and 
attiitude* variations may have caused the noncon- 
stancy. 


In j practice one usually cannot find a direction 
In jspace along with B^ is_ strictly constant,. In- 
stead one uses, as an approximation, that direc- 
tion which yields minimum variance in the corres- 
ponding field component (Ref. 2). This d-'rection 
is along the eigenvector N 3 corresponding to the 
smallest eigen value (= variance X 3 ) of the matrix 


Mi^ = - BiBj 


( 1 ) 


which can be formed from the magnetopause magnetic 
vector data set by appropriate averaging (denoted 
by overhead bars). The vector _N 3 is taken to be 
directed away from the earth. The two other eigen- 
vectors and 2 ^; corresponding to the largest 
and the intermediate eigenvalues, and X 2 , are 
tangential to the local magnetopause surface and 
due approximately north and west, respectively. 

The right-handed orthogonal set (Nj, _N 2 » pro- 
vides a convenient natural coordinate system in 
which to examine the magnetic field data. Error 
estimate formulas for the normal vector and 
the normal component F 3 may be found in Ref, 8 
and 11. These estimates represent lower limits 
because they do not include any systematic effects 
such as attitude or time changes during the cros- 

s^^^C- 


It is sometimes found that the matrix Mj^j is 
nearly degenerate with A 3 and A 2 almost equal. In 





that case an inaccurate determination of results 
and the error estimates, containing terms of the 
form (X 2 ~ ^ 3 )'"^, «aro correspondingly large. It is 
easy to show that such a state of affairs arises 
when the magnetopause current-density vectors are 
uhidirectional or nearly i nidirectional , As an 
illustration, consider the case of a two-dimensional 
tearing structure in the magnetopause. The current 
is unidirectional but not uniform, being concentrat- 
ed to the centers of the tearing islands. The min- 
imum variance analysis applied to such a structure 
gives X 3 = 0, X 2 small, and Xj^ large. The vector 
is along the current, i.e. , it is tangential, 
rather than normal^ to the magnetopause surface! 

When the irlnimum variance analysis is applied to 
an optimal regment of the magnetopause crossing in 
Figure 1, the eigen values (X^, X 2 ^ X’ 3 ) are found 
to be (1890, 338, 45) showing that i the varianee in 
the field component along ^3 was large but still 
much smaller than the variances corresponding to 
N? and 1^. Thus, the magnetopause can be said to 
have been one^dimensional in an average sense, but 
with large superimposed two or three dimensional 
fluctuations. 

Tile estimated normal vector points in a reason- 
able direction, given the location of the space- 
craft, and has an estimated error of ± 6 ®, The 
normal magnetic field component. Bp = B ’3 = 7.6 ±6.7®, 
is not significantly different from zero. The 
crossing is shown in polar, or hodogram, form in 
Figure 3. The plot on the left represents the 
behavior of the tangential field (B^^, B?) , the one 
on the right shows the normal field B 3 . A constant 
value corresponds to a vertical trace' in the 
latter plot. The high level of turbulence peaking 
in the 0,01-0.10 Hz range is perhaps the most 
striking feature of the diagram. 

Figure 4 shows a hodogram pair for the crossing in 
Figure 2. The eigen values are (480, 22, 5). The 
estimated error in N 3 is only ± 6 ® but it is seen 
that, apart from several minor excursions, the 
tangential hodogram indicates a fairly constant 
value of B 2 , corresponding to a nearly unidirection- 
al current along ^^ 2 . Thus, the normal direction 
is strongly influenced by the properties of the 
aforementioned minor field excursions and should 
not be trusted. Use of different data segments 
for the analysis also tends to give inconsistent 
results. This example illustrates the importance 
of examining the hodogram pairs before one draws 
conclusions about the reliabilifrv of the calculated 
normal vector and normal field uOmponent , 

Figure 5 shows an example of a high quality deter- 
mination of the normal vector, with an estimated 
angular error of only ±2.4®, and of the normal field 
component B^^ = ^3 == 0.1 ± 0.4 nt. The set of eigen 
values for this case is (747, 82, 2) indicating 
that the layer is one-dimensional to a good approx- 
imation. Note that the fields on the two sides of 
the magnetopause are very nearly anf:iparallel , yet 
the magnetic trademark of reconnection, a B^^ value 
substantially different from zero, is entirely 
absent. It seems clear that factors otlier than the 
orientation of the magnetosheath field relative to 
the earths field play an important role in control- 
ling reconnection. 

The nature of the tangential field in this outbound 
magnetopause crossing is remarkable. The field 
reverses direction by a rotation rather than by a 
decrease of the northward field to zero followed 
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Figure 3. Hodogram representation, with field . 

components in nt, of a portion of the ■ ? 
magnetopause crossing in Figure 1. The f | 
magnetosphere has > 0. GSM projec- ^ ' 
tions of (N-j , N»> , 2 ^) on the lower left. 

(X -- 48130, Y = - 20240, 2 - 3330 km.) ; ] 
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Figure 4. Hodogram representation of magnetopause . I 
crossing in Figure 2. (X = 54050, 

Y = - 45540, 2 =3590 km.) r- 

by an increasing southward field as the magneto- - ■ 
sheiath side of the magnetopause is approached. 

Such behavior indicates the presence of substantial | 
field-aligned, i.e,, highly muItidireGtional , cur- ^ I 
retits in the magnetopause. It is this field rota- ^ ! 
tibn, and the corresponding large variance X 2 that 
helps produce an accurate normal vector de termina - y 1 
tion. ■■ ■ . ^ 


Thd crossing in Figure 5 was extremely rapid, oc- 
curring in a time span of only 9 sec. Thus one 
may think that a snapshot picture of the magneto- 
pause structure was obtained. However, the cros- 
sing Was followed within less than 10 sec by a 
second even more rapid one back into the magneto- 
sphere (see Ref. 11, Figs 4 2 and 3). The structure 
of the second crossing was sufficiently different 
from Figure 5 to suggest that structural changes, 
ipresumably associated with rapid convection of 
magnetic stVuctures along the magnetopause, occur 
on very short time scales. In other words, the 
hodogram obtained during a magnetopause crossing 
Is often severely time aliased. It should be 
added that the normal vector orientations for the 
pair of crossings discussed above suppor.t the 
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Figure 6 * Hodogram representation of magnetopause 
crossing with F 3 0, (X = 33550, 

Y «= - 10950, Z = 17730 km.) 


Figure 5. Hodogram representation of magnetopause 
dressing for which an accurate deterini^ 
rtation^of N 3 and ¥3 is obtained, and for 
which B 3 - 0. (X = 56480, Y = - 30080, 

Z « 6590 km.) 

Figure 6 shows another example of a good normal 
vector determination and a vanishingly small normal 
field component == 3 * 2 = 1.2 ± 1.3 nt. Again, the 
fields on the two sides are nearly antiparallel, 
but this time the structure of the tangential hodo- 
gram is rather more complicated. 

As Reported in Ref. 11, values of |E* 2 | less than 
the error estimate are obtained in approximately 
25% of the crossings. Values greater than three; 
times the error estimate are also found about 25% 
of the time. T\sro examples of the latter situation 
are shown in Figures 7 and 8 . In the former case 
F 3 = - 5.1 ±0.9 nt, in the latter ~ + 3. 1 ± 0 .4 nt 
Results of the minimum variance analysis for five 
J data segments of the ^econd crossing (Ref. 7, 13) 
j are sho^^m in Table 1. The consistency of the re- 
1 suits, except for the shortest segment, the flatness 
' of the B 1 B 3 hodogram, and the extremely large sep- 
aration between X 3 and X 2 ^ the optimal segment, 
provide convincing evidence Chat the magnetopause 
does Indeed on occasion develop a substantial normal 
magnetic field component. This crossing occurred 
In a time spa.n of only 4 seconds so that an approx- 
imate snapshot of the structure may have been 
obtained. 


5. TANGENTIAL DISCONTINUITIES 

Current sheets having 0 are referred to as 

taru'ential discontinuities (Ref. 34) . Self con- 
sistent bne-dimensional Vlasov equilibria of such 
sheets can. be generated by allowing the distribu- 
tion function f to depend on the constants of motion 
of a particle: the total energy and the two 

generaj-ized linear momenta tangential to the sheet. 


and ? 2 . 

f «= f(f, Pi, P 2 ) i ( 2 ) 

£ = >5 mv^ + e<!> (3) 

*= + eAjL (4) 

?2 = tnv2 + eA 2 (5) 


Here v = ^3^ velocity of a particle 

of mass m and charge e. Also, ‘K^b) ^nd A= (Ai(x 3 ) , 
A 2 (x 3 ) , 0 ) are the scalar and vector potentials , 
respectively, X 3 being the coordinate along the 
direction normal to the sheet. 

During ^he last 20 years large numbers of equilibria 
of this type have been generated (for reviews see 
Ref. 35-37). The most recent work by Lee and Kan 
(Ref. 38) and by Roth (this confereaca) shows prom- 
ise of being able to account for several observed 
magnetopause features. In general, by a sufficient- 
ly clever choice of the distribution functions, one 


Table 1. 
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Number of 
B Vectors 


Eigen Values 

Y2 

Normal Field 

Error 

in N 3 

Segment 


^2 

^3 

B'3 ± 683 , Y 

ANs'Ni 

^3* N 2 

1 

48 

76 

1.5 

0.15 

3,4 ± 3.4 

0.02 

0.13 

2 

96 

267 

23 

0.58 

10.0 ± 1.2 

0.02 

0.06 

3 

144 

463 

83 

1.3 

8.1 ± 0,7 

0.02 

0.04 

4 

192 

570 

124 

1.6 

8,0 ± 0.5 

0.02 

0.03 

5 

240 

722 

130 

1.9 

8 . 1 -±0.4 

0.01 

0.03 

The five 
Flf ty-slx 

nested data segments are centered at 
vector samples per second . 

65938239 ms UT. 








lit rates df 50-100 km/s^*. Such accelerations^ even 
If opcraL|lvo only durlnj* the time It takes a space- 
craft Co cross the map,ne.topausQ , which would be an 
unlikely coincidence* leads to unacceptably high 
plasma' velocities. The remaining possibilities are 
1:0 balance by changing the tangential momentum 
of magnetos heath plasma convecting across the cur- 
jfent sheet or by utilising the tangential component 
pf the plasma pressure tensor which develops when 
tlhe pressures p parallel and perpendicular to 
ihe magnetic field are unequal. This latter stress, 
incorporates effects of particles reflected at the 
current sheet, live convective stress is jr^ ® -pVnAvj. , 
vhern p and are the plasma density and flow com- 
ponent along N^,; respectively, while is the 
change in tangential flow velocity as the plasma 
crosses the sheet. The net tangential stress as^ 
sociated with pu is a = 

(P|j'-Pj^)(Pq/b’) is the pressure anisotropy factor. 

The net tangential momentum balance now becomes 
(Ref. 40): 

pVjjAvj. = (6) 

Tills equation applies to fast and slow shocks, for 
which ^ 0 and At t»elng the angle formed 

between the tangential magnetic field components on 
the two sides of the sheet. It also applies to 
rbtational discontinuities (v^ 0, 0<A^<7 t) and to 

contact discontinuities (v^^ = 0, At 0, it) . The 
rotational discontinuity is of key importance for 
the magnetopause application, because it is the only 
type of current layer with r 0 that permits of a 
tangential field rotation by an arbitrary angle Aj.. 
This important point was first recognized by Levy 
et al, (Ref, 41) . 

For the rotational discontinuity it may be shown 
(Ref. 40) that 

Pl/P2 

where the subscripts 1 and 2 refer to conditions on 
the magnetoshca th and magrietosphere sides of the 
magnetopause, respeGtively. Purther, the normal 
flow speed is the Alfvdn speed, modified by the 
pressure anisotropy factor: 

|vnl = ( (8) 

When plasma as well as magnetic field data are 
available for a magnetopause crossing, these formu- 
las may be used to check the hypothesis that mag- 
netosheath plasma flows across the magnatopause in 
the presence of 0. First, without knowledge' 

of Vn and one may check whether, the two vectors 
AZt are parallel (or antiparallel) 

asj required by Equation (6). If that is the case, 
one may then obtain the ratio Vn/B^ from Equation (6) 
and compare it to the ratio obtained from Equation 
(8), Finally, if can be reliably determined 
from the nwignctic data by minimum variance analysis, 
onb may check whetlier its sign agrees with that 
required by Equation (6). This type of analysis of 
ISEE plasma and magnetometer data may be found in 
the review by G. Pas chmann (this conference). 

It may be shown (Ref. 40) that for a rotational 
discontinuity the balance of normal stresses leads 
to the requirement 

A{p^+B^/2Mq} - 0 , (9) 


The general energy equatioh is of the form 

where £ is the pressure tensor, £ the heat-flow 
vectorT* and 1 the total magnetopause current. Thi 
latter equation shows that the electromagnetic 
power • !_ associated with reconnection at the 
magnetopause may be converted into kinetic or in- 
ternal energy associated with the plasma jets 
Inside the magnetopause (first term on the left) 
or into flow work (second term) or into a net heat 
flow away from the layer (third term). All of 
these terms may be of relevance in resolving 
Heikkila’s (Ref. 42) energy crisis. 

ijn summary, it appears that an unambiguous identi- 
fication of a rotational discontinuity at the mag- 
netopause may be performed solely on the basis of 
the magnetic field data. The requirement is 
^ 0 and 0 < A^ < r. In a noniso tropic plasma 
(pjl ^) it is not in general required that 
I I 15 2 I ; changes in field magnitude can occur 
as a result of changes in ct, provided Equation (9) 
remains satisfied. The energy equation and the 
condition that the entropy cannot decrease across 
the layer do, however, place certain constraints 
on possible. changes across the discontinuity, as 
‘discussed by Hudson (Ref. 43, 44). In particular, 
when a has the same value on« the two sides of the 
discontinuity then^ in addition to (Equa- 

tion 7), one also finds, Pn 2 - Pjjj. 9 = p i 
and |^2l ~ ISfl (these latter results are obtained 
only if the heat flaw: term in Equation (10) can 
be neglected). 

7 . R0TATI0N/\L AND CONTACT 

DISCONTINUITIES: STRUCTURE 

Vlasov equilibria for one-dimensional structures 
having Bn ^ 0 are difficult to construct analyt- 
ically because the tv^o generalized tangential 
momenta, Equations (4) and (5) , now contain vector 
potential components which depend on the tangen- 
tial coordinates x;j^ and xo. Thus any distribution 
function containing these momenta will lead to 
equilibria that are not one dimensional. Addi- 
tionally, the construction of equilibria contain- 
ing a net flow across the sheet require the use 
of a constant of motion C(v 3 ) which is linear in 
the normal particle velocity V3. In geuiral such 
a constant is not readily available. 

Exact symmetric hot-plasma equilibria with Aj, = tt 
and B 2 = 0 have been constructed numerically by 
Eastwood (Ref. 45) 'with application to the geo- 
magnetic tail current sheet in mind. These solu- 
tions may perhaps be viewed as rotational dis- 
continuities in the firehose limit (a » 1) but 
more appro.priately they should be classified as 
contact distontinuities / 

‘Hot-plasma non-linear wave solutions exist 
(Ref. 46-4S) corresponding to infinite wave trains 
of circularly polarized large-amplitude electron 
and ion whistlers. For this case a simple con- 
stant of motion C(v 3 ) does exist. 

Cold plasma solutions with a net flow across the 
layer yield either infinite wave trains or elec- 
tron-polarizcd solitons (Kef , 49-52) , the latter 
Involving a net rotation of the tangential field 
of Aj. « 2tt. These results raise the fundam 
question whether laminar time-indepencrent equilib- 
rium structures of the rotational discontinuicy 


type do in fact exist. 

Finally. > analytical cold -plasma Vlasov equilibria 
for asymmetric contact discontinuities have beep 
found recently (Ref. 53). In these solutions, a 
C 9 ld plasma beam is incident from Che magnetosheath 
aiid is magnetically re flee ted at the current layer, 
njie existence of such equilibria serves to emphasize 
tlie possible importance of particle reflections in 
raagnetopause structures with 0 . 

A question of considerable interest is the sense of 
pollarization of magnetopause rotational disc^utinu- 
itjics. Bnt;h the electron and the ion polarizations 
h^ve been seen, but with a preference for the former 
(I^ef. 4, II). First-order orbit theory of the ro*- 
tationai discontinuity (Ref. 54) indicates that only 
the electron polarisation should occur if the current 
layer is sufficiently thin. A more recent simple 
model (K,ef . 13) suggests that Che electron polari*^ 
zation should be expected in layers that are suf- 
fltiently thin so that electrons, but not ions, 
flpw across the layer by sliding along the magnetic 
field lines. In this model, the electrons provide 
the fie Id -aligned current that causes the ro'Cation 
of 1 the tangential field, while the ions provide ‘an 
additional current that is more uniformly spread 
over the width of the layer. This simple idea was 
usqd in Reference 13 to account for the circular 
potation of the tangential hodogram in Figure S, In 
general, for the earth *s magnetopause, the electron 
polarization corresponds to a tangential hodogram 
with negative B 2 when B 3 is positive, as in Figure 
8 ; with positive B 7 when B 3 is negative, ns in Fig- 
ure 7. For Jupiter’s magnetopause the electron 
polarization yields B 2 < 0 for B 3 < 0 and B 2 >0 
for B 3 > 0 . 


boundary he tween magnotoshcath-llke (”1 J on ’ s roars 
and magnetosphe rebuke (’’chorus”) KIF noise, with i| 
occasional indications of emissions characteristidi I 
of the magnetopause itself. Kecentiy, Ciurne.tt qt 
al. (Ref. 17; also, this conference) have found 
evidence of enhanced noise levels in tiie magneto- 5 | 
payse and boundary layer at frequcncios from a fcwjtsf 
Hz up to well beyond the electron plasma frequency, 
with power generally decreasing rapidly with in- 
creasing frequency (except for the occasional ap^ 
pearance of narrow-band electrostatic amissions 
near the icloctron plasma frequency) . 

Examination lof high time resolution magnetometer 
data at the imagnetopause iudicates that most of 
the fluctuation power occurs at frequencies below 
the range measured by Guinett, and probably below 
tlie piroton gyrof requency as well. A substantial 'I 
contributioiV: to this low f requency noise may come ^ 
from radial oscillations in the magnetopause posi- 
tion and from waves on the magnetopause surface 
‘ (see review byD. Southwood, this conference). r| 
However, the possible importance of kinetic Alfven* ^ 
waves (Ref. 60 ) a nd o f the tearing mo d e n e e d s to 
be examined in detail. In the following paragraphV| 
we comment briefly on the signatures of the latter f 
at the magnetopause. ^ * 

The tearing mode may in principle be detected from J 
a single satellite, as a result of the convection | 
of tearing islands along the magnetopause. The ' 
resulting path of the satellite relative to the 
islands is shown by the slanted dashed line on the 1 
lower left in Figure 10. For a convection speed i | 
of 200 km/s and a fastest growing wave length 
L = 4 ttA, ZZ - 200 km being the characteris tic width 

of the magnetopause, we find a frequency of 0.16 Hi 

: I 


8 . FLUCTUATIONS 

Table 2 contains a list of relevant natural frequen- 
cies for conditions typical of the magnetopause and 
boundary layer plasma, 'nieoretical suggestions for 

Table 2. 


Electron plasma frequency 
Electron gyro frequency 
Proton plasma frequency 
Lower hybrid frequency 
Proton gyrof requency 
Tearing convection frequency 

h « 15 cm“^; B » 35 nt 


fpe = 34.7 X 10 3 Hz 
fee “ 978 Hz 
fpi =810 Hz 
fLH « 22.8 Hz 
fpi = 0.53 Hz 
ft =0.16 Hz 


the generation of noise at these frequencies include 
the clcctron-cyclotroa drift instability, (fca; > 
55) the ion acoustic instability (fpil Ref. 56), the 
iower hybrid drift instability (fLil^ v ^57) , the 
tiu’o-s trennv or current-driven ion-cyclotron instabil- 
llties (fci> 58, 16) and the tearing instability 

(hi Ref. 59 ) . The two former require current-layer 
widths of the order of the electron inertial length 
(Xq ^ 1.4 km) which Is unrealistic; the latter may 
operate when the width is of the order of the ion 
inertial length (\;j^ » 59 km) or more. 


Experimental information concerning electric and 
magnetic oscillations in the magnotopause remains 
scarce. Ne ugob a tie r e t al . ( Ref. 9 ) a nd F a i r f 1 el d , 
(Ref. 16 ) l\ a ve re po r t o d 1 0 w- f r 0 q ue n cy m agn e tic 
oscillations near the pro con gyrof requency which 
they interpret ns ion cyclotron waves. At higher 
frequencies (Ref. 9) the magnetopause appears as a 


In the local coordinate system (Xr >S z) shov^m in " ^ 
Figure 10 the tearing mode produces no fluctuations 
in By. From the condition ^ = 0 one may furth I 
dqduCQ that, away from the islands themselves, the j 
ratio of the fluctuation amplitudes in and B^ 
should be about L/2il - 2ir. In other words, the 
prediction is that tlie largest fluctuations should | 
occur in Bg, much smaller ones in and none in I 

By. The magnetopause .crossing in Figure 11 (Ref . 

1 6 ) d i s pla y s t h r ee s ub s egme nt a o f p s c il la t i o n s 
having essentially this property, the peak power f | 
occurring at a frequency of O.ll Hz. ^ | 

If the magnetopause contains a single mode (which , , 
was not the case in Figure 11), then the tearing i | 
iDode oscillations also display a characteristic ^ | 
polarization pattern as illustrated on the right 
in Figure 10. ‘ . , 

, ■ ^ ■ ."‘I 

This discussion indicates that a systematic search, j 
for tearing signatures at. the magnetopause may be 
a worthwhile enterprise . A first step in this . ^ 
direction may be found in Reference 61. j 

f ■ 

9. CONCLUSION 

In sumnmry, the local electromagnetic structure of ; 
the tnagnetopauae appears to be highly time var lab 1 i 
ranging from nearly laminar states to highly tur- 
bulent ones. To a cco un t f 0 r t h i s situ a t i 0 n , 0 ne 
may imagine a magnetopauso surfuce cqnslstiag of I 
many intermixed laminar and turbulent patciics, of | 
various sizes and in various proiTortions, wluxh 
are being convected along tlie surface away from . 
the subsolar region. Some of the patches may con- J 
tain persistent normal field components, and it 
remains pos.sible that occasionally the main featurfes 




of Dungey's classical reconnection model appear at 
the front side maRnctopanse . The field fluctuations, 
at all frequencies, in the maaootopause and boundary 
layer provide a promising arq,Vof further research, 
because these fluctuations contain essential infor- 
mation about local micro- and macroprocesscs. 
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Figure 10 


Signatures of the tearing mode. Upper 
left-hand diagram shows the perturbation 
vector potential A. On the lower -left, 
the correspondiag magnetic island struc- 
ture is shown. Slanted dashed line 
represents satellite trajectory relative 
to the convectirig islands. On the right, 
expected polarization properties of the 
perturbation magnetic field for different 
X values. 
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Abstract 

After a review of the time sequence of crossings of Jupiter's 
magnetopause by Pioneers 10 and 11, as identified in the plasma 
and magnetic field data, the detailed magnetic structures of a 
number of crossings are presented. Normal vectors are determined 
by minimum variance analysis, and the behavior of the normal and 
tangential field components is discussed. The normal vectors are 
compatible with a blunt model of Jupiter's magnetosphere. The 

i ■' . , ■ . ’ ; 

normal magnetic field component is usually small. The Observed 
structures and apparent thicknesses are compared to those of the 
earth's magnetopause. 


♦This report is based on a paper presented at S. Chapman Conference 
on Magnetospheric Boundary Layers, June 11-15, 1979, Alpbach , 
Austria. 


1 . Introduction 


It is the purpose of this report to discuss the structure of Jupiter's 
ipagnetopause and to compare it with the terrestrial magnetopause. 

Data from the JPL helium vector magnetometers on board Pioneers 10 
dnd 11 (Smith et al,1976) will be used along with supporting evidence 

I ' ' . 
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from the AMES plasma analyzer' (In triliga tor and Wolfe, 1976). The 
terrestrial magnetopause structures shown for comparative purposes 
were obtained from the GSFC flux-gate magnetometers onboard OGO-5 
(J. Heppner, principal investigator). 

2. Pioneer 10 Observations 

Figure 1 shows the time sequence of magnetopause crossings during 
the Pioneer 10 encounter with Jupiter's magnetosphere. During the 
inbound leg of the trajectory, three’ magnetopause crossings were 
identified in the plasma data: November 27, 1973, at 20:35 UT 

(ground receipt time) ; December 1 at 03:18 and 14:20 (Intriligator 
and Wolfe, 1976) . In the magnetic data the corresponding crossings 
appeared at 20:32, 03:20, and 13:56, respectively, indicating that 
the plasma was present somewhat inside the magnetopause. Several 
additional crossings, not apparent in the plasma data, may have 
occurred, as shown in the figure, in this report, we shall discuss 
the initial sequence of crossings (November 27 19:25, 19:40, 20:32) 
an<3 the two crossings on December 1 at 03:20 and 13:56. Also, the 
final outbound crossing on December 14 at 19:36 will be examined. 

The other outbound crossings, identified in the plasma data, have^^^^^^^^ 
yet to be uniquely associated with corresponding magnetic field^^ ^^^ ^^^^^^ 
structures. 


Figure 2 shows a time record of the magnetic field during the initial 
three crossings. The field. is represented by its magnitude (dashed 
line) its latitude angle 0 in SJ coordinates (solid line) and its 
longitude angle (dotted line) . The field outside the magnetopause 
is due slightly north. As the magnetopause is entered, the field 
first swings in the "wrong" direction, i.e. further north, and then 
turns south over a time span of 3 minutes to take on the direction 
typical of Jupiter's magnetosphere. Hov/ever, the magnetopause starts 
moving back over the spacecraft almost immediately, and the field 
returns to the magnetosheath direction. The return is much slower 
and requires ~ 11 minutes. The excursion into the magneto sheath is 
brief and the third crossing is again a. rapid one (•” 3 minutes) . 

This time it is evident that an entry into the magnetosphere proper 

■i ■''' ■ 

has occurred because the inner edge of the plasma boundary layer, 
mhrked by an abrupt field increase, is crossed at 20; 41. This time 
should replace the approximate time 20:35 given by Intriligator and 

t ■ 

Woilfe (1976) . Firom the change in field magnitude, with an assumed 
tdmperature T = 5 x 10®°K, we conclude that the plasma density in 
the boundary layer was n ^ 1.3 cm^ and that the beta value was 

■ 1 ■ ■ ■■.M 

3 i> 1.6 (6 is the ratio of plasma to magnetic pressure) . 

\ ' ■ 

{ ■... . . 

DOi the two first "sorcalled" crossings really represent traversals 
of the magnetopause or are they magnetic-field structures located 
in the magnetosheath outside the magnetopause? Figure 3 illustrates 
that both the normal vector orientations and the duration of the 
crossings are compatible with the interpretation of the time; record 
as three distinct crossings, the multiplicity being caused by wave 
motioh on the magnetopause , The normal vectors for the three cros- 
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sings, determined from minimuin. variance analysis, are given in terms 
of their solar- jovian X, Y, and Z components in Table 1. 

Table 1 Normal vectors for the Pioneer 10 magnetopause 
crossings on- Nov 27, 19;25, 19:40, 20:32 UT. 

0JJ is the latitude angle of the normal vector, 
and At is- the duration of the crossing. 


Crossing 

% 

Ny ■ 

“z 


At 

(sec) 

1 (19:25) 

.785 

-.574 

i -:‘235 

|. -13.6 

180 

2 (19:40) 

.855 

-.509 

-.104 

- 6.0 

700 

3 (20:32) 

.943 

-.224 

-.247 

-14 .3 

190 


All three vectors point approximately in the direction expected 

I 

for a forenoon crossing occurring slightly below the equatorial 

I 

plane (longitude 5 ^ = 325°, latitude A = - 6 °) . The monotonic in- 
crease in Ny time progresses may be the result of a very long 
wave length ondulation of the magnetopause surface. The change 
in Ng is not monotonic. It is much greater than estimated errors, 
and it is compatible with the wave model shown in figure 3 on the 
left. The wave travels south and is responsible for the triple 
crossing. of the magnetopause. Assuming -the wave to sit still and 

the spacecraft to move north along the magnetopause instead, with 

i 

the wave speed, v, and across it with the spacecraft speed (v ~ 10 

I 

km/s) , we see that the first crossing should have a relatively 

I 

lairge negative Ng and a short crossing time At. The second cros- 
sing should have a smaller negative N^ and a long crossing time, 

and the third one again a larger negative N and a short crossing 

'Z " 

time. As shown in the table, this is exactly what occurs. 
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The right-hand part of figure 3 shows that. the model can also be 
used to estimate the actual magnetopause thickness h^p. Assuming 
small angles, three equations relating h^p/v to the angles 0 „, 6 


and 6 are obtained. 


N' 


MP" 


h^/(v4t.) = - (- 0 ^- 6 ) 


hjjp/CvAtz) = 


” ®N3 


1 

"MP ^ 

(-e_ - 6)J 


0JJ2” 


h„„= 32.lv ~ 3850 km 
MP 


h™ = 37. 8 v 4540 km 
MP 


Here 6„ = sin ^N„ is the measured tilt of the normal vector, 6 ~ v /v 
is the angle associated with the skew traversal of the spacecraft 
aqross the .wave, and G^^p is the average tilt of the magnetopause. 
Also, At is the time duration of a crossing, measured fro|n peak 
northward to peak southward field, or vice versa. From the three 
equations, by elimination of ” 5)/ two estimates for h^p are 

obtained. Given the crude nature of the model, the agreement be- 
tween the two estimates is satisfactory. Assuming wave speed of 
V — 120 km/s, the thickness is h^p ~ 4200 km, i.e., about ^ 

(R.^ . = 1077 km for a 500 Volt particle in a 3 y field) . On the 
pther hand, for n = 1.3 cm~^ the ion inertial length = (m^/y^ne^) ^ 
200 km so that h^^p ~ 21 A^. This •result supports the view that the 
relevant scale length for the thickness of a tangential-discontinuity 
type magnetopause is Rj^^ rather than A^. As a comparison, for a 
typical terrestrial magnetopause crossing the thickness may be 


b,,T. ~ 400 km while R^ • = 90 km, A. = 4 0 km (e = 500 Volt, B = 35y, 
MP Lx X 

n = 30 cm” . Since the jovian magnetosphere on this occasion was 
more than 100 times larger than the terrestrial one, it is also 


-5- 


evident that the magnetopause was not proportional to the size of 
the magnetosphere. 

The thickness hgj^ of the plasma boundary layer may also be estimated. 

i 

For V = 120 km/s one finds hgj^ ~ 6600 km. The ratio hgj^/h^p is 
comparable to the terrestrial values at corresponding locations on ' 
the magnetopause. Finally, the average tilt of the magnetopause 
relative to the Z axis may be calculated. For v ~ 120 km/s it is 
®MP ~ the latitude of the crossing (A = - 6 °) this result 

corresponds to a relatively blunt magnetopause shape. 

The sequence of three crossings can also be accounted for by assum- 
ing periodic radial motion of the magnetopause. In such a case too, 
the middle crossing is expected to have a longer penetration time. 

In a simple model with magnetopause motion either inwards or outwards 
with speed Vj^p one may then use the observed penetration times to 
find Vjyp = km/s, h^^p = 5000 km. The nature of the model is such 

that it overestimates the thickness h,,». The conclusion that R„ . 

MP ill 

rather than X-, or a constant fraction of the size of the magneto- 
sphere, is the appropriate scale length remains vaJid. 

We shall now examine the structure of the three crossings more in 
detail. Figure 4 shows a polar plot of the tangentlaJ .Held com- 
ponents (B 1 B 2 ) on the left, and of the normal field component B 3 
on the right, for. the first of the three crossings. Note that the 
tangential component B 2 remains negative in the crossing. Note 
also that the average normal field B 3 is nearly zero (although the 
ring-shaped structure in the middle of the crossing apparently has 
an orientation different from that of the main magnetopause) . The 
structure is, on the whole, consistent with a tangential disconti- 


n^ity. The second crossing, shown in figure 5, displays a com- 
plicated turbulent structure. In reality it may not have been 
more turbulent than the first one. Rather it is the much longer 
crossing time in the second crossing that allows one to see more 
of the turbulent changes in the structure. Note that B 2 remains 
negative as in the first crossing, and that the average normal 
field component B 3 remains essentially zero. The third crossing, 

! . ^ I . 

i : i 

shown in figure 6 , is again relatively rapid and displays less 
turbulence. The component B 2 is negative and B 3 ~ 0, as in the 
previous cases. The internal consistency (B 2 < 0 , B 3 ~ 0) of the 
three observed structures supports the interpretation of the outer 
two structures as true magnetopause crossings. 

Figure 7 shows a comparison of the time record for the last of the 
three Jupiter crossings above and for a terrestrial magnetopause 
crossing observed by the GSFC magnetometer on board OGO-5. There 
is a striking qualitative similarity between the two records. 

First, the inner edge of the plasma boundary layer is seen clearly 

i ■ ■ , . ^ ■ ■ ' ■ ■■ ■ 

in both records (although in the terrestrial case with T = 10®°K 

the boundary layer density was n ^ 225 cm”®) . The field direction 
changes can be compared only after one of the two 6 traces has been 
mirror imaged in the time axis to compensate for the fact that in 
the equatorial plane Jupiter's field is due south, while the earth' 

I 

field is due north. When this has been done, it is seen that both 
records display an initial direction change in the "wrong" sense 
(ais mentioned already in connection with figure 2 ) near the outer 
edge of the magnetopause . This behavior is seen frequently at the 
earth's magnetopause and can probably be accounted for in terms of 
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a finite gyroradius effect. Glearly, Jupiter’s and the earth's 

I ' ; I ^ ■ ' 

magnetopause are occasionally remarkably similar. 

Figure 8 shows, in polar form, the magnetopause structure during V 
the next exit of Pioneer 10 from the magnetosphere. In this case 
too there may have been a sequence of three crossings, but the 

case for such an interpretation is less convincing than for the 

I 

earlier sequence. Thus, only the innermost crossing is shown here. 

It had B 2 > 0 and contained a fairly large amount of turbulence. 

On the average the normal field component B 3 remained approximately 
zero. The penetration time was about 10 minutes. 

Figure 9 shows the polar plots for the last crossing on the inbound 
leg of the Pioneer 10 orbit. Its structure was complicated, but 
again one notes a tendency for the normal component to average to 
zero. The penetration time was about 5 minutes. 

Finally, figure 10 shows the last outbound crossing observed by 
Pioneer 10. This was a very rapid crossing with a penetration time 
o^ less than one minute. The structure was remarkably simple and 
contained a portion that we identify as a rotational discontinuity* 
The normal magnetic field component is - 0. 82y ± 6.26 y; in other 
words, it is significantly different from zero. The negative sign 
is what is expected in an open magnetosphere model for a crossing 
south of the reconnection line. The spacecraft latitude was small 
(A = + 8 . 2 °) so that a location south of a hypothetical reconnection I 
line is entirely possible. The sense of rotation of the tangential 
field vector is remarkable. In a model discussed by Sonnerup and 
Lediey (1979) tlie field rotation in a thin rotational discontinuity 
is produced by the electrons in the magnetosheath plasma flowing 
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across the discontinuity. These electrons move adiabatically and 
essentially along the field thus causing a field-aligned current. 

In the present case, this model works only if the flow is from the 
magnetosphere into the magnetosheath, rather than vice versa. There 
is anotlier piece of evidence that suggests the same thing, namely 
the field magnitude increase at the outer edge of the magnetopause, 
followed by a small field rotation. In the OGO-5 crossing of the 
earth's magnetopause shown in figure 11 (Sonnerup and Ledley, 1979), 
the same structure is present, but now at the inner edge of the 
magnetopause instead. This feature has been discussed by Sonnerup 
and Ledley in terms of a slow expansion fan, as proposed by Levy 
et al (1964) , followed by a field direction change of the type 
discussed by Parker (1967 a, b) . The general Gonfiguration shown 
in figure 11 is compatible with reconnection, with flow from the 
magnetosheath into the magnetosphere. Does the Jupiter crossing 
indicate reconnection with a reverse flow direction? We do not 

insist that this must be so. But the remarkable similarity between 

i _ ^ ' ■ ' 

the Jupiter crossing and the terrestrial one indicates the importance 
of developing a better understanding of rotational field structures 
at the magnetopause. 

3 . Pioneer 11 Observations 

We turn now to a brief discussion of the Pioneer 11 magnetopause 
crossings. The time sequence of possible crossings on both the in- 
bound and outbound legs is shown in figure 12 for plasma as well 
as magnetic field data. The following comments should be made. 

The first inbound crossing at 02:08 on November 27, 1974, is gradual 
and poorly defined in the magnetic data. A second crossing ( from 


the magnetopause back into the magnetosheath) is found at 08:19. 

This crossing is well defined (although of 20 - 30 minutes duration) 

..I ■ ■ ■ . ■ i 

but, paradoxically, it occurs substantially before the crossing 

time indicated by the plasma data. A reexamination of plasma and 

magnetic field data is underway for this case. The final entry 

into the magnetosphere at 13:30 on November 29 was in all likelihood 

preceded by a brief entry of 1 hour duration. We shall discuss only 

the final crossing (figure 13) , 

On the outbound leg, the situation is initially complicated in that 

i , ; . 

there is much magnetic field activity in a region identified via 
the plasma results as being the magnetosphere. Included is one 
field change at 06:42 on December 6 which has many, of the features 
of a regular magnetopause crossing (figure 15). The crossings at 
08:47 (figure 16) and at 19:07 on the same day are well defined 
but involve only small field direction changes . The final outbound 
crossing is also well defined i (figure 17). 

Figure 13 shows the polar representation of the magnetic field during 
the last crossing on the inbound leg. The duration of the crossing 
was 20 minutes, yet it does not display a grossly turbulent structure 
Rather, small scale turbulence seems to occur in patches. This turbu 
lence has considerable structure with alternating regions of circular 
and linear polarization. The frequency range centers around 0.05 Hz. 
One possible interpretation is in terms of tearing bubbles being 
Gonvected past the spacecraft (Greenly and Sonnerup, 1979). The 
gross structure of the crossing is that of a rotational discontinuity 
The normal field component .is + 0.71y ± 0.43y,, and the relationship 
between the sense of rotation an^ the sign of the normal field com- 



ponent is that expected from the model by Sonnerup and Ledley (1979), 
p this time assuming flow across the magnetopause from the magneto- 
sheath into the magnetosphere . 


i'i For purposes of comparison, figure 14 shows an OGO- 5 crossing of 
I ; the earth * s magnetopause which displays essentially the same features 
ais the crossing in figure 13, In this case too, an interpretation 

• : ; - ■ i ■ r: 

fi ; in terms of a rotational discontinuity with flow from the magneto- 

sheath into the magnetosphere and with convecting tearing structures 

;■ ■ ; : ■ : • r • ^ 

• appears reasonable (Greenly and Sonnerup, 1979) . 



Figure 15 shows a magnetic field structure observed on December 6, 

1974 at 06:42 during the outbound portion of the Pioneer 11 trajectory. 
The sti'ueture occurred! in a region identified via the plasma data • 
as the magnetosphere. However, it is hard to believe that it is 
not a bona fide magnetopause crossing. It has the basic features 

• ' “ f ■ \ ■ i 

of a rotational discontinuity. The normal field component, which 
appears positive on the slide, is in fact found to be slightly neg- 
ative (- Q.20y ± 0.26y) when the minimiim variance analysis is per- ^ ^ 

fbrmed on the main field rotation data segment alone. Thus, the 
main rotation is in the right sense, assuihing plasma flow from the 
magnetosheath into the magnetosphere. Fuithermore, the magnetic 

: ^ ■ - , ■ :• ■; , , : 1: - I ■ 

field incirease and final field rotation at the inner edge of the 
crossing appear to be the exact counterparts of the features shown 
in figure 11 for an OGO-5 crossing of the earth's magnetopause. 

Figure 16 shows a time record of the butbound magnetopause Crossing 
at 08:47 on December 6 , 1974 . Here the field direction change is 
quite small; yet it is clear that the magnetopause' is an effective 
barrier for the plasma. With T = 5 x 10®°K, the field magnitude 
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change indicates a density change An ~ 3. 9 cm*’^ across the magneto- 
pause. The normal magnetic field component is near zero for this 
crossing. 

: ■ ■ 'I 

The last outbound crossing by Pioneer 11 is shown in figure 17. It 
has a duration of ~ 7 minutes, and the structure is complicated. The 
BiBi diagram suggests that the magnetopause attitude may have changed 
during the crossing. For this reason the Y component of the normal 
vector N^y (which should be ~-0.12 given the location: 353°, 

A = 31.2°) is unreliable . 

4 . Shape of Jupiter's Magnetopause Surface 

Figure 18 summarizes information concerning the gross shape of 
Jupiter's dayside magnetopause, extracted from the normal vector 
calculations. The angles 6^^ and^^^^ are the latitude and longitude 
angles of the normal vectors, while A and are the corresponding 
quantities for the spacecraft position veGtor. The 45° line in 
the two figures helps one compare the data to a model consisting^ 
of a centered sphere . 

The most striking feature of the diagrams, when compared to cor- 
responding ones for the earth's magnetosphere, is the large spread 
in the data. Presxmably, this larger scatter is the result of a 
jovian magnetosphere that is far more "floppy" than the terrestrial 
one. The data do not provide support for a sharp-nosed (in the^ 
noon-midnight meridional plane) model of Jupiter's magnetopause. 

On the average, the shape may be about as blunt as a sphere . 
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5 . Conclusion 

We have shown that there is a striking similarity between the mag- 
netic structures of the jovian and the terresti:ial magnetopause. 

Our estimate of the magnetopause thickness supports the view that 
normally the ion gyroradius R_ . is the relevant scale size, rather 
than the ion inertial length, or some constant fraction of the mag- 
netosphere radius. Thus, one should not expect to find thick- 
nesses less than one or two The possibility that the magneto- 

pause thickness sometimes is many . and perhaps even proportional 

ijX 

to the size of the magnetosphere cannot be excluded, since such 
extremely thick structures have been found at the earth's magneto- 
pause on rare occasions. 
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8. Figure Captions 

Fig. 1 Pioneer 10 crossings of Jupiter's magnetopause in 1973, as 

recorded by the JPL magnetometer and the AMES plasma analyzer 

Fig. 2 Time record of magnetic field observations during the first, 
triple, inbound magnetopause crossing by Pioneer 10. Theta 
and phi are the SJ latitude and longitude angles of the 
magnetic field vector, respectively. Jupiter's magnetosphere 
has 6 ~-90°. The field magnitude is measured in nT (lnT=lY) 

Fig. 3 Schematic of magnetopause wave motion model to account for 
the normal vector variation and the penetration times for 
the three crossings in figure 2. 

Fig. 4 Polar representation (hodograms) of the magnetic field 
during the first of the three magnetopause crossings in 
figure 2. The field components Bj and Ba , measured in nT, 
are tangential to the magnetopause surface, with the former 
due approximately north and the latter approximately west. 
Thus B 1 is negative in Jupiter ' s magnetosphere . The 
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component B 3 is perpendicular to the magnetopause and is 
positive when directed away from the planet. The matrix 

in the lower left-hand corner of the figure gives the 

i ■ ' ' ^ 

components of the corresponding unit vectors Ni , N 2 » and 
Ns , along these three directions / in SJ coordinates. N 3 
is the magnetopause normal vector. 

Fig. 5 Hodogram for the second crossing in figure 2. 

Pig. 6 Hodogram for the third crossing in^ f^^^^ 

Fig. 7 Comparison of time record of the third crossing in figure 2 
(top) with a terrestrial magnetopause crossing by the sat- 
ell ’ te OGO-5 (bottom) . 

Pig. 8 Hodogram of Pioneei 10 inbound magnetopause crossing on 
December 01/ 1973,. 03:20 UT. 

Fig. 9 Hodogram of Pioneer 10 final inbound magnetopause crossing 
on December 01, 1973, 13:56 UT. 

Fig. 10 Hodogram of Pioneer 10 final outbound magnetopause crossing 
on December 14, 1973, 19:36 UT. 

Pig. 11 Hodogram of OGO-5 crossing of the earth's magnetopause. 

Note that Bx > 0 in the earth's magnetosphere while Bi < 0 
in Jupiter's magnetosphere. Compare to figures 10 and 15. 

Fig. 12 Pioneer 11 crossings of Jupiter's magnetopause in 1974, as 

recorded by the JPL magnetometer and the AMES plasma analyzer 

Fig. 13 Hodogram of Pioneer 11 final inbound magnetopause crossing 
on November 2,9, 1974, 13:30 UT. 

Fig. 14 Hodogram of OGO-5 crossing of the earth's magnetopause. 
Compare to figure 13 . 
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Fig. 15 Field structure observed by Pioneei: 11 on December 06, 1974, 
06:42 UT. Compare to figures 10 and 11. 

Fig. 16 Time record of Pioneer 11 outbound magnetopause crossing 
I on December 06, 1974, 08:47 UT. 

Fig. 17 Hodogram of Pioneer 11 final outbound magnetopause crossing 
on December 08, 1974, 01:22 UT. 


Pig. 18 Shape of Jupiter's magnetopause surface* The angles A and 

y^are the SJ latitude and longitude of the spacecraft during 
p|enetration of the magnetopause . The angles 9^^^ and are 
the SJ latitude and longitude of the magnetopause nornjal 

i ■ ■ ^ ^ ' 

vector Ns, determined from minimum variance ahalys is. The 
terms "sharper" and "blunter" refer to the magnetopause 
shape relative to a centered sphere (= 45° line) . 
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TEARING MODES AT THE MAGNETOPAUSE 

J. B. Greenly* and B. U. 6. Sonnerup 
Dartmouth College, Hanover, New Hampshire 03755 

Abstract 


This paper examines the possible occurrence of collisionless 
tearing modes in the dayside .magnetopause. The expected magnetic 
signature of tearing is obtained from existing theory. Magnetometer 
data from one terrestrial magnetopause crossing, and one crossing of 
Jupiter's magnetopause are then examined in detail. Magnetic field 
oscillations are found in three subsegments of the terrestrial cross- 
ing at a frequency of 0.1-0. 2 Hz and with peak amplitudes of 5-10 
nanotesla (nt), and: in one segment of the Jovian crossing at .05-.1 Hz 
and vjith 2 nt amplitude . . The frequency range as well as the orienta- 
tion of the magnetic field perturbation vectors agree with a model 
in which tearing-produced magnetic islands are convected past the 
satellite with the plasma flow in the current layer. In both casdsr^ 
the magnetopause structure was of the rotational discontinuity type 
with a nonvanishing normal magnetic field component . Hence, if the 
tearing structures were active, i.e., growing, at the observation 
site, ion tearing must be invoked. But it is also possible that the 
structures were passive, consisting of "debris" from active tearing 
elsewhere on the magnetopause surface, this debris being convected 
along the magnetopause past the observation site. 
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1. Introduction 

The structure of the magnetopause current layer has been the 
subject of many Investigations, both theoretical and observational 
(see, je.'g. Willis, 1971, 1975, 1978; Sonnerup, 1976; Sonnerup andLedley, 
1979 ). The most obvious and striking characteristic of the magneto- 
paluse as observed during spacecraft crossings is its variability. 
Steady-state one-dimensional models may help explain certain average 
basic features of the current layer, although their success to date - 
has been limited. But, in order to account for observations in 
detail, one must include temporal as. well as two and three dimensional 

spatial variations , * 

It is our purpose here to examine the possibility that the col- 
lision-free tearing mode may be operative in the magnetopause as 
suggested by Galeev and Zeleny (1978) . To do so we first describe 
some of the basic magnetic field signatures associated with tearing 
and then compare these signatures with substructures in the earth ’ s 
magne t op aus e , observed by the Goddard Space Plight Center magnetometer 
onboard OGO-5, and in Jupiter’s magnetopause, observed by the Jet 

Propulsion Laboratory vector helium magnetometer onboard Pioneer ‘1. 

! ^ ■ . ■ . " ' ' ^ ■ , ■ ; , - , ■ , 

2. Tearing Mode Signatures 

The basic unperturbed magnetic field geometry assumed in tearing 
mode analyses is of the form: 

w'-;; : ■ i'' 

■ : , ' ■ ' ' ^ ! ' 

£ ' Mox iyV _ (1) - 

Here £. is the characteristic width of the region of the magnetopause 
in which tearing occurs. It is referred to as the shear length. 

Also, B = B is the constant magnetic field component perpendicular 


to the magnetopause i B = B _ a constant field component in the plane 

V ‘ i ' ■ ' 

of the magnetopause and parallel to the current. iPinally, B„ is the 
reversing field component which takes on the values when jx] >> Z 

The tearing mode perturbation of the above field has a vector 
potential of the form 

P A = iy (2) 

wheire y and k are the growth rate and wave number,, respectively, and 
A(x) is an unknown function of .x, the coordinate perpendicular to 
the layer, to be determined along with the growth rate from the 

Vlasov-Maxwell ejquations; When such a field is superposed upon the 

unperturbed field,, the resulting field line configuration is as shown 
in Figure 1, The behavior of A(x) is qualitatively as shown in the 

figure. Note that, apart from the phase factor, B ~ kA(x) while 

' 

Bg, ~ dA/dx . 

Detailed expressions for A(x) and , for the growth rate y have 
been derived under a variety of assumptions, none of which are 
entii’ely applicable to the present situation. However, the following 
two results are of Interest : 

(i ) In a colli Sion- free plasma with the unperturbed magnetic 

field component perpendicular to the current layer, equal to zero, 
Galeev and Zeleny (1978) found a linear growth rate 

''' ''«V. B -V ^ 

^ (3) 

■ ---■ /ir „ oy : , ' l' 

In this formula v, , is the electron thermal velocity and!' X = (m /y ne 

I I ■ G 0 O 

i 0 i 

is the electron inertial length (m^ being the electron mass and n the 
particle density). This formula describes a purely growing mode , 


! 7 ^’ 




I: 


-^1- 


dominated the electron dynamics in a region of width comparable 
to X around x =0, and therefore Referred to as electron tearing. 

G f ' ^ ! 

' ,r' . i ' i i- - . I i. 

In arriving at Equation (3) the assumption B^/B >> (R. .'/£) ^ was 
made, where the ion gyroradius in the field This assump- 

tion is at best marginally satisfied at the magnetopause where, in 
usual circumstances, is cpmparaible to B^^ and Z to Even 

more seriously, the normal magnetic field component B at the mag- 

O A 

netopause is often sufficiently large so as to suppress the electron 
tearing mode , 

When Y is evaluated from (3) at the earth's magnetopause, where 
®oz ~ ®oy* ^ “10 cm“®, kT^ = 25eV, and Z = 50 km, one obtains a 
maximum growth rate = 1/50 sec“V for a wave length X = 27r/k = 

irioLX 

= 600 km. • 


i-~ 

li 




I •/.. 


(il) In a collision-free plasma with B^_ / 0 but B = 0 Galeev 

O A oy 

(1978) finds i 


Y - + 1^) (1 - k^Z^) 

^ /^Z ^ f ^ i 


(» 


In this^^^^^^^ m^^ referred to as ion tearing, the ions provide the es^ 

sential dynamics and it is required that B^^/B^^ < y m^/eB^^. (A 
lower limit on B./B also exists; see Gaieev (1978) , ) The asSump- 

O A O Z 


tion B = 0 is usually not satisfied at the magnetopause. 

oy j 

In a proton electron plasma at the earth's magnetopause with 
kTj^ = 170 eV, B^^ = 35 nt, T^/T^ = 0.15, Z = 50 km, and k = 1/2Z one 
finds from (4): Y = 1 sec“V, and B^^/B^^ < Yi^j^/aBQg = 0.28. 


Neither of the two cases discussed above is strictly applicable 
to the magnetopause situation. For example, the influence of velocity 


r= 

il': 
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shear in the magnetopause is not Included in either case. Addition- 
ally, for magnetopause crossings with a nonvanishing normal field 
component, plasma flow across the magnetopause is likely to occur, 
and it is not known how such an effect influences the ion tearing. 
Killeen and Shestakov (1978) have shown, in the high-conductivity 
limit of resistive MHD tearing theory, that nonsymmetric modes may 
be strongly destabilized by diffusive flow across the current layer. 
Their analysis, however, is again not directly applicable to the 
magnetopause. Furthermore, in most observed magnetopause structures 
B is not zero and does not remain constant throughout the layer. 

if 

Thus, the results given above must be considered as illustrative 
rather than as quantitatively applicable, They nevertheless suggest 
that the following statements are likely to be relevant to magneto- 
pause tearing: 


(a) Ion as well as electron tearing leads to the formation of mag- 
netic islands which do not propagate relative to. the plasma. The 
fastest growth occurs at a wavelength of the order of where £ 
is comparable to the width of the current layer (compare Eq . (l) ) . 


(b) At least for thin current layers , the growth ratds may be suf- 
ficient to permit the tearing modes to yield finite magnetic per- 
turbations before the island structures are swept away from the front 
side magnetopause by the solar wind* For the terrestrial magneto- 
sphere, ’the time scale for the latter process is of the order of 
5-io minutes. However, the grdwth rate decreases rapidly with 
increasing layer thickness, so that active tearing may not be seen 
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in thick magnetopause structures. At any rate, tearing should be 
more easily observable away from the sub-solar point than near it , 

The results of npnlinear tearing analysis are limited at present, 
but several investigators have found, by analytical means or by 
numerical simulation (Dlckman et al, 1969; Biskamp et al,^ 1970; 
Rutherford:, 1973; Schindler, 197.^; Dralce and Lee, 1977; White et al, 
1977; Galeev and Zeleny, 1977), that single modes can grow to large 
amplitudes. One recent analysis (Galeev et al,.l978) even indicates 
explosive growth in the nonlinear stage of ion tearing. On the other 
hand, Pellat (1978) has argued that collision- free Ion tearing is not 
possible except in very special, circumstances . However, resistive 
tearing in a current sheet with a finite normal magnetic field 
remains possible, the resistivity being provided by for example the 
lower hybrid drift instability (Huba et al, 1977). On the whole, 
there appears to be no strong theoretical reason to doubt that mag- 
netopause tearing modes may occur and may grow to amplitudes of at 
least several nanotesla. A quantitative i description of such large 
amplitude tearing structures in the magnetbpause field is not yet 
available, a fact that limits the possibilities for detailed compari- 
son between theory and observed signatures (see point (e) below). 

(c) The convection of magnetic tearing structures along the mag- 
netopause can makb the tearing mode observable from a single satel*^ 
lite crossing the magnetopause. If the convection speed is denoted 

by V , a Doppler shifted frequency of f = (v k/2Tr)cos6 is obtained. 

o ■ ■ c, 
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Here 6 is the angle between the convection velocity vector and the 
z axis (the direction of the "propagation" vector k ) . For v = 200 km/s 
(el«g., Aubry et al^ 1971) and with k = l/2£ = 10“^m“* (corresponding 
to .the most rapidly growing wavelength) we find f = (0.3 cos0) Hz 
for A = 50 km. Because of the short wavelength cutoff at k = 1/5- 

(see (3) and (4)), frequencies above 0.6 Hz should not occur. 

^ i : 

The angle 6 also requires comment. If 0 = ir/2 no Doppler shift 
occurs and the tearing mode remains unobservable from a single rate 1- 
lite. For example,' if the convection velocity vector is taken to point 
approximately away from the sub solar point, 0 - it /2 occurs in the 
equatorial plane of the magnetosphere when the magnetosheath field 
is antiparallel to the earth’s field. 

The yelpcity of a satellite relative to the magnetopause may vary 
greatly because of radial motion of the magnetopause Itself. At the 
terrestrial magnetopause, relative velocities of the orcier of v =5-10 
km/sec (normal to the magnetopause) are typical. Thus, with a con^ 
vection speed of 200 km/s, the motion of the satellite relative to 
the tearing structure in Figure 1 is at a rather skew angle as indi- 
cated in the figure .During a crossing a satellite would be expected 
to transverse N = (v.cos0/v )(£k/2F) wavelengths. For v./v = 40, 

0 = 0, and k = 1/2JI we find N = 3. The conclusion is that under 
favorable circumstances, i.e., when v /v : is large and 0 is suffi- 
ciently small, tearing should be observable in the magnetic field 
records taken during magnetopause crossings , 
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(d) All tearing theories to date are such that they leave the B 

y 

component entirely unperturbed. In other words, the magnetic field 
perturbations are confined to the xz plane. This characteristic 
signature should be helpful in identifying tearing qsoillatlons during 
magnetopause crossings. It is perhaps somewhat surprising that B 

y 

should remain constant in the tearing mode. Such a result requires 
the macroscopic ion and electron motions In the xz plane to be iden- 
tical so that no Hall current loops develop in that plane. A more 
refined analysis is likely to reveal the presence of such loops and 
associated small variations in B . 

y 


(e) For a single mode, the tearing perturbations in B and B are 

X z 

out of phase by ir/2. Furthermore the amplitude ratio is 


B 


B, 


|kA| IkA 

1 dA/dx 1 Ta7T 


= kJl . 


(5) 


Since kX < 1 always, the B^^ oscillation should^ have larger amplitude 

z 


than the B„ oscillation. However j, as shown In Figure 1, B_ vanishes 

I X z 

at X = 0 and at x = ±d. An observer at fixed x position, watching 


in the z direction, would see 


the tearing structure convected 
various polarizations of the B , B_ field oscillations, as shown in 

I X z 

the figurel For |x| > d, the polarization is elliptical with [B 1 > ]B 
but for |xj £ d, a complicated polarization dependence on x occurs. 

The scale length d is Independent of mode amplitude and depends on 
k for small amplitude (linear) modes, and' d -»■ 0 as kA -»■ 1. The size 
of d for large amplitude modes at the magnetopause is not knoxm, but 
should be of the order of the Island width , as shown in the figure . 

In principle , it appears possible, then, that the complex polarization 
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changes for |x| < d could be observed. In practice the situation 
may be more complicated for two reasons. First, the x value of the 
satellite changes as it crosses the magnetopause, leading to a 

■ I : - ' , . 

relatively rapid transition from .one polarization pattern to another. 
Second, more than one wavelength is likely to be- present. Thus one 
would not expect to easily observe the polarization ellipses discussed 
above. However, since these ellipses have their major axis along 
the z direction except in a layer, probably narrow compared to Jl, 
around x = 0, one would expect the perturbations to have |B ] > |B ] 
except in that layer. 


The phase (polarization) relation between B and B_ is antlsym- 

. . : I ■ X ' Z 

i . . ■ : 

metric about X =0, A spacecraft traversing the x = 0 plane of a 
single tearing mode must observe elliptical polarization of opposite 
sense at equal distances on either side of x= 0. This antisymmetry 
would be strong evidence of having passed through x = 0 even if the 
complex behavior for [xj < d is not seen. Conversely, if at the 
beginning and end of an o)?servatlon of oscillations in B^ and B^, 
the sense of the polarization is not reversed, then the observer i 
has not passed through x - 0. The polarization antisymmetry is -not 
urique to tearing modes, but is a necessary property of any pertur- 
bation with the same symmetry as the unperturbed magnetopause field 
of Eq. (1), that is, symmetry under l80° rotation about the y axis. 

thus have a criterion for traversal of the k • B = 0 (x = 0) plane 
that is Independent of the assumed mode type. Unless this overall 
polarization antisymmetry is observed , one should not expect to see 
any of the detailed structure near x = 0. 
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(f) In actual observed magnetopause crossings, a constant or nearly .1 
constant value of B (corresponding to a magnetopause current that is 

J : R 

directed entirely along the y axis) is seldom observed. Rather the "* 

i 

B value may change in a systematic manner as the satellite passes [ 

through the magnetopause . For example, in a rotational form the 

quantity (B / + B remains approximately constant (see e.g., Sonnerup .1 
y ^ 

and Ledley, 197^). In order to apply existing tearing mode analyses ”1 
and much of the previous discussion to such crossings it is necessary 
to examine subsegments of a crossing, selected such that one magnetic f 
field component tangential to the magnetopause remains nearly constant 
during the subsegment. A local Cartesian coordinate system is then 
constructed with the y axis along the constant field, with the x axis " ] 
remaining perpendicular to the magnetopause surface, and with the 
z axis completing the right-handed orthogonal triad. It is in this 
local system that Figure 1 is to be applied. The angle 9 under point 
(c) above is then the angle between the local convection velocity .1 

(which may differ from the magnetosheath flow velocity) and the local 'i 

I 

z axis. Note also that the characteristic width Z is associated with 
the subsegment rather than with the entire magnetopause. J 

As is evident from the previous discussion, it is important for ,, 

the observational study of magnetopause tearing to have precise ^ 

information about the direction (x) normal to the magnetopause. Such ■' | 
Information is obtained from minimum variance analysis df the magneto- 

: : : ' '‘'I 

pause magnetic field vector data set (Sonnerup and Cahill, 1967). I 

It is known that, because of degeneracy of the variance ellipsoid, 
this method falls to give an accurate determination df the normal 
direction when the current in the magnetopause is unidirectional as ; I 

it is in the model in Figure 1. Indeed, if minimum-variance analysis 
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were applied to ^ crossing in which the structure is exactly as in Figure 
1,| the normal (i.e., the minimum variance) direction would be selected 
aliDng the y axis rather than the x axis. Thus, it is desirable to 
examine actual magnetopause crossings in which the current is not 
unidirectional. Rotational forms would appear ideal for this purpose. 


(g) It would of course be of interest to examine not only the mag- 
netic signature of magnetopause tearing but the particle signature 
as well. However, we have found that existing collisionless tearing 
theory has for the most part not reached a stage of development where 
reliable or reasonably reliable predictions may be made concerning 
the speetra of energized electrons and ions. For example, while 
simple tearing analysis leads to a purely growlngmode, more detailed 
studies are likely to yield overstability instead (e.g., Drake and 
Lee, 1977) • Such a change may Influence the particle acceleration 
in a profound way. Thus, the observed layer of energetic electrons 
at the magnetopause (Meng and Andeisdn, 1970; Baker and Stone, 1977) 
may possibly provide important information concerning tearing, but 
at present we have no way of taking advantage of the Information. 


3. OGO-5 Grossing on March 10, 1968 at 02; 20 UT 

In order to Illustrate the ideas discussed in the previous sec- 
tion, we now examine in detail one OGO-5 crossing of the earth's 
magnetopause where we think tearing mode signatufes may have been 
p|re sent. Its gross magnetic structure is that of a rotational form 
so that a good determination of the normal direction is obtained in 
spite of the fluctuation in the normal field component, presumably 
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associated .with the tearing. The crossing also occurred at a sub- 

I""' ■" ■ ■ I ' 

stantial distance from the subsolar point: the solar magnetospheric 

latitude and longitude were 23.5 °N and 7^. ^°VT, respectively. This 
crossing has been studied previously by Sonne r up and Ledley (197^)* 
Information about the results of the minimum variance analysis may 


be. found in their paper. 


A polar (hodogram) representation of the magnetic field during 
the crossing is shown in Figure 2. In this figure the axes labeled 

B1 , B2 , and B3 are aligned with the directions of maximumj inter- 

^ r ; ; ' i : 

mediate, land minimum variance, respectively. The latter (B3) repre- 
sents the direction normal to the magnetopause. Thus , the field 
component B3 is identical with The right-hand hodogram indicates 

that an average magnet! e field component B3 = B « ~ 8nt is present, 
with a large amount of superimposed noise. The two field compbnents 
B1 and B2 are tangential to the magnetopause so that the left-hand 
hodogram shows the rotational behavior of the tangential field 


It = ^ - 

We observe three subsegments in this crossing, denoted by I, 

II and III, where a substantial amount of magnetic noise is present. 
By direct Inspection, it is evident that the noise is not isotropic: 


the amplitudes in the direction along the tangential hodogram trace 
and in the B3 direction are much larger than those in the direction 
perpendicular to the tangential hodogram trace. Since the electric 
current in a rotational discontinuity is directed essentially per- 
pendicular to the tangential hodogram trace (l.e., it is approxi- 
mately radial in the left-hand hodogram) , it is evident that the 
magnetic perturbation along the current direction is small for each 
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subsegment, a result that agrees with the constancy of B In the 
tearing mode model of Figure 1. Since the normal field component 
is different from zero, the tearing, if it is locally active, should 
be of the ion type. ^ But we cannot exclude the possibility that tear- 
ing "debris”, produced elsewhere, is being convected past the satellite. 
The ratio /B^„ is of the*order of 0.2 which may be in the unstable 

OX OZ 

i ' ‘ 

range (compare (ii) in the previous section). i , 

For each of the subsegments, I, II, III we have selected a local 

: j ; 

(y, 2 ) coordinate system such that the fluctuation in B is a 

y 

minimum. These local systems are shown in Figure 2. In each case, 
the third coordinate (x) remains along B3. 

The orientation of the convection velocity vector on the 
magnetosheath side of the layer is assumed to be tangent to the curve 
of Intersection of the magnetopause surface and the plane containing 


the earth’s center, the sun and the satellite. Within the magneto-; 
pause, structure, the convection velocity is presumably directed some- 

I L : — 

what differently and has different magnitude. If the magnetosheath 
plasma flows across the magnetopause, as one must assume for a 
rotational discontinuity, then the tangential velocity and magnetic 
field changes in the magnetopause are proportional, i.e., Av^ = ±AB^//ii^p 
where p is the mass density. The choice of sign depends on the sign 
of B3. For each of the subsegments in Figure 2, we have added a 
corresponding Av. to v (assuming v = 200 km/s and a density of : | 

U CG GO I 

10 protons/cm^) in order to obtain the local convection velocity v^. 

The orientations of v^ and the angle 0 between v^ and the local z 
axis (section 2c) are shown in Figure 2. It is noted that in a few 
recent cases the velocity change Av^. of the plasma has been directly 
measured. (Paschmann et al, 1979.) 




SI 


Table 1 gives a summary of results of detailed analysis of the 

three Bubsegments. A Past. Fourier Transform (PPT) was made of the 

data for each field compohent of the subsegment, giving: a discrete 

power spectrum whose points are multiple^ of the basic frequency 

fQi ''= (NT)“V, where t is i.he sampling interval and H the number of 

data points^ Spectra: of the .three components for subsegment I are 

shown in Figure 3. The large peak at f„ in B and smaller peaks ih 

and B correspond to the shear in the magnetopause field across 
^ y 

the layer. It is not of interest here. The other major peak, at 
.ll-il375 Hz in all components , cohtAins essentially all of the power 
in the field oscillations evident in the hodogram. 'No other signifi- 
cant peaks appear up to the Nyquist frequency fj^ = (2t)“* = 3*5 Hz. 

The; ratiop of areas under, the major peaks in the three components 
give the power ratios among the field components (see Table 1) . For 

j ’ ^ ^ 

subsegment I the corresponding amplitude ratios are [B, I/]B„|/iB ] = 10/ 
^.6/1.^. The angle 0 for segment I is about ^0®, and the local con- 
vection speed is v_ = 275 km/s. Thus the inferred wavelength is 


mi 


lit! 


I i 

i 


X = V cos 
c 


observed 


- (275 km/s) (cos40®)/0 . 125 Hz 1685 km 


When the observed amplitude ratio |B^|/|B^| - 0 .kS is substituted 
into Eq. (5) there results 


0.46/k = 0.jt6X/2Tr = 125 km 


The interpretation of this magnetopause substructure in tc:: s of a 
. ^ ; ■ i ' ' ■ ^ " ■ 

tearing mode thus gives a realistic estimate for the scale length 

of the magnetic field shear in this portion of the layer. Similar 
calculations for subsegments II (0 - 65®, v_ = ^00 km/s) and III 








- 15 - 

(0 e 8o°, Vg ~ H55 km/s) yield i values of 60 km and 35-65 km, re- 
spectively (see Table 1) . 

None of the subsegments has significant phase coherence in the 
plane so that no conclusion can be drawn from detailed comparison 
with the tearing-mode polarization structure described earlier in 
the paper. • 

Pioneer 11 Crossing on Nov 29, ig'?^ at 13:30 UT 

i This crossing of Jupiter’s magnetopause is the last one on the 
inbound leg of the Pioneer 11 trajectory. This crossing too occurred 
a substantial distance av;ay from the subsolar point: the solar-jovian 

latitude and longitude were 9.5°S and respectively, and R = 61|.7R* 

The hodogram of the crossing is shown in Figure ^ . Except for 
the field magnitude, which is much lower in the jovian case, the 
overall features of this crossing are similar to the terrestrial 
crossing in Figure 2. The Jupiter magnetopause is a clear case of 
a rotationai form with a nonzero normal field component and it con- 
tains several subsegments of magnetic fluctuations similar in char- 
acter to those in Figure 2. 

Table 1 gives results of the analysis of the subsegment I in 
Figure These; results show the presence of a mode of the same 
type as above, with a somewhat lovirer frequency, and a smaller ampli- 
tude, but similar ratio of amplitude to total field strength. In 
this subsegment there is also definite phase coherence be tvM'een B : 

A ; 

and 'as shown in Figure 5* B^ generally leads B by ^ during the 

■ . Z , A Z c. 

'V. , , : I : . ,, ■ ■■■■■ ^ ■' 

whole subsegment, although clearly there are significant fluctuations 

in this relation, and especially complex behavior near the center 
of the subsegment. Although the time plot of Figure 5 might Suggest 
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that the observer had passed through x - 0 near the middle of the 
•subsegment, there is not an overall polarization antisymmetry, and 
the frequency is actually slightly higher in the second half of the 
subsegment. It is therefore probable that two separate structures, 
rather than the two sides- of a single mode, are crossed here . This 
segment is nevertheless the clearest example we have yet found of 
the type of magnetic oscillation we are describing here. 

For the segment I the local convection velocity v^, directed as 

i ■ • - ; 

shown in Figure has a magnitude of 190 km/s and forms an angle 
0 s; 58 ° with the s axis of the local coordinate system. In obtain- 
ing these results we used v^ = 200 km/s and a density of 2 protons/cm^. 

-C 

Performing the same: calculation as before we then find the wavelength 
X « l680 km and the shear length i - I 70 km. These numbers appear 
somewhat too small when compared, to an estimated thickness of Jupiter’s 
magnetopause of ~^{Q0O kni (Tsurutanl et al, 1979) / However, it should 
be remembered that a is associated with the width of the substructures, 
and not with the total magnetopause thickness. 

.other segments of this crossing also show magnetic oscillations 
of the same character, in the frequency range .03-. 1 Hz, but. have 
several superposed frequencies, or are very short intervals, so they 
are not examined in detail . 

5 . Pis cusslon 

j An important common feature of the four data segments we have 
analyzed is that the z axis of the local coordinate system, chosen 
for minimum variation in B , is always approximately perpendicular 
to the local tangential field direction at the time of observation 
of the magnetic oscillations . This means that, although we have not 
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found evidence of actually crossing the k • B = 0, i.e. the x = 0, 
plane of any of the structures, the oscillations, if produced by 
tearing, are all observed fairly close to that plane. It is a basic 
feature of the tearing perturbations that they decay exponentially 
with increasing distance away from x = 0. Thus they should be easily 
observable only relatively near x =0. 

Examination of five other 0G0^5 crossings has shown morb than 
twenty subsegments of large amplitude oscillations similar to those 
described here. The local coordinate system for each subsegment is 
always such that the z axis is nearly perpendicular to the local field 
Although we have not been able ^t show clear evidence of the 
polarization behavior (Fig. 1) associated with tearing-mode magnetic 
islands it appears that the observed oscillations are consistent 
with tearing Irt the following ways: 

i) The magnetic oscillations are very nearly two-dimensional with 
virtually no oscillations in the B component. 

• y 

ii) The observations are Comparable with the hypothesis that the 

i ■■ 

oscillations are confined to a relatively narrow region about 
the k • B = 0 plane . 

iii) The oscillations fall in a frequency range ( 0.0 5-0. Hz) which, 
with typical magnetopause convection velocities , gives reason- 
able inferred values of the wavelength. 

iv) Tlie ratio of the amplitudes of the two oscillating field compo- 
nents is consistent with reasonable values of the shear length 
which is representative of the thickness of the tearing substruc- 
tures . 

Two alternative interpretations should be mentioned. First, 

oscillatbry radial motion of the magne topause would produce B os cll- 

z 
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lations in a local coordinate system with its z axis perpendicular 
to the iocal tangential field. But it would not generate an associated 
oscillation in the norinal component B . This latter objection may 

, , A. 

be removed in a second model in which the radial oscillation is the 

^ : i i I , • ; , ' 

I : ■ ■ : . :-■'■■■ . ^ ' 

result of small amplitude wave motipn of the entire magnetopause, 

j ' ' ' ^ 

with the wave vector Ic having a substantial component along the local 

tangential field in the magnetopause. Without actual observations 

of the polarization reVersal at x = 0 in the tearing mode it is not 

eajsy to discriminate between tearing structures and this type of 

wave motion. However, if the wave; propagation direction k remains 

fixed (along v say)', the wave model predicts that the amplitude 

of, the B„ oscillation should be a maximum for subsegments where the 

i ^ 

local magnetic field is parallel to k. No such effect is seen either . 
in Figure 2 or Thus we feel that the tearing mode interpretation 
is stronger than the wave interpretation for the cases analyzed here. 
However, study of a much larger number of events is required in ordor 
to discriminate in a reliable way between the two possibilities. - 
It also remains possible that some other instability or wave 
motion might be responsible for these magnetic oscillations . Progress 
must be made in the theory of both tearing and other instabilities 
alnd wave modes in strongly sheared, thin collision-free current layers 
with plasma flow, as well as in the observations of su structures 
at planetary magnetopauses, before definite conclusions may be drawn.; 
The most Important observational improvement will be simultaneous 
measurements from separated points in space, which the ISEE mission 
is now providing. Independent, and unambiguous information about the 
magnetopause thickness and other basic parameters will be important 
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in deciding among possible modes. In particular, the growth of 
tearing modes is controlled .by the layer thickness (the shear length 
fi.) and by the size of the normal component. Since it is somewhat 
unlikely that sufflGiently clear signatures of phase and polarization 
will be found to permit of an unambiguous Identification of the mode, 
correlation of its occurrence^ with magnetopause parameters may be 
the strongest evidence that can be added toward identification. 

In summary, we have shown evidence of the existence of a partic- 
ular two-dimensional low-frequency oscillation of the magnetic field 
in the dayside magnetopause current layer of the earth and of Jupiter. 
This oscillation has been shown to be consistent, in its location, 
orientation, amplitude and frequency with our present understanding 
of the tearing mode. However, the tearing mode interpretation is not 
necessarily unique . • 
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Table-1-.-— J^req^iency- and-waveiength^^-^^^ three 060-5 data 

segments and one Pioneer 11 segment. 


Segment 

Duration 
( sec) 

fo=(NT)- 

Hz 

Spectral 

peaks 

(ffo)Hz 

■ B‘ " B 

X y 


Peak 

amplitude 
. _|BzInt 

Power” in 
spectral 
peaks 

Inferred 

X 

(km) 

Inferred 

Z 

• (km) 

060-5 

36 .6 

0.0275 

0.11- 

0.1375 

0.11 

10-15 

10:2.1:0.2 
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0.22 


10 

10: 1.5: <0.1 
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65 

H i 
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O.ll 

0.11 

0.11 

5 

10:3:1.1 
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65 

III 
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- i 

0.165 

0.165 

' ]■ 


10:2:<0.1 

il80 

35 

Pioneer 11 

j. 

96 

0.01 

0 .06 ■ 

0.06 

0.06 

1.5-2 : 

10:^1 :0.^ 

1680 

170 
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Flgure Captions 

Pig. 1. Schematic representation of tearing mode perturbation. 

Central x = 0 (k‘B=0) planb is shown by dotted line. 

Diagonal dotted line represents path of a spacecraft 
crossing the layer as the structure is convected past it. 
Plots of Bg vs Bx indicate polarizations of oscillations 
which would be observed at the fixed x positions indicated. 

A is the vector potential. 

• . ■ ; 

Pig. 2. Polar (hodogram) representation of OGO-5 magnetopause 
crossing. The left hand diagram shows the behavior of 
t>he magnetic field components B1 and B2 bangentlal to the 
magnetopause; the right-hand diagram shows the normal 
component B3.| The axes are the principal axes obtained 
from minimum variance analysis and shown in GSM coordi- 
nates in lower left-hand corber. Three subsegments^ I> : 

. II, and III, of field oscillations are present. For each 
of these local x, y, z Coordinate system, corresponding 
to the one in. Figure 1 is identified. Also shown is the 
local convection velocity vector and angle 9 for each sub- 
segment . 

Pig. 3. Past Fourier transform power spectra of the three components 
for OGO-5 segment I. Vertical scale is arbitrary, horizontal 
scales are in units of .11 Hz. All spectra continue to 
■ decrease beyond the frequency range shown, up to the.Nyqulst 
frequency at 3 • 5 Hz. 

Fig. ^ . Hodogram representation of Pioneer 11 crossing of Jupiter ' s 
magnetopause . See Figure 2 for details . The principal 
axes are given in the solar jovlan (SJ) system. 

Fig. 5. B^ (solid line) and B„ (dashed line ) oscillations for a 

150 sec. interval of uhe crossing of Jupiter’s magnetopause. 
Polar plots are for the three subintervals indicated . 
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